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Disclaimer 

This  report  is  intended  as  a record  of  the  presentations  and  discussions  that  took  place  at  a 
NIST  Metallurgy  Division  and  NIST  Center  for  Theoretical  and  Computational  Materials 
Science  sponsored  workshop.  The  opinions,  conclusions,  or  recommendations  that  are  expressed 
herein  are  those  of  the  organizers  or  individual  presenters  and  do  not  necessarily  reflect  the 
views  of  NIST.  All  references  to  commercial  products  in  this  report  are  for  identification 
purposes  only  and  do  not  constitute  any  endorsement  by  NIST.  The  policy  of  NIST  is  to  use  the 
International  System  of  Units  (SI,  metric  units)  in  all  its  publications.  In  this  document  however, 
works  of  authors  outside  NIST  are  sited  which  describe  properties  in  certain  non-SI  units. 
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Abstract 


Databases  for  computational  thermodynamics  and  diffusion  modeling  can  be  applied  to 
predict  phase  diagrams  for  alloy  design  and  alloy  behavior  during  processing  and  service. 
Databases  that  are  currently  available  to  scientists,  engineers  and  students  need  to  be  expanded 
and  improved.  The  approach  of  the  workshop  was  to  first  identify  the  database  and  information 
delivery  tool  needs  of  industry  and  education.  Improved  modeling  capabilities  result  in 
decreased  new  alloy  and  process  development  time  and  ultimately  in  cost  savings.  A roadmap 
will  be  developed  for  how  these  needs  can  be  met  during  the  next  decade  in  a cost  effective  way 
through  expanded  collaborative  efforts  in  education,  basic  research  and  database  development. 

The  workshop  format  was  a series  of  invited  talks  given  to  the  group  as  a whole  followed 
by  general  discussions  of  needs  and  benefits  to  provide  a roadmap  of  future  activities. 
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Program 


Thursday,  March  21 

8:00  Opening  Remarks 

Ursula  R.  Kattner,  NIST,  Gaithersburg,  MD 

John  E.  Morral,  University  of  Connecticut,  Storrs,  CT 

8:10  Air  Force  Programs  in  Materials  Engineering  and  Design 

Craig  S.  Hartley,  Air  Force  Office  of  Scientific  Research  (AFOSR),  Arlington,  VA 

8:40  Application  and  Needs  Assessment  of  Thermodynamic  Databases  at  Howmet  Castings 
Ty  W.  Hansen  , Howmet  Research  Corporation,  Whitehall,  MI 

9: 10  Computer  Simulations  and  Database  Needs  for  Energy  Savings,  Zero  Environmental 

Impact  and  Quality  Improvements  For  the  Heat  Treating  Industry 
Rick  Sisson,  Worcester  Polytechnic  Institute  (WPI),  Worcester,  MA 

9:40  Materials  Programs  at  NSF 

K.L.  Murty,  National  Science  Foundation  (NSF),  Arlington,  VA 

9:50  Break 

10:00  Microstructure  - Process  - Property  Models  for  the  Aluminum  Industry 
Joanne  L.  Murray,  Alcoa  Technical  Center,  Alcoa  Center,  PA 

10:30  Applications  of  Computational  Thermodynamics  to  Virtual  Aluminum  Castings 
Ravi  Vijayaraghavan,  Ford  Motor  Company,  Dearborn,  MI 

1 1 :00  Use  of  Quasi-Thermodynamical  Models  for  Simulation  of  Epitaxy  of  III-V  Compound 
Semiconductors  in  Electronic  Industry 

Yuri  Makarov,  Semiconductor  Technology  Research  Inc.,  Richmond,  VA 

11:30  Lunch 

12:30  NSF  Supported  Computational  Education  Program  at  Penn  State 

Zi-Kui  Liu,  Penn  State  University,  University  Park,  PA 

1 :00  Thermodynamic  Measurements 

Philip  Nash,  Illinois  Institute  of  Technology,  Chicago,  IL 

1 :30  Configurational  Entropies  in  Real  Alloy  Phase  Diagram  Calculations 
Fan  Zhang,  CompuTherm  LLC,  Madison,  WT 
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2:00  Current  Capabilities  in  First-Principles  Modeling  of  Alloy  Thermodynamics 
Patrice  E.A.  Turchi,  LLNL,  Livermore,  CA 

2:30  Break 

2:45  Database  Development  and  Industrial  Applications 

Arthur  D.  Pelton,  Centre  de  Recherche  en  Calcul  Thermochimique  (CRTC),  Montreal, 
Quebec 

3:15  Development  of  Thermodynamic  Databases 

Philip  J.  Spencer,  The  Spencer  Group,  Ithaca,  NY 

3:45  Application  of  Thermodynamics  to  Alloy  Design  and  Development 
Charles  J.  Kuehmann,  QuesTek  Innovations  LLC,  Evanston,  IL 

4:15  Construction  and  Application  of  a Diffusion  Mobility  Database  for  Ni-Base 
Superalloys 

Carelyn  E.  Campbell,  NIST,  Gaithersburg,  MD 

4:45  Incorporation  of  CALPHAD  Calculations  into  Phase  Field  Modeling 

William  J.  Boettinger,  NIST,  Gaithersburg,  MD 

6:00  Dinner 

7:00  Current  and  Future  Applications  of  CALPHAD  Technology 

Larry  Kaufman,  Brookline,  MA 


Friday,  March  22 

8:00  Thoughts  on  the  Acceleration  of  Phase  Equilibria  Research 

J.-C.  Zhao,  General  Electric  Company,  Schenectady,  NY 

8:30  Discussion  and  Roadmap  Oversight 

10:30  Summary  and  Conclusion 

11:30  Lunch 


12:30  Adjourn 
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Presentation  Summary 

The  workshop  program  consisted  of  three  parts.  The  program  on  Thursday  morning  was 
devoted  to  the  identification  of  national  needs  and  benefits  featuring  speakers  from  industry, 
funding  agencies  and  academia.  The  program  on  Thursday  afternoon  featured  speakers  from 
academia,  national  laboratories  and  consulting  firms  presenting  current  capabilities  and 
software.  The  program  on  Friday  morning  was  dedicated  to  roadmap  plan  development. 

The  workshop  began  with  opening  remarks  by  J.E.  Morral  (University  of  Connecticut) 
and  U.R.  Kattner  (NIST  Metallurgy  Division)  in  which  the  objectives  and  the  planned  outcome 
for  the  workshop  were  presented.  The  first  two  presentations  focused  on  the  data  need  of  the 
aerospace  industry.  C.S.  Hartley  (Air  Force  Office  of  Scientific  Research  (AFOSR)) 
emphasized  the  immediate  need  for  improved  performance  of  engineered  materials.  He 
presented  the  MEANS  (Materials  Engineering  for  Affordable  New  Systems)  program  at  AFOSR 
that  fosters  the  development  of  scientific  principles  and  understanding  which  are  the  basis  for 
computational  materials  engineering,  data  needs  and  necessary  experimental  validation. 

T.W.  Hansen  (Howmet  Research  Corporation)  described  the  products  and  product  requirements 
at  Howmet.  He  discussed  the  application  of  current  phase  equilibria  modeling  capabilities  and 
their  limitations.  He  then  presented  future  data  and  software  needs  to  improve  and  expand  the 
modeling  efforts.  The  presentation  of  R.D.  Sisson  (Worcester  Polytechnic  Institute  (WPI), 
Center  for  Heat  Treating  Excellence)  showed  how  computer  simulations  can  be  employed  for 
energy  savings,  environmental  and  quality  improvements  for  the  heat  treating  industry.  He 
pointed  out  that  heat-treating  modeling  is  limited  by  the  available  databases,  i.e.,  look-up  tables 
and  that  computational  thermodynamics  could  solve  part  of  the  data  problem. 

K.L.  Murty  (National  Science  Foundation  (NSF))  gave  an  overview  of  funding  provided 
by  the  Division  of  Materials,  Metal  Program.  He  stated  that  the  majority  of  the  2001/2002 
funding  increase  went  to  bio-,  nano-  and  information  technology  programs. 

The  next  two  presentations  focused  on  data  needs  for  aluminum  and  other  light  metal 
alloys.  J.L.  Murray  (Alcoa)  gave  an  overview  of  commercial  aluminum  alloys  and  properties 
which  are  important  for  the  manufacturing  process.  She  pointed  out  that,  although  many  of  the 
binary  phase  diagrams  with  aluminum  are  known,  not  all  phase  boundaries  are  known  with 
sufficient  accuracy  and  identified  the  systems  that  need  further  refinement.  R.  Vijayaraghavan 
(Ford  Motor  Company)  presented  the  modeling  efforts  in  the  design  and  processing  of  cast 
aluminum-alloy  engine  parts.  For  this  modeling,  a suite  of  software  codes  is  employed.  The 
codes  include  property  calculations,  as  well  as  structural  analysis.  The  data  needs  for  the 
modeling  tools  were  highlighted. 

In  last  talk  of  the  Thursday  morning  session,  Yu.N.  Makarov  (Semiconductor 
Technology  Research,  Inc.)  talked  about  the  modeling  requirements  for  the  manufacturing  of 
IH-V  compound  semiconductors.  He  presented  quasi-thermodynamical  models  that  also  take 
into  account  kinetic  processes  during  the  epitaxial  growth  of  the  semiconductors. 

The  presentations  of  the  Thursday  morning  session  clearly  demonstrated  the  power  of 
modeling  approaches,  but  also  underscored  the  urgent  need  for  reliable  databases,  such  as 
thermodynamic,  kinetic  and  physical  property  data.  Although  the  construction  of  databases  has 
made  significant  progress  in  recent  years,  faster  progress  is  needed. 
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In  the  afternoon  the  focus  was  on  current  capabilities  and  programs  in  data  generation 
and  assessment.  Z.-K.  Liu  (Penn  State  University)  gave  an  overview  of  the  redesigned 
Materials  Science  and  Engineering  program  at  Penn  State  University.  This  program  integrates 
traditional  education  in  thermodynamics  and  kinetics  of  materials  with  new  courses  that  employ 
computational  models.  The  goal  was  not  only  to  include  modern  material  development  tools  into 
the  curriculum,  but  also  to  stimulate  student  interest. 

Although  thermodynamic  modeling  can  reduce  the  amount  of  experimental  data  needed, 
measured  data  are  crucial  for  the  development  of  databases  and  verification  of  the  predictions 
obtained  by  these  databases.  P.  Nash  (Illinois  Institute  of  Technology  (DT))  presented  an 
overview  of  the  experimental  methods  that  are  available  for  thermodynamic  measurements  at 
IIT.  At  the  end  of  the  presentation  an  overview  of  centers  in  the  U.S.A.  and  worldwide  carrying 
out  experimental  research  activities  in  this  field  was  given.  There  are  so  few  such  centers  that  a 
critical  mass  barely  exists  to  ensure  continued  competence  in  the  U.S.A.  in  the  measurement  of 
thermodynamic  data  in  the  future. 

The  following  two  presentations  focussed  on  the  contribution  that  ab  initio  modeling  can 
make  to  the  modeling  of  real  materials.  F.  Zhang  (CompuTherm,  LLC)  showed  how  results 
from  first  principles  ( ab  initio ) calculations  can  be  used  to  improve  the  model  descriptions  that 
are  used  in  phenomenological  modeling  approaches,  such  as  the  Calphad  method. 

P.E.A.  Turchi  (Lawrence  Livermore  National  Laboratory  (LLNL))  gave  an  overview  of  the 
available  ab  initio  methods,  software  packages  and  property  data  that  can  be  obtained  from  these 
calculations.  (See  also  “Some  Notes  on  First-Principles  Electronic  Structure  Methods  and 
Calculations,”  page  17.)  He  also  gave  examples  on  how  the  results  from  ab  initio  calculations 
can  provide  data  for  other  computational  modeling  approaches,  i.e.,  the  Calphad  method. 

These  presentations  were  followed  by  two  talks  on  database  development.  A.D.  Pelton 
(Centre  de  Recherche  en  Calcul  Thermochimique  (CRTC))  presented  an  overview  of  the  steps  in 
developing  databases  and  a description  of  the  decision-making  process  in  selecting  models.  He 
presented  the  databases  that  were  obtained  from  this  strategy  and  gave  examples  of  their 
application  to  industrial  processes.  P.J.  Spencer  (The  Spencer  Group)  gave  an  overview  of 
available  thermodynamic  databases  for  pure  inorganic  substances  and  alloy  systems.  He  pointed 
out  the  needs  and  benefits  of  database  construction  and  emphasized  that  self-consistent  databases 
require  careful  construction. 

The  remaining  presentations  of  Thursday  afternoon  were  dedicated  to  the  application  and 
coupling  of  thermodynamic  calculations  with  kinetic  modeling.  C.J.  Kuehmann  (QuesTek 
Innovations,  LLC)  gave  examples  for  the  incorporation  of  thermodynamic  calculations  into  a 
larger  scheme  of  materials  design.  Different  software  suites  are  used  depending  on  the 
dimensionality  of  the  material  (atomic,  nano,  micro,  ...).  Results  from  the  software  package 
PrecipiCalc  were  demonstrated.  C.E.  Campbell  (NIST  Metallurgy  Division)  reported  the  results 
of  the  construction  of  a mulicomponent  diffusion  database  for  superalloys.  This  database  was 
constructed  using  the  same  principles  as  for  CALPHAD  type  thermodynamic  databases  and  it 
was  shown  that  it  has  the  same  predictive  power.  W.J.  Boettinger  (NIST  Metallurgy  Division) 
discussed  a few  approaches  to  phase  field  modeling.  Examples  were  shown  for  solidification 
and  precipitation  simulations.  He  showed  how  thermodynamic  and  mobility  data  fit  naturally 
into  the  phase  field  modeling  approach.  Although  phase  field  modeling  is  a powerful  tool  for  the 
simulation  of  growth  processes  its  application  is  still  limited  by  computational  power. 
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The  final  presentation  of  the  first  day  was  given  by  L.  Kaufman.  He  presented  a series 
of  examples  of  analysis  of  multicomponent  materials  to  illustrate  the  wide  variety  of  applications 
for  computational  thermodynamics  and  diffusion  modeling. 

The  discussions  on  Friday  morning  began  with  a short  talk  by  J.-C.  Zhao  (General 
Electric  Company)  in  which  he  presented  his  thoughts  on  the  acceleration  of  phase  equilibria 
research  using  concepts  of  high  throughput  or  combinatorial  methods.  He  discussed  the  use  of 
samples  constructed  of  diffusion  multiples  to  accelerate  experimental  phase  diagram 
investigations.  He  also  showed  how  the  results  were  used  to  verify  the  databases.  The 
remaining  time  of  Friday  morning  was  dedicated  to  the  discussion  of  what  was  learned  from  the 
presentations  and  to  the  development  of  a draft  of  a roadmap.  The  discussion  was  lead  by  W.J. 
Boettinger  (NIST  Metallurgy  Division).  The  results  of  the  discussion  and  roadmap  oversight 
are  summarized  in  the  following  chapter. 
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Computational  Thermodynamics  and  Diffusion  Modeling  Roadmap 

The  roadmap  is  the  result  of  the  discussions  following  the  presentations  at  the  workshop  and  has 
been  agreed  on  by  the  majority  of  the  workshop  participants. 


Vision  2010 

By  the  year  2010  all  Materials  Science  and  Engineering  students  will  learn  to  solve  real 
world  materials  problems  using  thermodynamic  and  diffusion  modeling  software,  while 
materials  engineers  in  industry  will  have  the  computational  tools  and  the  background  needed  to 
lead  the  world  in  the  design  and  optimization  of  both  materials  and  materials  processing. 


Benefits 

The  benefits  of  providing  materials  engineers  with  the  knowledge,  software  and 
databases  to  model  materials  and  processing  are  savings  in  costs,  energy  and  natural  resources 
via  optimization  studies.  Also,  properties  can  be  improved,  new  materials  can  be  discovered  and 
development  time  can  be  reduced.  Specific  benefits  are: 

• Computational  thermodynamics  has  already  been  successfully  used  in  the  development  of 
new  alloys  and  the  understanding  of  microstructures. 

• Reduced  product  cycle  time  via  Federal  programs,  such  as  MEANS  (Materials  Engineering 
for  Affordable  New  Systems)  and  AIM  (Accelerated  Insertion  of  Materials),  will  reduce 
costs  sharply  and  make  USA  industry  more  competitive;  e.g.,  savings  to  the  auto  industry  on 
engine  blocks  and  cylinder  heads  could  top  $ 100  M. 

ic  i '■y 

• Reduction  in  energy  use  by  the  heat  treating  industry  alone  is  estimated  as  7x10  ' J (7x10 
BTU)  over  a ten-year  period  due  to  improved  efficiency. 


Approach 

• Organize  a series  of  “work  groups”  consisting  of  representatives  from  Academia,  Industry 
and  Federal  Agencies  that  will  oversee  and  expand  plans  to  create  new  teaching  tools  for 
Universities  and  Technical  Colleges. 

• Build  an  experimental  infrastructure  that  can  make  ultra-high  precision  thermodynamic  and 
diffusion  measurements. 

• Build  an  educational  infrastructure  for  the  teaching  and  application  of  computational 
thermodynamics  and  diffusion  modeling. 

• Prepare  a comprehensive  public  library  of  thermodynamic  and  diffusion  property  databases 
that  apply  to  common  and  emerging  alloy  systems. 

• Increase  processing-structure-property  modeling  efforts,  i.e.,  improve  application  modeling. 
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Objectives,  Strategies  and  Action  Items 

Coordinators  are  suggested  for  each  of  the  action  items. 


I.  Education 

1.  Interest  faculty  in  teaching  applied  thermodynamics  and  kinetics  with: 

a.  New  course  materials  (traditional  and  web  based). 

b.  Industrial  contacts. 

c.  Research  opportunities  (e.g.,  via  joint  projects). 

2.  Interest  students  in  applied  thermodynamics  and  kinetics  with: 

a.  Real  world  problems. 

b.  Hands  on  experience  doing  computational  thermodynamics  and  diffusion 
modeling. 

c.  Job  opportunities  in  computational  materials  design. 

3.  Increase  the  technical  ability  of  students  with  more  emphasis  on: 

a.  Better  appreciation  of  the  underlying  physics  as  well  as  classical  concepts. 

b.  The  understanding  of  the  behavior  of  multicomponent  alloys,  not  just  binary  alloys. 

4.  Encourage  institutions  to  support  teaching  initiatives: 

a.  Provide  funding. 

b.  Support  outside  funding  initiatives. 

5.  Provide  training  for  those  already  in  industry. 

Action  Items 

1 . Prepare  a collection  of  useful  applied  problems. 

- Karl  Spear  (Penn  State) 

2.  Publish  an  article  about  the  NSF  program  at  Penn  State  on  computational 
thermodynamics  and  kinetics. 

- Zi-Kiu  Liu  (Penn  State) 

3.  Prepare  a catalog  of  currently  available  course  materials. 

- TBD 

4.  Explore  the  possibility  of  a course  on  computational  methods. 

- Bill  Scott  (ASM) 


II.  Thermodynamic  and  Diffusion  Data 

1 . Experimental  data: 

a.  Promote  a national  initiative  to  rebuild  the  thermodynamic  and  diffusivity 
measurement  infrastructure. 
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b.  Establish  goals  for  improving  the  accuracy  of  phase  diagram  measurements  by  an 
order  of  magnitude. 

c.  Establish  “best  practice  guide”  for  using  assessed  data. 

d.  Make  assessment  software  more  user  friendly,  more  robust,  and  with  less 
complexity  (i.e.,  make  it  run  faster). 

2.  Virtual  data: 

a.  Encourage  First  Principles  calculations  of  interaction  energies  and  defect  formation 
energies. 

b.  Encourage  Molecular  Dynamics  calculations  of  atomic  mobilities. 

c.  Develop  estimation  software  for  unknown  thermodynamic  parameters  (empirical, 
semi-empirical,  first  principles),  for  example  Miedema's  model  for  the  estimation 
of  enthalpies  of  formation. 

d.  Develop  estimation  methods/software  for  atomic  mobilities/diffusion  data. 

Action  Items 

1 . Review  the  status  of  estimating  molar  volume. 

- Phil  Nash  (HT). 

2.  Expand  the  scope  of  the  CALPHAD  community  to  include  diffusion  modeling. 

- Bill  Boettinger  (NIST),  Zi-Kiu  Liu  (Penn  State). 

3.  Capture  the  diffusion  data  in  NIST  archives. 

- TBD 

4.  Organize  diffusion  assessments  and  modeling  sessions  at  the  2003  CALPHAD  meeting. 

- John  Morral  (UConn),  Afina  Lupulescu  (RPI),  Arthur  Pelton  (CRTC) 

5.  Organize  diffusion  modeling  and  assessment  sessions  at  ASM  meetings  through  the 
Atomic  Transport  Committee. 

- Rick  Sisson  (WPI),  John  Morral  (UConn). 

6.  Prepare  a “wish  list”  of  needed  thermodynamic  and  kinetic  information  for  physicists  to 
model. 

- TBD 

7.  Prepare  a list  of  physicists  doing  First  Principles  and  Molecular  Dynamics  studies  and 
distribute  the  “wish  list”  to  them. 

- Patrice  Turchi  (LLNL) 

8.  Develop  standard  problems  for  establishing  “best  practices”. 

- Ray  Thompson  (UAB),  Charlie  Kuehmann  (QuesTek). 


III.  Public  Databases 

1 . Encourage  the  publishing  of  a full  set  of  parameters  used  in  phase  diagram  and 
diffusivity  assessments  that  appear  in  the  open  literature. 

2.  Develop  public  databases  that  are  comprehensive  for  common  commercial  alloys 
systems  and  emerging  alloy  systems. 

3.  Develop  databases  for  other  data:  molar  volume,  viscosity,  etc. 

4.  Encourage  a standard  format  be  used  in  public  databases. 
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5.  Develop  methods  of  delivering  public  databases  into  macrocodes;  encourage 
standardized  interfaces. 

Action  Items 

1 . Write  letters  to  journal  editors  encouraging  policies  that  require  authors  to  list 
parameters  used  in  assessment  articles. 

- Everyone 

2.  Encourage  authors  of  proprietary  databases  to  publish  descriptions  of  constituent 
subsystems. 

- Everyone 

3.  Recommend  and  publish  the  SGTE  format  for  data  (CALPHAD  web  site  and  journal). 

- Zi-Kiu  Liu  (Penn  State). 

4.  Encourage  North  American  participation/membership  in  SGTE,  which  is  now  open  to 
non-European  members. 

- Phil  Spencer  (The  Spencer  Group). 

5.  Compile  database  and  software  needs  for  practical  applications. 

- Charlie  Kuehmann  (QuesTek). 

6.  Devote  an  issue  of  the  CALPHAD  Journal  to  data  delivery,  coupling  of  micro  and 
macro  computer  codes. 

- Zi-Kiu  Liu  (Penn  State). 

7.  Encourage  links  to  CALPHAD  web  site  (http://www.calphad.org). 

- Everyone 


IV.  Funding 

1.  Make  use  of  NSF  undergraduate  initiatives  to  obtain  funding  for  education  programs. 

2.  Obtain  public  funding  for  experimental  measurement  of  data. 

3.  Promote  the  value  of  thermodynamic  and  diffusion  modeling. 

Action  Items 

1.  Distribute  present  roadmap  to  contract  monitors,  industrial  research  laboratories, 
funding  agencies. 

- Everyone 

2.  Create  a list  of  materials  that  have  been  certified  to  meet  specified  criteria  with  the  aid 
of  computational  thermodynamics  and  diffusion  modeling. 

- TDB 

3.  Obtain  a grant  to  study  the  effect  of  thermodynamic  and  diffusion  modeling  on  the 
economy. 

- TDB 

4.  Approach  researchers  in  industry  to  supply  data  on  financial  benefits  of  using 
computational  thermodynamics. 

- Everyone 


Workshop  Report 


11 


V.  Work  Groups 

1.  Educational  Programs  and  Instructional  Materials. 

- Zi-Kiu  Liu  (Penn  State). 

2.  Experimental  Measurement  of  Thermodynamic  Data. 

- Phil  Nash  (HT) 

3.  Diffusivity  Assessments  and  Modeling. 

- Bill  Boettinger  (NIST),  Carrie  Campbell  (NIST),  John  Morral  (UConn) 

4.  Thermodynamic  Assessments. 

- Ursula  Kattner  (NIST) 
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The  following  section  was  contributed  for  the  report  after  the  workshop. 


Some  Notes  on  First-Principles  Electronic  Structure  Methods  and  Calculations 

by  Patrice  E.  A.  Turchi 

LLNL  (L-353),  P.O.  Box  808,  Livermore  CA  94551 
E-mail  Turchil@Unl.gov 


I.  When  performing  so-called  first-principles  (or  ab  initio ) electronic  structure  calculations,  it  is 

always  wise  to  keep  in  mind  the  following  constraints: 

A.  Constraints  within  Density  Functional  Theory  (DFT) 

1.  Local  Density  Approximation  (LDA)  (or  LSDA:  Local  Spin  Density  Approximation) 

Exchange  Potential,  e.g. : 

- von  Barth  & Hedin 

- Vosko 

- Ceperley  & Alder  (Perdew  & Zunger) 

The  selection  of  an  exchange  potential  will  impact  the  results,  e.g.,  the  total  energy. 

2.  Beyond  LDA 

LDA+U 

LDA++ 

GGA  (Generalized  Gradient  Approximation) 

SIC  (Self-Interaction  Correction) 

GWM  (GutzwiUer  Wave  Method) 

DMFT  (Dynamical  Mean-Field  Theory) 

Going  beyond  the  LDA  may  be  necessary  in  some  instances  when  electron  correlation  play  an 
important  role,  e.g.,  to  reproduce  the  bandgap  in  semiconductors,  the  electronic  behavior  of 
oxides  such  as  NiO  or  MnO  or  of  strongly  correlated  systems  such  as  those  based  on  Ce  and  Pu, 
or  when  magnetism  is  involved. 

B.  Other  Constraints 

1.  Shape  of  the  Potential  Function: 

- Muffin-Tin 

- ASA  (Atomic  Sphere  Approximation;  with  equal  or  non  equal  sphere  radii) 

- FuU  Potential 
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2.  Core  versus  Valence  Electrons  (This  selection  is  usually  not  talked  much  about  but  may 
be  important  if  calculations  are  carried  out  as  functions  of  pressure). 

3.  Scalar  Relativistic  versus  Fully  Relativistic  (For  magnetic  systems  and  heavy  elements, 
including  the  5d  transition  metals  and  beyond,  a fully  relativistic  treatment  of  the  electrons  is 
recommended). 

4.  Magnetism:  Co-linear  (Ising-like),  Non  Co-linear  (Heisenberg-like)  (In  the  case  of 
magnetic  systems,  it  may  be  necessary  in  some  instances  to  carry  out  large  supercell  calculations 
to  predict  the  proper  type  of  magnetism,  such  as  in  the  case  of  chromium). 

5.  Atomic  Positions  and  Site  Occupancy  (A  calculation  for  an  ordered  compound  requires 
the  knowledge  of  the  atomic  positions;  in  the  case  of  an  off-stoichiometric  or  chemically  random 
alloys,  additional  approximations  are  necessary). 

6.  Bom-Oppenheimer  approximation  (electron-phonon  coupling)  (This  approximation 
may  be  severe  when  performing  molecular  dynamics-type  calculations). 


II.  Two  main  categories  of  methodologies  are  defined  depending  on  the  desired  type  of 

calculations: 


A.  Methods  Limited  to  Ordered  Structures  (Wave  Methods) 

- Pseudo-Potential 

- LAPW  or  FP-LAPW  (Full  Potential  Linear  Augmented  Plane-Wave) 

- LMTO  or  FP-LMTO  (Full  Potential  Linear  Muffin-Tin  Orbital) 

- ASW  (Augmented  Spherical  Wave) 

- LASTO  (Linearized  Augmented  Slater-Type  Orbital) 

B.  Methods  that  can  also  handle  Disordered  Structures  (Green’s  function  Methods) 

- KKR  (Korringa-Kohn-Rostoker) 

- TB-LMTO  (Tight-Binding  Linear  Muffin-Tin  Orbital) 

These  two  methods  can  be  used  within  the  mean-field  CPA  (Coherent  Potential  Approximation), 
or  beyond,  to  treat  in  an  approximate  way  chemical  disorder. 


III.  Based  on  electronic  structure  calculations  the  following  properties  can  be  obtained: 

• Equilibrium  Properties 

- lattice  parameter(s) 

- structural  energy  differences 

- crystal  structure  minimization 
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• Elastic  Properties 

- bulk  modulus 

- elastic  constants 

- phonon  spectrum 

• Electronic  Properties 

- densities  of  states  (y  of  heat  capacity  data) 

- spectroscopic  data  (ARPES,  UPS,  XPS,  PAS) 

• Transport 

- conductivity  (DC  and  AC) 

- reflectivity 

• Magnetic  Properties 

- magnetic  moments  and  magnetic  ordering 

and  Alloying  Effects  on  Properties,  in  particular: 

- Heat  of  formation  (heat  of  mixing) 

- Heat  of  transformation 

- Ordering  energy 

IV.  A series  of  software  is  now  available  to  perform  electronic  structure-based  calculations  of 

materials  properties,  e.g.: 

Pseudo-potential-based  Methods: 

VASP  (http://ems.mpi.imivie.ae.at/vasp) 

Pseudo-potentials  and  plane-wave  basis  set 

Siesta  (http://www.uam.es/depatamentos/ciencias/fismateriac/siesta) 

Pseudo-potentials  and  LCAO  basis  set 
ABINIT  (http  ://www.abinit.org) 

Pseudo-potentials  and  plane-wave  basis  set 

Full-potential  Methods: 

Wien97(http://www.tuwien.ac.at/theochem/wien97),  replaced  now  by: 

Wien2k  (http  ://www,  wien2k.  at) 

Linear  Augmented  Plane  Wave  (LAPW)  and  local  orbital  method 
LmtART  (http://www.mpi-stuttgart.mpg.de/andersen/LMTOMAN/lmtman.pdf) 
Full-Potential  Linear  Muffin-Tin  Orbital  (PLMTO) 

TB-LMTO-ASA  (http://www.mpi-stuttgart.mpg.de/andersen/LMTODOC/LMTODQC.html) 
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V.  Limitations  and  Challenges  for  Ab  Initio  Approaches: 

There  is  still  a number  of  challenges  that  ab  initio  approaches  have  to  face.  Notable  examples 

include  a full  ab  initio  description  of: 

- The  liquid  phase  of  multi-component  alloys  and  the  amorphous  state  of  matter  for  which  the 
interaction  between  the  fluctuations  of  alloy  composition  with  topological  disorder  has  to  be 
properly  accounted  for. 

- Some  systems  such  as  hydrides,  carbides  and  nitrides  for  which  large  lattice  distortions  have 
to  be  accounted  for. 

- Charge  transfer  effect  in  alloys  and  compounds  (e.g.,  oxides). 

- The  un-bias  description  of  magnetic  order  in  alloys. 

- Correlated  electron  systems  for  which  the  current  approximations  based  on  the  local-density 
approximation  (LDA)  and  beyond  are  not  applicable.  Recent  progress  made  in  this  field 
based  on  the  dynamical  mean-field  theory  (DMFT)  is  worth  noting  although  the 
implementation  has  not  been  carried  out  in  “user  friendly”  codes. 
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Presentations 

View  graphs,  additional  remarks  from  the  speakers  and  comments/questions  of  the  workshop  attendees 
are  given. 


Opening  Remarks 

Ursula  R.  Kattner,  NIST,  Gaithersburg,  MD 

John  E.  Morral,  University  of  Connecticut,  Storrs,  CT 

Excerpts  from  a recent  editorial  by  J.  Agren,  Physical  Metallurgy,  Department  of  Materials  Science  and 
Engineering,  Royal  Institute  of  Technology,  Stockholm,  Sweden  (J.  Phase  Equilibria  22  (2002)  2-3) 
where  he  describes  his  views  of  the  CALPHAD  future  were  shown: 

Today,  Calphad  assessments  usually  come  as  spin-offs  from  other  projects. 

Nevertheless,  calculation  of  phase  equilibria  and  phase  diagrams  by  means  of  the  Calphad  technique  has 
been  tremendously  successful  from  all  points  of  view  over  the  last  two  decades.  . . . 

. . . But,  what  are  the  future  needs?  What  are  the  challenges? 

All  experts  in  the  field  realize  that  good  databases  are  needed,  but  it  does  not  seem  likely  that  the 
funding  agencies  will  change  their  mind  about  databases  ...  Of  course  some  databases,  covering 
limited  subsystems,  are  of  industrial  interest  and  can  be  developed  on  a commercial  basis.  In  general 
these  databases  will  not  be  available  to  the  general  public. 

The  only  way  to  demonstrate  the  importance  of  thermodynamic  databases  is  to  apply  thermodynamic 
calculations  to  other  areas  of  science  and  technology.  And  perhaps,  that  is  the  best  way  to  get  the 
inspiration  and  stimulus  to  develop  the  field  further  and  avoid  getting  old-fashioned. 


Thus  the  challenge  of  the  next  decade  is  to  demonstrate  that  Calphad  has  the  strength  to  solve  various 
problems,  problems  that  would  be  very  difficult  to  solve  otherwise.  The  existence  of  new  software 
interfaces  that  allow  the  user  to  include  advanced  thermodynamic  calculations  in  her  or  his  software  will 
make  this  possible.” 

Agren  also  emphasized  the  future  of  applied  thermodynamics  in  such  as  Scheil  solidification, 
paraequilibrium  calculations,  and  DICTRA  simulations. 


The  objectives  of  Workshop  were  identified  as: 

• Needs  (most  important) 

• Benefits  (savings  in  cost,  energy,  environment) 

• Abilities  (current  abilities) 

with  the  goal  to  develop  a roadmap  for  obtaining  databases  with  improved  accuracy  for  practical 
applications. 
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Air  Force  Programs  in  Materials  Engineering  and  Design 

Craig  S.  Hartley,  AFOSR  (Air  Force  Office  of  Scientific  Research),  Arlington,  VA 


THE  CHALLENGE  FOR  MATERIALS 
DESIGN 
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The  issue 


Beyond  the  One  Hoss  Shay 


NIST  Workshop  on 

Databases  for  Computational  Thermodynamics  and  Diffusion 
Modeling 


Or.  Craig  S.  Hartley 
AFOSR  Program  Manager 
Air  Fore©  Research  Laboratory 
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^Materials  Science  and  Engineering  has  not  progressed 
as  rapidly  as  other  disciplines  in  contributing  to  the 
reduction  in  the  product  cycle  during  the  last  decade. 

& New  developments  in  materials  are  not  being  exploited 
as  rapidly  as  desirable  because  of  the  time  and  cost 
necessary  to  obtain  information  on  material  properties  and 
characteristics. 

('Products  are  being  designed  arid  fabricated  with 
existing  materials  having  verified  design  properties, 
resulting  in  the  use  of  less  than  optimal  materials  for  many 
applications. 


W The  Goal 


A Bit  of  History 


Unlike  The  Deacon’s  Masterpiece,  which  was  built  of  the 
best  available  materials  to  last  as  long  as  it  could  . . . 


During  the  test  forty  years,  engineer  teg  disciplines  based  on 
strusturs-prspesty  relationships  m metals  and  ceramics  have  been 
merged  with  elements  of  solid  stars  physics  and  chemistry  and 
polymer '.chemistry  to  term  She  discipline  we  call  “Materials  Science 
and  {engineering"'.  The  field  is  characterized  by  the  tetrahedron: 


STRUCTURE 


we  strive  to  make  things  using  optimized  processes  and 


materials  designed  to  meet  specifications  for  reliability 
and  performance. 


...BUT  toe  “Engineering”  component  has  nut  kept  pace  with  rapid 
developments  to  the  use  of  modem  computational  power  to  compress 
the  design  cycle  m the  development  arte  use  ©?  new  materials. 


W 


The  Present  Condition 


WHAT  DO  WE  MEAN  BY  f* 

+ PERFORMANCE?  % 

— ■» - — 


THE  SCIENCE  BASE  M THE  PLANE  OF  STRUCTURE  PRO  PERT ;£$>•• 

PROCESSING  IS  STILL  DEVELOPING,  COUPLING  EXPERIMENT  STRENGTH  IS  ONLY  ONE  OF  MANY  PROPERTIES 

Am  THEORY,  A HO  ACCELERATING  M THi  LAST  TEN  YEARS,  REQUIRED  FOR  SUCCESS : 


AUGMENTED  BY  DRAMATIC  ADVANCES  IN  COMPUTATION 


HOWEVER.  THE  PATH  TO  PERFORMANCE.  THE  ULTIMATE  GOAL 
DESIRED  BY  DESIGNERS.  IS  STILL  MOSTLY  EMPIRICAL,  EXPENSIVE. 
AND  VERY  TIME  CONSUMING 


MANUFACTURABILITY 
.><  TOUGHNESS 
SPECIFIC  STIFFNESS 
«•  FATIGUE  RESISTANCE 
& AFFORDABILITY 
w REUABIUTY 
ENVIRONMENTAL  STABILITY 

OUR  PAST  STRATEGY  HAS  BEEN  TO  DEVELOP  ONE  OR 
MORE  OF  THESE  CAPABILITIES,  OFTEN  ATTEMPTING  TO 
OPTIMIZE  IN  SEQUENCE.  AND  THEN  MEASURE  ALL  OF 
THE  REQUISITE  PROPERTIES  TO  ESTABLISH  A DATA 
BASE  FOR  DESIGERS  TO  USE... CAN  WE  DO  BETTER? 
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t J A NEW  AFOSR  APPROACH  FOR  j£f 

W STRUCTURAL  MATERIALS  \ ' 

4 j p| 

THE  MEANS  TO  A DESIRABLE  END kr 

MATERIALS  DEVELOPMENT  IS  ONE  OF  THE  GREAT 
SUCCESSES  OF  AIR  FORCE  RDT&E.  WHAT  HAS 

CHANGED  TO  JUSTIFY  A NEW  APPROACH? 

Materials  Engineering  for  Affordable 

New  Systems 

* MEANS  WILL  DEVELOP  THE  SCIENTIFIC 

PRINCIPLES  AND  UNDERSTANDING  TO  ENABLE. 
MATERIALS  ENGINEERING 

* MEANS  WILL  BE  COMPUTATIONALLY 

INTENSIVE,  BUT  WILL  BE  TIED  TO 

EXPERIMENTAL  VALIDATION 

*r  MILITARY  PROCUREMENT  HAS  DECREASED 
< INDUSTRY  IS  FOCUSED  ON  CIVILIAN  APPLICATIONS 

• INVESTMENT  IN  NEW  MATERIAL  DEVELOPMENT  IS  REDUCED 

- SYSTEM  DESIGN  TIME  HAS  DRAMATICALLY  SHORTENED 
- LONG  DEVELOPMENT  TIMES  CONTINUE  FOR  NEW  MATERIALS 

• NEW  MATERIALS  LOSE  OUT  TO  THOSE  ALREADY  CERTIFIED 

PROBLEMS  HAVE  BECOME  MORE  COMPLEX 

• WE  RE  PUSHING  THE  ENVELOPE  IN  HIGH  TEMPERATURE, 

STEALTH,  LOW  WEIGHT.  MULTIFUNCTIONALITY,  ETC. 

’ AFFORDABILITY  HAS  BECOME  A MAJOR  ISSUE  FOR  THE  AF 

*■:  COMPUTATIONAL  ALTERNATIVES  ARE  NOW  FEASIBLE 

• THE  UNDERSTANDING  OF  MATERIALS  SCIENCE  HAS  MATURED 

• COMPUTATIONAL  SPEED  HAS  REACHED  A PRACTICAL  RANGE 

ELEMENTS  OF  MEANS  J|q 

/nferoperafoiV/fy  is  the  Key  to  fT3? 

Constructive  Collaboration 

SOME  OF  THE  THINGS  WE  WILL  NEED  FOR  SUCCESS: 

MODELS  AND  EXPERIMENTS  THAT  ARE 

ADEQUATE  TO  THE  TASK 

STRATEGIES  FOR  LINKING  MODELS  IN 

DIFFERENT  SPATIAL  AND  TEMPORAL  REGIMES  - 
SO-CALLED  “MULTISCALE  MODELING” 

- OPTIMIZE  THE  COMBINATIONS  OF  AVAILABLE 

MODELS,  EXPERIMENTS  AND  PRQBABALISTIC 

DATA  BASES  TO  MINIMIZE  DEVELOPMENT  TIME 

AND  COST 

No  . - - - More 

(1  i 

- 

'mm  > - ' i|||p5x 

Spherical  . . Chickens 

- ...  - 

LINKS  TO  OTHER  PROGRAMS 

The  Challenge 

- MEANS  BUILDS  ON  PROGRAMS  IN  MANY  AGENCIES 
FOCUSED  ON  “MATERIALS  BY  DESIGN" 

-'MEANS  BUILDS  ON  PROGRAMS  IN  MANY  AGENCIES 
FOCUSED  ON  “MULTISCALE  MODELING” 

- MEANS  WILL  DEVELOP  THE  SCIENTIFIC  BASE  FOR  AND 
DERIVE  INSIGHT  AND  DIRECTION  FROM  THE 

6.2  PROGRAM  IN  DARPA  CALLED  AIM 
(ACCELERATED  INSERTION  OF  MATERIALS) 

- MEANS  WILL  FOCUS  THE  6,1  COMMUNITY  ON 
KNOWLEDGE  BASE  CONSTRUCTION.  PROPAGATION 

OF  ERRORS.  UNKING  OF  SCALES.  AND  NEW,  EFFICIENT 
EXPERIMENTAL  APPROACHES 

if. 

Exploit  computational  Materials  Science  and 
Engineering  to  develop  techniques  for  coupling 
models  of  material  behavior  to  design  software, 
enabling  materials  design  to  be  an  integral  part 

of  the  global  design  process. 
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CURRENT  STATUS  OF 
COMPUTATIONAL  METHODS 


S&Jmeans  Projects  in 

9S 


Metallic  Materials 
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* COMPUTATONAL  DESIGN  OF  ADVANCED 
AEROTURBINE  MATERIALS:  PI  -G.  Olson, 
Northwestern  U. 

^DEVELOPMENT  OF  A PHYSICALLY  BASED 
METHODOLOGY  FOR  PREDICTING  MATERIAL 
VARIABILITY  IN  FATIGUE  CRACK  INITIATION  AND 
GROWTH:  PI  - K.  Chan.  Southwest  Research  Institute 

^ MICROSTRUCTURE-BASED  MODELING  FOR  LIFE- 
LIMITED  COMPONENTS:  PI  - H.  Fraser.  Ohio  State  U. 

^DEVELOPMENT  OF  AN  ACCELERATED 

METHODOLOGY  FOR  THE  EVALUATION  OF  CRITICAL 
MECHANICAL  PROPERTIES  OF  POLYPHASE 
ALLOYS:  PI  - P.  Dawson  & M.  Miller,  Cornell  U. 

,x:MURI  on  Design  of  Multifunctional  Materials:  TBA 


Computational  Design  of  Advanced 
«+*  Aeroturbine  Materials  (Olson,  NWU)  - 


v J DEVELOPMENT  OF  A PHYSICS-BASED  METHODOLOGY 
FOR  PREDICTING  MATERIAL  VARIABILITY  IN  FATIGUE 
CRACK  INITIATION  AND  GROWTH  (Chan.  SWRI) 


I? 


A hierarchy  of 
computational  models 
will  be  integrated 
through  computational 
thermodynamics  to 
design  a metal/ceramic 
system  that  addresses 
control  of  oxygen 
behavior  both  in  a BCC 
Nb-based  matrix  as 
well  as  in  stable  oxide 
films  with  controlled  ™ 
expansion  and  v>3 
adherence. 


Objectives: 

* To  develop  physics-based  fatigue  crack 
initiation  and  growth  models 

* To  develop  a probabilistic  approach  for  linking 
physically  based  models  into  a continuum 
framework 

* To  demonstrate  the  utilities  of  the 
methodology  in  a probabilistic  design  setting 


ATOMIC  SCALE 


DISLOCATION  SCALE 


1}  The  friction  stress  for  irreversible  slip  during 
fatigue  will  be  modeled  at  the  atomistic  scale 
using  the  Peierls-Nabarro  model  and  the  thermally 
activated  flow  approach. 


2)  The  formation  of  dislocation  cell  substructure  will 
be  modeled  by  considering  interactions  of 
dislocation  pile-ups  and  cel!  walls  at  the  dislocation 
ceil  size  level  using  results  from  irreversible  slip 
modeling. 


Moth-da* 
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MICROSTRUCTURE  SCALE 

«*§*■  v%^ 


3)  ft/Sicrostructure-based  fatigue  crack  initiation  and 
growth  models  will  be  developed  at  the  grain  level 
using  results  from  dislocation  structure  modeiing; 

krwfeshdEtdir  MioMiug 


.Vlx  iNwipK  fete 

kiiiyvK  Cicidt  Jreiw.xu  %k«£rf 

i Man 

| l-ilivu 

Ci*9&  Umik  kkd  v'l 

5 

\ 

S v ?<lictW«*Kwt 

1 M 

| dN 

//  s 

{S:  Stress  Range:  N,:  Fatigue  Life;  da/dN:  Crack 
Growth  Rate;  ? K;  Stress  Intensity  Factor  Range) 


applications 

Potential  applications  of  the  probabilistic  fatigue 
models  include  component  design  and  life-prediction 
analyses  that  include  material  variability  and 
confidence  limits  for  the  fatigue  properties. 


Mlfjosl. iH\ttjrr~haf.ini  D .ni.-r-  -■ 


Qumbbstive  cbatacfeittabon  wsed  Ut  describe  fflicioslrijclufe 


Maura!  ratorotfes  uaari  ta  reveal  fuftetiotvri  dependencies  ot 
properties  on  Jtticrosffuciur&l  features 


Production  oJ  vartirttofis  of  osfcrostr  octree  torprfbej  with 
property  assessment  provides  aatassary  databases 


This  iftfenristxm  wiff  t>e  used  to  devakip  phy sically-bssed 
motteis  for  predic  tkm  of  fatigue  propertfes 


CONTINUUM  SCALE 


W Microstructure-based  Modeling  for  Life- Urn  tied 

Components  {Haroish  L,  Frsser.  Somnath  Ghosh,  James 
Larsen".  Yunzhi  Wang.  James  C.  Williams;  CAMM-QSU 
^ >APRU 
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APPROACH  (CONT'D) 


re 


Compiwivon  at  Phase  Fh#k3  smuriailon 
piederiitm  ffrrffjwth  oKpoftmental  observation 
jiircffti.  courtesy  oft.  Bendersky*  ctf 
mkTcstfucturo  termed  during  the  D015  ? O* 
phase  transformation  in  Ti-Al-Nh 


Phase  FlaW  model  wfff  *>t»  used  as  the  primary 

coriiptitatiwKjf  rtret  hod  to  tMtriMtate  the 
mtcrostructiiral  OYOfutton 

fHo drtinq  vrill  indudfc  tire  coupling  t>l  growth  to  IIk- 
cSffoslon  fia&fs  erf  minority  alteytag  elements, 
ttccotrjmodgfion  et  coherency  strain,  and  anisotropy 
in  mtorfociaf  A plain  boundary  energy  & mobility 
Oevefoptnpm  of  constitutive  oqtMSHWsMfk.Jem 
re&reed-ortier  models  for  grate  growth  and  overall 
transfoiraation  kinetics  biciudmg  micfsatiaa  growth 
and  coarsttnunfj  ^ 


4)  The  microstructure-based  fatigue  crack  initiation 
and  growth  models  will  be  integrated  into  a 
probabilistic  framework  at  the  continuum  level. 


Aim* 

the  davafeoptmirt  of  mtiuosJrticfui  databases  for  the  Alloys  Ti-fc-4  and  1 pfj-2-4-2 

A dewrcmnatfon  erf  the  proj*erty<pftltolti«g 

mimosrf  wot  i»al  feafuies  tnfluencmg  low  syefe  fatigue 
afld  fatigue  crack  growth  fn  Tt-C-4  and  TI-&-2-4-2 
A robust  roefbtMfofogy  tc<r  quantitetryety  tfctemtmwig 

nticr  structural  features  and  th«r  repr&santation  in 
modeiing  and  sfrrariatfett 

Quaterttitve  satiulatioit  tnerftevrlofoniev  tor  the 

prediction  ef  th#  development  of  micros tiructurad 
Mtrvti,  which  are  key  to  influencing  tCF  emit  f&v’ttkl, 
as  a fttnctteb  of  heet-treatmem 

Tlte  dsw-tepment  o-f « serf  of tinseostHic fure-'-beoett models  for  hie  pre-diction  of  LCF  and 
dUvdN  m TMvJ  and  Tf-b-2-4-2 

The  proviatan  or  a act  ot  r£M-baa*d  lords  for  fh«  anaTyata  oftita-Iriniting  foatrsas  itt  turbine 
rc-ior* 


PrYveiott-mant  ot  Physfe-.Hvriy-foaa>>d  Models  for  Fat^oo 

Acfo«v*a  capadslfify  fw  ptodictfng  tho  vartahon?  in 
firtkjtw  performance; 

Lrnk  these  vw  rations  to  miorosbuctural  pammeiers 
Hw  ouyh  the  rnfeyratKKi  a I mociuioistK^fly  rc-airstk. 
aub-moefef*. 

Buiid  cm  the  nourei  networks  rcrfaJin^  microstrucfuje 
to  Mtrp.ie  & Creek  yaowlh  rate  iietiavlof  to  develop 
phy*ic^ifty  haserl  morlefc 


Adapt ivo  Finite  Fieragnt  fee-detinn:  Analysis. <»f  Ufe- 

ti»;-r»K>  Festufes 

Coupled  multifile  scale  sartuSatifst  of  Jhe  rfefonrialmn 

and  fraetute  pioc^sses  of  mufti-coiony. 

Ti  afktys  wifi  bo  developed 
Tlte  nwif>*s>:ale  compumfoaal  system  vdll  create  a 
fwarchy  ot  cnmpniatKKKiJ  sub-domains  providing 
ne«e*cary  rwolotten  in  proefletitki  defonnatioft  and  the 
evolution  of  damage  drtd  Iraclwre 


dtttpw:  Fi'diuWtMi*:  IxPmwc  I*»«4kiw# 

toclwilJJJ  S1k»  iKHBWrrBXr  1 iHIHiafcll 
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Accelerated  Methedoiogy  tor  Evaluation  Of  Critical 
Properties  in  Polyphase  Alloys 
IP.  Dawson  A HI  Milter,  Corn  oil  U j 


^Objectives 

-Rapid  evaluation  of  stiffness  and  strength  ot  polyphase  metallic 
systems 

-Reduced  time  for  insertion  of  alternative  materials  in 
mechanical  design 

- ^Methodology: 

•Develop  protocols  for  a suite  of  simulations  and  experiments  to 
assess  elastic  moduli  and  anisotropic  yield  surfaces 
•Deploy  around  the  Digital  Materia!  framework 
•Interface  required  diagnostic  tools  via  the  Digital  Material 

-Diagnostic  Tools: 

•Simulation:  grain-by-grain  finite  element  models  of  poiycrystals 
•Experiments:  mechanical  tests  and  in  Situ  diffraction 
measurements 

•Visualization:  advanced  graphics  as  interpretation  aids 


Experiments  and  Simulations 


mww.  | 

X 

j *ieKVrtdk>ni}T-.'ifls 



; ] firam  j 

.1  (Til 

1 

tSE 

Digital  Material  - 
Collaboration  Environment 


The  Vision 


| Geomd.nc 

AUrilwt?* 

| - comjxjweai; 

- uor  MODI- 

j - intuit* 

\ - jxmictex. 

- COm|W$ttKHVS&:  { 

l f r,o  if  v)  nr’- 


liibk 

M* 
■ ■ > 


^ with 

|§  work! 

j Simulation 
i i&peromsRi  j 


. J | . , MODF 

Grain  size 

mm  r 

IbtiJums 

i 1H) 

V J 

***  i 1 00) 

/ribs  1 « • r x tv  : 

* The  periodic  table  is  21. 
the  ultimate  data 
base 

* Imagine  a design 
space  that  extends 
from  the  periodic 
table  to  input  into 
current  design 
software. 

* What  do  we  need  to 
do  to  fill  the  gaps  in 
this  space? 


Ps&&c  T«sli  o#  Seretrtj 


Following  the  Vision 


The  Payoff 


Optimal  utilization  of  materials  and  processes  to 
produce  affordable,  reliable  and  durable  products  for 
military  and  civilian  applications. 
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Additional  points  made  by  speaker: 

• To  include  computational  materials  as  a design  component  the  ability  to  compute  design  properties 
is  needed 

• Materials  computational  power  is  not  up  to  speed  to  use  to  replace  existing  materials 

• Engineering  part  of  materials  science  has  fallen  behind;  need  to  make  connection  to  performance 
(strength,  toughness,  fatigue  resistance,  etc. ) 

• MEANS:  Materials  Engineering  for  Affordable  New  Systems 

♦ Tied  to  experimental  validation 

♦ Need  models  (robust) 

♦ Need  methods  to  link  models 

♦ Need  methods  to  optimize  models 

♦ Interoperability  is  key 

♦ Focus  on  basic  science  on  knowledge  base 

♦ Make  materials  science  an  integral  part  of  global  design  process 

Questions/comments  from  workshop  attendees: 

• Is  OOF  related  to  the  OSU  program?  — No,  not  formally,  yes  philosophically. 

• Performance  - structure  relationship  not  well  correlated  yet 

• Need  to  get  manufacturers  involved  in  model  development  or  they  will  not  trust  design  process  and 
use  the  new  materials. 
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Application  and  Needs  Assessment  of  Thermodynamic  Databases  at  Howmet  Castings 
Ty  W.  Hansen  , Howmet  Research  Corporation,  Whitehall,  MI 


Application  and  Needs  Assessment 
of  Thermodynamic  Databases  at 
Howmet  Castings 


Application  and  Needs  Assessment  JsL 

of  Thermodynamic  Databases  at  Howmet  Castings 


Introduction 


Hawmci  Castings  Background 


Tyros  W.  Hansen 
Howmet  Research  Corporation 
3/21/02 


Alcoa  Industrial  Components  Group 


^Current  Capabilities 
^ Howmet  Applications 
X Howmet  Thermodynamic  Needs 


Howmet 

Castings 


a 


Overview 


rtorffnei 

Casluwjs 


SL 


M Howmet  Part  Alcoa  Industrial  Components  Group 

H Based  in  Salt  Lake  City,  Former  Cordant  Headquarter’s 
f&  Group  Also  Includes: 

M Kuck  Fasteners 
^ Alcoa  Forgings. 

M Alcoa  Automotive 


Total  Revenues  - $3.48 
42  Locations 
1 7,030  Employees 


Airframe  Market 


Howmet 
Cast*  Tigs 


fH  World  Market  Leader  in  2 
Very  Attractive  Markets 

M Leadership  Position 

Differentiated  & Deferrable 

H Technology  Leader 
Si  Operational  Sxccitence 

B Bxcetot  Customer 
Relationships 

H 29  Manufacturing  Facilities  in 
5 Countries 

H Over  It, 000  Employees 


Major  Customers 


1^93  Total  Revenues;  Million 


m 

4€% 


Other 

4% 

8% 


\ 

Aero  Engine 


Hcwmet  ** 
Castings. 


Aircraft  Components:  Ai  & Ti 
Doors,  Nacelle.  Strut,  Right  Control. 
APU  and  Environmental  Systems. 

Missile  Components.  A!  £ Ti 
Missile  Bodies.  Pins.  Optical  Seeker 
amt  Fuel  Component*.  Electronic 
Boxes/frames. 

Rocket  Components:  5iA 

High  Pressure  Fuel  and  LQX  Pump 

and  Turbina  Housings* 


Aft  Engit  -m  Mount 


Aero  Engine 

Wi  Pratt  & Whitney  Aircraft 
W.  Pratt  A.  Whitney  Canada 
M.  Rolls  f?»y cp 

f§  General  Electric  Asrcratt 
Engine 

M Honeywell 


Industrial  Gas  Turbine 


Siemens  Westinghouse 
Alstom 

General  Electric  Power 
Generation 

Solar 
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Application  and  Needs  Assessment 
of  Thermodynamic  Databases  at  Howmet  Castings 


Current  Capabilities 

3i  Howmet  utilizes  three  so i vers; 
B t»roCAST  * UL 
B JMatFro  - m 
§§  Panda!  - Ml 


Howmci  S5 
CAAtJtigs 


Application  and  Needs  Assessment  CmUbj*  ®- 

of  Thermodynamic  Databases  at  Howmet  Castings 


1N71H 


Application  and  Needs  Assessment  cSng» 

of  Thermodynamic  Databases  at  Howmet  Castings  ' , 

IN71S 


Application  and  Needs  Assessment  C«Wng*  -®~ 

of  Thermodynamic  Databases  at  Howmet  Castings 


INTIS 


I’amJal:  3 miutites 


Application  and  Needs  Assessment  cKngs  S. 

of  Thermodynamic  Databases  at  Howmet  Castings 


Howmet  Applications 

I Heat  T real  Optimization 
i solvus,  solidus 
: Ailoy  Development 

solvus, i volume,  liquidus/solidus.  phase 
evolution  and  stability 
I Attoy  Targeting  tor  Castabiliiy 

mf  solvus,  iiquiduslsolidus,  phase  evolution  and  stability 


Application  and  Needs  Assessment  Caffign  & 

of  Thermodynamic  Databases  at  Howmet  Castings 


Benchmark  SX  Allot 
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Application  and  Needs  Assessment 
of  Thermodynamic  Databases  at  Howrnet  Castings 


Howmet  G 


Application  and  Needs  Assessment 

of  Thermodynamic  Databases  at  Howmct  Castings  " '■  • ' 


II  \-tXU 


Temperature  fF) 


Howmet  Thermodynamic  Database  Needs 
S3  improved  Accuracy 

S Optimized  about  typical  alloy  compositional  limits 
m Density  Model 

® Model  diffusion  in  the  solid  state 


Summary  from  workshop  organizers: 

Products: 

• Turbine  blades  and  vanes  (IGT):  directionally  solidified,  single  crystal,  equi-axed 

• Airframes  (Al,  Ti) 

• 88%  of  the  production  are  turbine  parts,  35%  are  spare  parts 
Success  measures:  Price,  quality,  delivery 

Current  modeling  capabilities: 

• ProCast:  finite  element  software  package  for  simulation  of  casting  processes,  includes  lever  rule  and 
Scheil  solidification  calculations:  Ni,  Ti,  Fe,  Al  databases  from  Thermotech 

• JMatPro:  software  package  for  the  calculation  of  physical  and  mechanical  properties,  includes  lever 
rule  calculations:  Ni  database  (Thermotech) 

• Pandat:  phase  equilibria  calculations:  Ni  (not  using,  too  difficult,  long  calculation  time) 

• Test  phase  equilibria  calculations  for  IN718  had  different  run  times  and  gave  different  results 

• Heat  Treatment  Optimization 

♦ Calculation  of  y’  solvus,  solidus  (motivation:  sell  reduced  heat  cycle  to  customers,  need  more 
experimental  verification) 

• Alloy  Development 

♦ Calculation  of  y’  solvus,  y’  volume,  liquidus/solidus,  phase  evolution,  stability 

♦ Objective:  improve  castability  (customers  - small  base;  non-patented  alloys) 

♦ Use  lever  calculations,  trend  analysis  (not  accuracy) 

♦ Patent  conflicts 

♦ Stress-rupture  calculations 

• Concerns/  Needs 

♦ Are  we  using  the  tools  correctly 

♦ How  do  we  make  the  tools  better 
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♦ Database  needs 

■ Improved  accuracy 

■ Density  model  (liquid)  porosity 

■ Add  Pt,  Si,  Pd  to  database  for  coatings 

■ Model  diffusion  in  solid  state 

Questions/comments  from  workshop  attendees: 

UES:  Si  added  in  latest  version  of  Ni  database 

Alloy  development  restricted  by  composition  limits  on  thermodynamic  database 
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Computer  Simulations  and  Database  Needs  for  Energy  Savings,  Zero  Environmental  Impact  and 
Quality  Improvements  for  the  Heat  Treating  Industry 

Richard  D.  Sisson,  Worcester  Polytechnic  Institute,  Center  for  Heat  Treating  Excellence  (WPI,  CHTE) 
Worcester,  MA 


Computer  Simulations  and  Database  Needs  for  Energy 

Savings,  Zero  Environmental  Impact  and  Quality 

Improvements  in  the  Heat  Treating  Industry. 

Research  Needs  identified  in  the  Heat  Treating  Technology 

Roadmap  Workshop  -6-7  February  1997 

Quenching  Technology 

R*Ck  Sis&oti 

Moh«mm«d  ftUmtru£Si*rmin 

Modi tis  lor  hast  transfer  hehav far  in  euenett  bath s iiat  ivti!  ensure  uniform 
costing  of  a rang*;  sf  loads 

Process  Modeling 

Materials  Science  & Ensioeeirirtg  Program 

Mechanical  Engineering  Department 

Worcester  Polytechnic  Institute 

Worcester,  MA  01608 

miaasmjeture  response  models  based  on  aiioy  composition,  atmosphere, 
temperature,  end  time 

models  lo  predict  distortion  and  residual  stress  profiles 

database  of  thermal  and  mechanical  properties  up  to  and  including  neat 

treatment  temperatures 

£>&&«*£  mpheitei 

models  for  continuous  cooling  transformations  (CCT)  of  tiers!.  trastabSa  aloys 
models  foressntinmus  beating  transformations  (Off)  of  near  treatable  alloys 

Software  Packages 

predicat'd  software  chat  heal  treaters  can  use  ti>  compare  ants  select  furnace 
fmm  diffemn t suppliers  (a  standard  method  of  predicting  fomaoe 
variability  is  incorporated) 

ttX-iK&X.  .< *r  te  .AWJaC KWvl  -•][— *~Mt — • Kf**^ 

»»5*<  » WBfWWI  Yfi  ' >(5$-T  K g g— 

Ws?  mrrswwe  ;«?:?««  A 

itfPi  CHTE 

w ITi  »««»»*  w»s»m  mrirsrt fmrrm it  -S- ■ 

General  Needs  for  Roadmap 

~ . .... 

Proposed  Research  from  Roadmap 

Heat  Treating  Processes  that  allow  a shorter  cycle  time  and  require 
tower  cost  equipment. 

Alternative  quenching  media  that  are  more  environmentally 
friendly,  and  that  transfer  heat  more  efficiently, 
t improved  process  sensors,  including  those  for  carbon  content, 
residual  stress,  and  cleanliness.  Also,  more  advanced  controls  that 
fully  exploit  there  and  other  sensors. 

Enhanced  computer  modeling  of  processes,  which  include 
composition,  distortion,  resultant  microstructures,  and  final 
properties. 

More  effective  dissemination  of  solutions  to  specific  problems. 

Goal  A:  integrated  process  models 

Qssendyng  Models 

Sectrerwagnstic  rnodels 

i Mechanical  Models 

Transfocmaissn  Databases 

Predict  Heat-up  parameters 

Predict  carburization,  retriding,  and  soluBonlring  behavior, 
t Predict  oosfag  behavior 

Materials  Databases  and  process  inputs 
beat  transfer  coefficients  (heating  & edging) 

physics!  property  data  (thermal  expansion,  speeoffc  tot,  thermal  conductivity, 
density} 

Bechtel  and  Magnetic  property  data 

Ht’TKHI*  r *»:>?»«  WitMK  /f 

Vc-^xwf- «-  #swj«  <*<e  »ws>>(  >/*  . I, u,  i. 

rto>»~vv»x> ACHTC 

Modeling  Heat  Treating  Processes 

An  Energy  Savings  Model  for  the  Heat  Treatment  of  Castings 

DOE  - funded  , 4 year  project,  WPI  & U.  Conn  & U.  Mass. 

Finite  Difference  Models  (FDM) 

Finite  Element  Models  (FEM) 

Computational  Fluid  Dynamics  Models  (CFD) 

Data  Needs  - 

thermal  properties  - k,  cy  «.  h f(T,  composition} 
diffusion  coefficients  - D.s  (T,  composition) 

phase  transformation  rates  f(T,  composition) 

elastic  and  plastic  properties  f(T,  composition) 

Develop,  verify  and  market  an  integrated  system  of  software, 
databases,  and  design  rules  to  enable  quantitative  prediction  and 
optimization  of  heat  treatment  of  aluminum  castings  to  increase 
quality,  increase  productivity,  reduce  heat  treatment  cycle  times 
and  reduce  energy  consumption. 

Seduce  cycle  times  for  saiutiontiiRg  cast  aluminum  aiteys  from  12  to  2 hours 
eeerease  energy  consumption  Say  more  than  56%. 
estimated  energy  savings  of  ? triiton  STU's  in  £0  years! 

Computer  Modules 

HTfTiSHACE  prethot  local  tempemiute  cycle  in  a past 

PRGAQ- optimizes  pari  tadsrsg 

TPCQNTBOi.  - optimizes  f/r.vyarrfme  eoniro!  in  fynm* 

■ QtitifCHCAlC  - predicts  caafej  rates  ami  «icK5SJsucti»«s  as  a function 
(ivendtiryj  process  parameters  and  position  in  She  part 

WPl^».n.  ^ - /-CHTE 

wm  JSsssaAszs^L,  CHTE 
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An  Energy  Savings  Mode!  for  the  Heat  Treatment  of  Castings 
DOE  - funded  , 4 year  project,  WPi  & U.  Conn  & U.  Mass. 

Sofutiomzing  Prediction 

PHASECAI.C  phase  diagram  for  industrial  alloys 
DIFFCALC  - difrusivities 

CASI5C6  - predicts  local  segregation  as  a function  of  alloy  composition 
and  thermal  cycle 

CASTSOlf*  presets  to!  segregation  after  soltitionfeing  treatment 

Property  Prediction 

PROPERTIES  ■ predicts  iocai  values  of  strength 
POSTPROCESS  • presents  results  to  facilitate  decision  malting 


Databases  for  heat  treating 
Databases 

Heat  Transfer  and  Quenching  performance  database 
heat  csosdty, 
tlvermal  conductivity, 
emissMdes, 

beat  transfer  coefficients 
: Phase  diagrams  for  selected  alloys 
: Diffusivities  versus  composition  and  temperature 
Properties  versus  microstructure 

(i.e.  porosity.  Si  particle  size,  umfissoivad  8 arid  grain  size) 


mm 


fSi»rfCt8tC  Xf 


iU 


XT7Z&SP? 


CHTE 


»38 ii&U*  HirifHId  : 


S/*!78CK? 


/•chte 


4140  Steel  IT  diagram  showing  two  quench  path 


Bubble  dynamics 

CHTE  4140  probe  in  mineral  oil 


WH, 


xxerxv?  ax  too 

*>x-syj«wxx>  ox*  6:Km»>  MHUi-yt  «i*7 


Characterization  of  Quenching  Fluid 


CHTE 


WH 


*S83f»Ttg  »<U«  fifth  <10 


940  ^ 

*n  4>.»rt»w  ftr  »««>■*«  fx« 


/CHTE 


Specific  Heat  vs.  Temperature  for  304 
Stainless  Steel 


* Heat  transfer 
coefficient  \$ 
determined  using 
bmpcd  parameter 
analysis 


dr  /?.t 

4t  pvc. 


= pftyr 

lZ-rA 


/-CHTE 


Sources:  Handbook  of  Stainless  Steel  and  TPRC 


0 200  400  5S5  500  1000  t£G0 

Temperature  *C 


/CHTE 
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Heat  transfer  coefficients  of  SS304  probe 
quenched  in  mineral  oil 


Cooling  rate  4140  steel  probe  as  a function 
of  temperature 

quoncfte-d  in  mineral  oh  (DHR8&A) 


/>|3HTE 


WH. 


vwtm 


CHTE 


Heat  transfer  coefficients  vs.  temperature 

SS30A  probe  qwt’nchotf  in  different  q went  hunts 


Heat  transfer  coefficients  vs.  temperature 


qwnehed  ♦ n miner#!  off  (HR88A) 


WM 


*»**!*»*  O'  ?‘;T* 


/CHTE 


CHTE 


QUENCHPAD  - CHTE  Quench  Database 


Quenching 
Database  - 
Flowchart 


]=> 


c=){ 


SR€JtU*  r*bT?lS«*tC  *I»':rvn 


Iftit/Mtot#  «•  to  i 

frwaatwita  »v5  x 

aammu 


CHTE 
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Additional  points  made  by  speaker: 

• Heat  treating  modeling  limited  by  available  databases 

• Software  uses  look-up  tables  as  databases,  source  of  potential  problems:  data  such  as  cp  are  sparse 
and  differences  between  different  sources  have  a trickle  down  effect  on  derived  properties 

Questions/comments  from  workshop  attendees: 

• Calculation  of  phase  equilibria  could  solve  part  of  the  database  problem 

• Quenching  modeling  will  be  geometry/equipment  specific 
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Materials  Programs  at  NSF 

K.L.  Murty,  NSF  (National  Science  Foundation),  Arlington,  VA 


NSF  Merit  Review  C riteria 

What  is  the  intellectual  merit  of  the  proposed 

activity '? 

What  are  the  broader  impacts  of  the  proposed 

activity  (including  specifically  the  value  added  by 
the  international  cooperation)'7 


OM!R -Office  of  the  director 


Advanced  Materials  & 


T«tw»o>0g?  Kwftb&ng 
Oudfr 


5* leave  Cluster 


MS¥  PY  2©§2  Budget  Request 

By  Priority  Areas 

Millions  of  Dollars 


FY  2001 
Plan 

FY  2002 
Request 

Percent 

Change 

Biocomp,  in  the  Env. 

$ 55 

$58 

5.9% 

Information  Tech. 

$259 

$273 

5.0% 

Nanoscaie  S&E 

$150 

$174 

16.1% 

Learning  for  the 
21st  Century 

$121 

$126 

3.3% 

NSF  Support  for  Materials 

approximately  $3®§  Million  Annually 


NSF  f Y 2H§2  Budget  Reguest 

By  Strategic  Goal 

Million^  of  Dollars 


FY  2001 

FY  2002 

Percent 

Plan 

Request 

Change 

People 

$888 

$1,002 

12.8% 

Ideas 

$2,251 

$2,220 

-1 .4% 

Tools 

$1,061 

SI. 024 

-3.5% 

Admin.&  Management 

$216 

$227 

5.0% 

Total,  NSF 

Total  i 

$4,416 

ovecl 

$4,473 

; $4,789 

1.3% 

8,4%  $ 

NSF  Goals 


'Metals  Research  at  NSF 


K.  Linga  <KL>  Mum 

Program  Director  flPA),  Metals  f&earcfe,. 
Nstimml  Science  FoasMahots 


\ > 1 Stuekcc  «Kvr»,v  .>.u: 

{m:  •A?)  d ,V  i 

vlNPy  N«rtfc  { ftndiou  Stole  I nivcrsilv 

Rifteigli  M'  l SA 


Yfi 
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Types  of  Individual  Investigator 
Grants 


Regular  grants 

Focussed  Research  Group  (FRG) 


Grant  Opportunities  for  Academic  Liaison 
Industry  (GOAL!) 


Small  Grant  for  Exploratory  Research  (SGER) 
Technical  Conference  Support 


MMitiafs  Pm§mm 


Emphasis  on  creating  a framework  for  the  development 
of  advanced  metallic  materials 


1 8 WSGCSSSSNO 
; ® CHARACTERIZATION 
f BMOOBJNS 


* New  Metallic  Systems 

* Bulk,  thin  films,  wires  sal  particles 
» Characterization 

» Processing 

* Environmental  Effects 


mm 


Emphasis  on  creating  a framework  for  the  development 

of  advanced  metallic  materials 



* Hew  l^ietaflsc  Systems 
» Bulk,  thin  films,  wires  a«cl  particles 

* Characterization 

* Processing 


B PROCESSING 
8 CHARACTERIZATION 
8 MODELING 
j 8 PROFESSES 


BMR  METALS  RESEARCH 
PROGRAM  FUNDING 


SiVJ 


12 

11.5 
11- 

10.5 
10 

9.5 
9- 

8.5 
8 


I I I I I I I I I 


i £ £ 


£ £ t £ £ £ 


FY  2001  METALS  RESEARCH 
PROGRAM 

—f i==— — 

Fasdisg  Actions  (budget  SI  I..1  HM) 

❖ 29  research  grants  (£U(2?tO 

<■  "0  amthtmiiom  ($7<498K) 

->  4 t-  (inferences  i S25K ) 

>5-  3ft  B.Et  Mudenis  (SI43K) 

<■  8 other  sopplwiiHftd*  {equip.,  vie./ ($J82K) 


Success  rales  for  proposals 

> 28%  regular  ! I \ proposals 

> 44%  CARPER.  proposals 


Materials  Research  and  Education 


# 
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HAPrograms 


Focussed  Research  Group  (FBg| 

Typically  $180  - $250K/yr  for  3 years  and  3 to  4 PfVl 

The  whole  is  greater  than  the  parts 

Mail  reviewed  (co-review  with  other  programs/divisioi 


i one,  but  August  - Oci 


HA  Programs 


HAPrograms 

— » — ■ — — — ”~4 

Grant  Opportunities  for  Academic  Liaison  with 

industry  (GdftLI) 

->  Typically  $90  - $125K/yr  for  3 or  4 yea^w; 

2 Pi’s  (one  from  the  industrial  company) 

* Strong  industrial  interaction  is  a must 
: Mail  reviewed 

v Deadline  - none,  hut  August  - October  preferred  foWMjBfc 
reserve  funds 

Small  Grant  for  Exploratory  Research  (sdfl| 

v Limited to$  100K  max.  for  1 year  V 

<-  Usually  high  risk,  innovative  ideas 
•>  PD  reviewed  / Target  date:  November  I 


FY  2002  Education  Proposals  IGERT 
(NSF00-78) 

X 

Integrative  Graduate  Education  and  Keseahshl'raitteeshtp  (IGERT) 


(sails  in  June) 


g for  new  models  of  graduate  training  witJf 
sdplinaty  research  themes 


Faculty  Early  CAREER  Develop! 

- Typically  $75  - $IQ0K/yr  for  5 years 
•:U  i >■  Usually  panel  reviewed  . 


FY  2002  DMR  Education  Grants 

OMR  supports  innovative  education  approaches  to 
materials-related  education 


•get  date:  November  1. 2001 


K.  Ltnga  <KL)  Muriy 
research  experience  - maim!  theme  > 


w v Ralston  - Defect  * 


Defects 

Characterization 

Kinetics 

& 

Microstrueture 


V.  km  yj  i 1 

';■■■■■  - "Cl 


1 1 jfe-IVedivtksu  and  I ill:-!  ^tension  of  Materials  in  .vemee  1 


kfH'tm-h  sttpfwtMt  hr  \St . iH)t.  t,£.  \IC\C.  ti'Hl , M tSSf] 
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Materials  Research  Science  and 
Engineering  Centers 
(MKSECs). 


o Sept  2601 ; a new  MRSEC  competition  began;  81 
prcproposals  were  received, 

❖ The  aim  is  to  select  36  candidates  for  full  proposals 
(due  January  14, 2602), 

❖ Reverse  site  visits  - May  2602 

❖ ~ 10  - 15  awards  (est.  825M  available) 


NSF  ABSTRACT  SEARCH 


**  ■■■  '■  ■■ 


he  N$T'-WER  site  Isas  a wealth  oHrffwrmatioii  at: 

For  jjjptraef  searefe?  go  u>  {})  scarclh (2)  fieldll  Ifcccfe  ant 
award  atetoeSs.,  (3) *  *et  date  Option  u>  ” d<»  r ltMN-9'3 
sstul  (4)  NSF  program  option  so  “e«o»ata$  METALS,*4 
ISft  trim*  arc  sliOWfe,  PIN,  institution*.  amhuats,  abstract* 


MEW  NSF  PROGRAM  ANNOUNCEMENTS 
Try  the  Custom  News  Service. 

To  sign  up  go  to: 

httpT/w  ww , ns  Ego  v/itome/cns''s  tart.]itm 


Instnunentation  ior  Materials  kfes^areh  (IMR) 
Smaller  groups:  avg.  award  size  $K©K. 
Deadline:  early  January  2002 
Major  Research  Instrumentation  (MR1) 

- Larger  multidisciplinary  groups:  avg,  award  H 
Deadline:  late  January’ 2002 


New  Programs  for  International  Coope 

US  International  Materials  Institute 

To  enhance  US  materials  research  and  educaticm 
interactions  with  foreign  investigators 
New  interactive  research  activities 

US/Europe  US/Canada  US/Mexico 


mm 


Additional  points  made  by  speaker: 

Division  of  Materials  Research  (DMR)/Metals  Program  at  NSF: 

• 01-02  funding  increased  by  8%  (majority  of  increase  to  bio-,  nano-,  and  information  technology) 

• Mission:  People  - Ideas  - Tools 

• Types  of  Grants: 

♦ Regular 

♦ FRG  (Focused  Research  Group) 

♦ GOALI  (Grant  Opportunities  for  Academic  Liaison  with  Industry) 

♦ SGER  (Small  Grants  for  Exploratory  Research) 

• Metals  funding  stays  at  the  same  level  at  around  1 1 million  dollars 
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Microstructure  - Process  - Property  Models  for  the  Aluminum  Industry 
Joanne  L.  Murray,  Alcoa  Technical  Center,  Alcoa  Center,  PA 


Microstructure  - Process  - 
Property  Models  for  the 
Aluminum  Industry 


Joanne  Murray 
Alcoa  Technical  Center 
2002  March  2 1 


Production 
of  Body  Stock  Alloys 

courtesy  of  Tom  Rouns,  Alcoa  Technical  Center 

Scalping  removes  the  as-cast  surface  of  the  ingot 

Preheating  heats  the  metal  for  hot  rolling  and 
produces  the  correct  microstructure  for  hot  rolling  and 
subsequent  operations 


Aluminum  alloys  & products 


Non-heat  treatable  - 

Work-hardened 

Heat  treatable  - 

Precipitation  hardened 

Other 

Ixxx  AJ-Fe-SI 

•Commercial  purity  A1 

•Litho  sheet 

•Bright  sheet 

6sxx  AJ-Mg-Si 

•6061  - Medium  strength 

structural  alloys 

•6063  - soft  alloy 

extrusions,  architectural 

applications 

•6262  - free  machining 

Casting 

alloys 

Powder 

met 

alloys 

3xxx  AJ-Mn 

•3003  foil 

•3004  can  sheet 

2xxx  Ai-Cu-X 

•2x24  - aerospace  alloy: 
high  strength,  damage 
tolerant 

•2x19  weldable  higher 
temperature  alloy 

Sxxx  Al-Mn-Mg 

•5 1 82  - end  stock 

•Lighting  sheet 

•Architectural  sheet 

?xxx  Al-Mg-Zn-X 

•7x75  high  strength, 
corrosions  resistant 
•7x50  high  strength,  better 
corrosion  resistance,  less 
quench  sensitivity 

Microstructural  Changes 
During  Preheating 

Many  complex  simultaneous  reactions: 

• Dissolution  of  Mg2Si 

• Precipitation  of  Al12Mn3Si  dispersoids 

• Transformation  of  Al6[Fe,Mn]  to 
Al12[Fe,Mn]3Si  insoluble  constituents 

• Partial  dissolution  of  Al12Mn3Si 
dispersoids 

• Growth  of  insoluble  constituents  by  Mn 
diffusion 

• Growth  of  Al12Mn3Si  dispersoids 
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Hot  Reversing  Mill 


5 Stand  Hot  Continuous 
Mill 


• Single  stand/ 
multiple  pass 

• Ingot  thickness 
reduced  to  ~4" 
in 

multiple  passes 


• Both  head  and 
tail  sheared  at 
final  slab 
thickness 


Multistand/single 
pass  rolling 
Slab  thickness 
reduced  to  hot  mill 
gauge  of  -0.1" 
Coiled  at  exit 
temperatures  that 
allow 

recrystallization 


Continuous  Cold  Rolling 


Hot  Continuous  Mill : 
exit  - cooling  - coil  storage 


• Two  unwind 
stations 

• Coils  ends  are 
butt  welded 
together 

• Sheet 
unwound  into 
accumulator 

• Accumulator 
pays-off  into 
cold  mill  while 
next  coil  is 
welded 
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Continuous  Cold  Rolling 


Oiler  and  Slitter 


• Multistand/single 
pass  rolling 

• hmg  thickness 
reduced  to  finish 
gauge  of  <0.01" 

• Coiled  at  exit 
temperatures  that 
allow 

stabilization 


oil  applied  to  sheet 


slit  to  width,  cut  to  length 


Alloying  additions  and/or  impurities 


Example  research  goals 


No  element  is  completely  soluble  in  aluminum, 
most  solubilities  are  quite  limited 


• Preheat  and  solution  treatment 
practice  improvement  and 
simplification 

• Reduce  end-to-end  variability 

• Continuous  casting 

• Grain  structure  control 

• Development  of  aging  practices 

• Reduce  through  thickness 
variability 


• Major:  Ag,  Cu,  Li,  Mg,  Mn,  Si,  Zn 


•Minor:  Cr,  Zr,  V,Mn,  Fe,  Si 
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Peritectic  systems  that  need  work 


Almost  all  systems  have  eutectic  reactions 

Melting  reaction  data: 


eutectic 

pentectsc  fHHI 
~ml  solubility 


H 

n 

| 

He 

Li 

Be 

i 

I 

B 

G 

!0 

F 

|Ne 

Na 

Mg 

j j 

f 

m 

Si 

W: 

s 

a 

Ar 

K 

Ca 

Sc 

Ill 

Mn 

Fe 

Co 

:Ni 

SCu 

;Zn 

Ga 

Ge 

As 

Se 

Br 

!Kr 

Rb 

Sr 

Y 

i m 

■ 

Tc 

Ru 

Rh 

Pd 

M 

iCd 

In 

Sn 

Sb 

Te 

I 

iXe 

Cs 

Ba 

La 

Bf  Ta 

jw 

Re 

57 

Ii 

Pt 

Au 

iT 

?b 

Bi 

Po“ 

At 

rT 

Fr 

Ra 

;Ac 

i | 

\ 

! 

Ce  Pr 

Nd 

Pm 

Sm 

Eu 

Gd 

Dy 

Tb 

Ho 

Er 

£ 

Yb 

Lu 

Th  Pa 

TJ 

Np 

Pu 

Am 

Cm 

Bk 

!Cf 

Es 

Fm 

Md 

No 

Lw 

Temary+  systems 


Solvus  data 


;jkn 


own 


Idiscrepancy 


H 

lie 

Li 

IBe 

|b 

!C 

IN 

To~ 

If- 

!Ne  : 

;Na 

Mg 

|ai 

Si 

iP 

is 

ia 

iAr 

!K 

iCa 

Sc 

Ti 

|v 

Cr 

!Mn 

iFe 

!Co 

Nl 

Cu 

Zn 

^Ga 

Ge 

As 

Se 

Br 

:Kr 

Rb 

iSr 

rv 

Zr 

|Nb 

Mo 

|Tc 

■Ru 

jRh 

Pd 

Ag 

Cd 

iln 

iSn 

:Sb 

Te 

ii 

jXei 

Cs 

:Ba 

La 

|Ta 

■w 

{Re 

j°s 

k.. 

jPt 

■Au 

•Hg 

in 

:Pb 

iBi 

iPo 

|At 

iRn 

;Fr 

'Ra 

Ac 

;Ce 

iPr 

Nd 

iPm 

Sm 

iEu 

Gd 

Dy 

Tb 

:ho 

Er 

Tm 

Yb 

Lu 

Th 

|Pa 

u 

|Np 

Pu 

i Am 

Cm 

;Bk 

|Cf 

.Es 

Fm 

Md 

No 

Lw 

Al-Mg-Si-Sn 

Al-Mg-Si-Pb 
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for  free  machining  alloys 


Al-Cu-Mg-Ag 

-Si-Zn-... 

to  follow  up  on  Polmear’s  work 


•Al-Cu-Li  solvus  along  polythermal  sections 

•Al-Cu-Mg-Li 

•Al-Li-Zn 

•Al-Li-Si ... 

for  preheat  of  Li  hearing  alloys 
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Elements  of  very  low  solubility 


Phase  diagrams  to  support 
modeling  of  continuous  castin 


• Na 

• Ca 

• S 

• P 

•Represent  interactions  in  the  Hall  bath 
•Important  in  quest  for  inert  anode 


Metallurgical  characteristics  of  interest 
are  through  thickness  variation  of  - 

- Grain  structure 

- Cell  structure 

- Dendrite  arm  spacing 

- Constituent  particle  size,  number 

- Solute  distribution 


•Na,  Ca,  Li  removed  by  gas  fluxing  (furnace 
or  in-line)  or  salt  fluxing 

•Low  levels  may  affect  alloy  properties 


Thermodynamic  properties  to  support 
modeling  of  aging  processes 
(esp.  nucleation) 

• Free  energies  of  the  non-equilibrium 
crystal  structures 

• Interfacial  strain  and  energies  for  2nd 
phase  particles 


Modeling  requires  (in  addition  to  the  usual 
suspects)  ike  phase  diagram  under 
pressure 


Impurity  diffusion  data 


Q koai'mol 


Bulk  modulus(T) 

(Also  - heats  of  formation  of  phases  for 
which  we  have  no  measurements) 

Via  VASP  or  other  first  principles 
techniques 


Impurity’  diffusion  in  aluminum 

O 


$ ♦>* 
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Diffusion:  needs 

• Grain  boundary  diffusivities 

• Effect  of  non-equilibrium  vacancy  distributions 

Additional  data  to  support  modeling 

Example;  surface  quality  and  bending  of 
autobody  sheet  as  a function  of  hot  mill  lav-on 
temperature 

-If  MgjSi  growth  occurs  while  the  sheet  is  on  the 

•Liquid  A1  - oxide/nitride/carbide  interface  energies 

hot  mill, 

-then  it  sees  a non-equilibrium  vacancy 
distribution  that  varies  through-thickness. 

•Average  grain  size  distribution  and  average  grain 

-because  most  of  the  deformation  occurs  at  the 
surface  of  the  sheet. 

boundary  mobility  in  a real  commercial  material 

-Variation  of  difiusivity  gives  rise  to  variation  of 

Mg,Si  particles, 

-which  in  turn  gives  rise  to  through-thickness 
variation  of  properties. 

Additional  points  made  by  speaker: 

Elements  of  interest  can  be  grouped  in: 

• Major  elements;  Ag,  Cu,  Li,  Mg,  Si,  Mn,  Sn 

• Minor  elements;  Cr,  Zr,  V,  Mn,  Fe,  Si 

Aluminum  industry  has  two  basic  alloy  groups: 

• Non-heat  treatable  alloys  - e.g.,  lxxx,  3xxx,  5xxx  series 

• Heat  treatable  alloys  - high-strength  2xxx,  7xxx  series 

An  example  of  a modeling  opportunity  in  the  area  of  non-heat  treatable  alloys  is  prediction  of  metastable 
phases  formed  during  solidification.  Different  phases  form  depending  on,  for  example,  solidification 
rate.  The  distribution  of  these  phases  affects  the  (important)  surface  finish  properties. 

Modeling  opportunities  abound  in  the  area  of  heat  treatable  alloys.  To  optimize  homogenization  heat 
treats,  one  needs  standard  phase  diagrams  and  relatively  simple  diffusion  models,  but  the  solvus  and 
melting  temperatures  are  needed  with  very  good  precision  (i.e.,  < 2.5  °C  (5  °F)).  More  sophisticated 
precipitation  models  are  needed  to  design  new  alloys  with  optimized  combinations  of  strength,  fracture 
toughness  and  corrosion  resistance. 

Some  other  microstructure  modeling  needs  are  associated  with  grain  structure,  texture,  pressure 
dependence  of  solidification  microstructure. 

In  summary,  the  success  of  computational  materials  science  in  the  steel  industry  demonstrates  the 
feasibility  of  using  a computational  approach  to  design  industrial  processes.  It  remains  to  develop  the 
corresponding  databases  for  aluminum  alloys  and  to  codify  our  practical  expertise  in  terms  of 
quantitative  microstructure  models. 
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Applications  of  Computational  Thermodynamics  to  Virtual  Aluminum  Castings 
Ravi  Vijayaraghavan,  Ford  Motor  Company,  Dearborn,  MI 


Virtual  Aluminum  Castlaj 
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to  Virtual  Aluminum  Castings 


Ravi  Vijayaraghavan 
(Ford  Motor  Company) 


Workshop  on  Databases  for  Computational 
Thermodynamics  and  Diffusion  Modeling 
NIST,  Gaithersburg,  MD 


March  21, 2002 


Reduction  of  process  and  prod  uct 
develop  meat  time  in  the  manufacture 
of  eagfac  blocks  anti  cylinder  heads 


« John  Allison 

* Chris  Wolverton 

* Xinyan  Van  - ALCOA 

* Mei  Li 

* Jacob  Zindel 

* Shannon  Weakley  (University  of  Michigan) 
•Larry  Godlewski 
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PmJk’iion  of  MicrostrucTure  and  Mrcrosegregatmn  m 319  Aluminum  alloys 
- X.  Van,  R.Vljayaraghavan,  S-L  Chen  and  Y»  A-  Chang 
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First-principles  free  energies  for  Al*Cu  tr  phase 
Existing  commercial  empirical  thermodynamics  code 
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future  \i,x-ds 


* Robust  *n(l  wcU-teslitl  Mg  alloys  databases 

* Extensions  to  currently  available  Al  alloys  databases 

■ Robust  and  validated  mobility  database  in  1)1  O KA  for  At  - alloys 

* An  automated  optimizer  for  multicomponent  databases 

* Extension  of  tools  for  first-principles  thermodynamics  to 
multicomponent  systems 

* Extension  of  phasc-ftclcj  mieroslmctural  prediction  tools 


to  multicomponent  systems. 

• First  principles  information  on  liquid  alloys 


Additional  points  made  by  speaker: 

The  design  of  alloy  and  process  for  a cast  aluminum-alloy  engine  parts  (blocks,  heads)  by  a coupled 
thermodynamic-kinetic-property- service  model  is  maturing  in  a program  at  Ford.  A suite  of  codes  has 
been  developed  and  is  being  tested.  These  codes  include  a thermodynamic  calculator  and  database  to 
predict  thermodynamic  properties,  effect  of  process  variables,  solidification  micro  structure,  material 
properties  (residual  stresses)  and  service  behavior.  Structural  analysis  is  carried  out  with  the  software 
package  Abacus. 

A major  effort  in  alloy  improvement  is  also  underway  as  part  of  the  modeling  effort.  Alloy  319  is  being 
examined  for  potential  5%  to  10%  cost  saving  through  increase  in  Fe  content.  Work  is  well  underway  to 
model  and  predict  various  properties  as  a function  of  Fe  content  and  Fe  plus  second-element  effects.  The 
modeling  program  is  closely  matched  with  an  experimental  verification  program  on  these  alloys. 

Solution  and  aging  heat  treatment  is  an  area  of  great  opportunity  in  the  cast  aluminum  systems.  Work  is 
beginning  with  kinetic  modeling  of  phase  transformation  using  the  DICTRA  software  program. 
Although  empirical  modeling  of  properties  such  as  hardness  exists  for  these  processes,  it  is  hoped  that 
more  robust  models,  based  on  phenomenological  descriptions,  will  lead  to  new  alloy  and  process 
developments. 

Future  needs: 

• Mg  and  expanded  Al  diffusion  mobility  databases 

• Automated  optimizer  for  multicomponent  data  descriptions 

• Comparison  of  results  from  first  principles  with  thermodynamics  databases 

• Extension  of  first  principles  and  phase  field  modeling  for  multicomponent  systems 

• First  principle  modeling  of  the  liquid  phase 
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Use  of  Quasi-Thermodynamical  Models  for  Simulation  of  Epitaxy  of  III-V  Compound 
Semiconductors  in  Electronic  Industry 

Yuri  N.Makarov,  Semiconductor  Technology  Research  Inc.,  Richmond,  VA 


Use  of  quasi“tSiermodyiiamic 
models  for 

simulation  of  epitaxy  of  III-V 
compound  semiconductors  in 
electronic  industry 

Yu.N.  Makarov, 

STR  Inc.,  Richmond,  VA 


Outline 


introduction 

Modeling  of  epitaxial  growth 
Chemical  processes  in  MGVPE 

♦ Gas-phase  chemistry 

♦ Surface  chemistry 

Quasi-thermodynamic  model  of  surface  chemistry 
Application  example:  growth  of  group  Ill-nitrides 

♦ Evaporation  of  group  Ul-nitrides 
4 GaN  MOVPE 

♦ InGaN  and  AlGaN  MOVPE 


♦ Conclusions 


Activities: 


What  is  Semiconductor  Technology 
Research,  Inc.? 

Pre-history: 

1993-1996  Group  for  modeling  of  crystal  growth  and 
epitaxy  at  University  of  Erlangen-Numberg,  Germany 


Research  projects  in  USA  and  Europe 

Consulting  in  the  area  of  crystal  growth  and  epitaxy 
EUROPE:  AIXTRON,  Wacker,  Osram,  Nukem,  SiCrystal, 
Forschungszentrum  Julich,  FBH,  etc. 


1996-  establishing  CompuSoft  International  Ltd.  in 
Savonlinna,  Finland 

Today:  Consortium  for  modeling  of  crystal  growth  and 
epitaxy: 

STR  Inc.,  Richmond,  USA 
-STR  GmbH,  Erlangen,  Germany; 

-CompuSoft  Int.,  Savonlinna,  Finland; 


Metal  organic  vapor  phase  epitaxv  (MOVPE) 

of  compound  semiconductors  -why  ? 


USA:  Emcore,  Fox  Group,  CFD  Res.  Corp.,  Cornell  University, 
Cree,  TDI,  Kyma,  Crystal  IS,  NCSU,  etc. 

Supply  of  software  for  modeling  of  crystal  growth 
and  epitaxy: 

-“Virtual  Reactor”  software  for  modeling  of  sublimation  growth 
of  SiC,  AIN; 

-CVD-Module  for  modeling  of  MOVPE/CVD  processes  on  the 
basis  of  commercial  code  CFD-ACE 

-Simple  1 D code  for  modeling  of  MOVPE  of  lll-V’s  and  nitrides 
in  vertical  and  horizontal  reactors 


Metal  organic  vapor  phase  epitaxy  (MOVPE) 

of  compound  semiconductors  -why  ? 


High  Efficiency 
Solar  Cells 


Broadband/ 

Communications 


InP.  InGaP 


GaAs,  InGaP,  AlGa.As 
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Metal  organic  vapor  phase  epitaxy  (MOVPE) 

of  compound  semiconductors  - how  ? 

General  scheme  of  deposition  process 

Gas  flow  (III  asid  V group  precursors  and  carrier  gas) 


What  should  describe  an  adequate 
model  of  MOVPE  of  lil-V  compound  ? 
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Metal  organic  vapor  phase  epitaxy  (MOVPE) 

of  compound  semiconductors  - how  ? 

Basic  reactor  types 


Horizontal  reactor 


Planetary  reactor 


Vertical  reactor 


Examples  of  flow  and  heat  transfer 
modeling 


The  $mmr  tspremm.  mt  of 


Susceptor 


Flow,  heat  transfer,  and  species  transport 
in  epitaxial 

reactors  can  be  modeled  using 
well-established 

Computational  Fluid  Dynamic  software. 


Gas-phase  reaction  mechanism: 
unimoleeuiar  decomposition  of  TMGa 


•§° 


<H 


kcal/mol 


■§°CH'  J° 

^ 77  ^ hrnl 


CH, 


35.4  kcal/mol  ^ 77 5 kcaUmol 

DMGa  MMGa  Ga 


°§0 


This  reaction  mechanism  is  typical  for  all  group  IH  precursors: 
Trimethylgallium  (TMGa),  trimethylaluminum  (TMA1), 
Trimethylindium  (TMln). 


Gas-phase  chemistry  in  AlGaN  growth 


TMA!  'jjtf.S*0 
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1.  Condensation  erf  adduct  TMAi:NfF,on  reactor  Inlet  and  cold  walls 
2.  Formation  of  oligomers  and  AIN  particles 
(based  on  T.G,  Mihopotios  et  aL,  ./.  Cry  si.  Growth,,  195, 733  (1998)) 
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Chemistry  in  MQVPE 


Full  kinetic  mode!  for  surface 
processes 


• Gas  phase  chemical  reactions  may  affect  the  deposition  rate 
and  uniformity  by  several  reasons:  a)  different  diffusion  rates 
of  precursors  and  reaction  products;  b ) parasitic  reaction 
leading  to  formation  of  particles  and  losses  of  material; 

• Rate  constants  of  reactions  are  usually  taken  from  experiments 
on  precursor  pyrolysis  in  flow  tube  reactors.  The  estimation  of 
molecular  structure  and  bond  energies  of  precursors  using  quantum 
chemistry  methods  may  be  very  useful  to  evaluate  the  validity1 

of  commonly  accepted  rate  parameters. 

• The  main  “bottlneck”  in  MOVPE  modeling  is  the  description 
of  surface  chemical  processes. 


Drawback: 

* lack  of  necessary  informa- 
tion on  mechanisms  and 
reaction  rate  constants 
requires  to  fit  parameters 


_.  A total  of  60  species  and  232  reactions  are 
considered  in  the  gas  together  with  19  species 
and  115  processes  at  the  surface  - 

51  AsH+  V <=>  AsH(a) 

52  AsH,  + V » AsHj  (a) 

53  AsH3  + V <=>  AsH3  (a) 

54  H + V<=>  H(a) 

55  CHj  +V<=>  Cl^ia) 

56  AsH  (a)  + H(a)  <=>  AsH,  +2V 

57  AsH,  (a)  + H (a)  <=>  AsH3  + V 

58  AsH,  (a)  + H (a)  <=>  AsH,  (a)  + V 


Kev  features  of  the  q uasi -thermod vnami c model 

chemical  processes  at  the  growing  surface 


of 


Advantages  of  quasi-thermodynamic 
mode!  for  surface  processes 


Basic  assumptions: 

1 . The  rates  of  atomic  incorporation  into  the  crystal  and  of  the  crystal 
decomposition  are  much  higher  than  their  net  difference 
(quasi-equilibrium  between  the  adsorption  layer  and  the  crystal  bulk). 

2.  Kinetic  effects  at  the  stage  of  adsorption/desorption  are  accounted  for 
by  the  sticking/evaporation  coefficients  of  individual  species. 

3.  Kinetic  effects  at  the  stage  of  species  transport  are  accounted  for 
by  solving  transport  equations  (Navier-Stokes  equations) 

Capabilities  important  for  modeling,  of  semi,con<ju;tQ/^o.M.h: 

- the  model  allows  the  calculation  of  species  desorption  rates  from  thermo- 
chemical  properties  of  gaseous  end  solid  substances; 

- the  accounts  of  the  clasSc  strain  effect  via  a contribution  to  the  Gibbs 
energy  of  solid  phase 


Group  Ill-nitrides  - metastable  materials 

with  unusual  properties 


Formation  of  kinetic  barrier  in  adsorption 

and  desorption  of  nitrogen 


• Nitrides  do  not  react  chemically  with  nitrogen; 
nitrogen  activation  is  required  for  growth. 

• GaN  evaporates  congruently  under  vacuum  conditions: 
higher  pressures  change  evaporation  mechanism. 

• Nitrides  are  stable  at  temperatures  > 300  °C,  whereas 
the  thermodynamic  estimations  predict  crystal 
decomposition 

• Liquid  phase  has  catalytic  effect  on  crystal 
decomposition 


54 


Workshop  Presentations 


Kinetics  of  nitrogen  evaporation  and 
condensation  - reason  for  the  mctastabilitv 


thermodynamic  1 

s'v'vx.  / 1 

experiment  : 

' AIN 

«v  i 
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Bcparimsnt:  L.H,  Drager,  et  all, 
XPhys.Chem.  68  {1963}  1556 


7.0  7.5  8.0  8.5 

lOVTemperature  (K) 

ZA,  Fvterssf,  A.W. 

Searcy,  j.Cfiem.Phys.  42  {1865}  4223 


IVIOVPE  growth  of  GaN 
(comparison  with  experiment) 
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Experiment:  C.H.  Chen  et  al. 
J.EIectrou.Mat.  25  (1996)  1004 
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2I>  modeling  of  InGaN  MOVPB  in  AIX  29#  RF  reacton 

comparison  to  eroerimeatel  data 
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Experimental  data  from:  M.  Schwambera  ei  al,  J.  Crystal  Growth,  203,  340  (1999) 

Operating  conditions:  total  flow  rate  4.5  SLM,  pressure  50  mbar, 
V/III  ratio  ~ 3000;  carrier  gas  nitrogen 


Surface  chemical  processes  daring  GaN  growth 


nh3  n g. 


^ ^ ( A 

@OfQQ  @©  OO  0 0 § cjocfo 


Substrate 


Rate- limiting  processes 

• Low  temperatures:  sites  for  Ga  adsorption  are  blocked  by  methyl  radicals 

• Intermediate  temperatures:  transport  of  Ga  containing  species 

to  the  growing  surface 

• High  temperatures:  Ga  desorption 

There  is  the  strong  experimental  trend  to  move  towards  the  high  GaN 
growth  temperatures  because  of  die  better  quality  of  the  material.  Accurate 
description  of  gaDiuin  desorption  front  the  growing  surface  is  necessary  for 
predictive  modeling. 


Surface  chemical  processes  during  InGaN  growth 


, NH,  N, 


( 


ft  MMIh  _ |n 
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Substrate 


Primary  reaction  pathways 

• Adsorption  of  gaseous  species 

• Ammonia  cracking  on  the  surface  giving  rise  to  N,  and  H2 

• Indium,  gallium,  and  nitrogen  incorporation  into  the  crystal 

Desorption  of  indium:  the  limiting  pr< 

iamqzo^ 


2P  modeling  of  InGaN  MO  WE  lit  AIX  200  RF  reactor; 

comparison  to  experimental  data 


® Experiment 

— Computations 

/.... 

0 10  20 

30  40 

Experimental  data  from:  MRS  Internet 
J.  Nitride  Semicond.  Res.,  4,  5 (1999) 


Operating  conditions: 

temperature  750  °C, 
total  flow  rate  4.5  SLM, 
pressure  50  mbar, 
V/m  ratio  20000^30000, 
carrier  gas  nitrogen 


In  content  in  the  gas  phase,  % 


lire  model  provides  a good  quantitative  reproduction  of  the  experimental 
is  very  important. 
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Effects  of  strain  and  T: 


Vertical  RDR  , 

Vg 

2Vg 

4Vg 

: 

700  800  900  1000 

Temperature,  °C 


700  800  900  1000 

Temperature,  °C 


I to.  content  decreases  gradually  f 
with  temperature  due  to  indium 
desorption. 

Elastic  straws  enhances 
indium  desorption 


Indium  desorption  is  suppressed 
by  increasing  growth  rate 


Effects  of  surface  chemistry  in  AlGaN  MOVPE; 

effect  of  growth  temperature 


900  1000  1100  1200 

Growth  temperature,  °C 


Experiment: 

S.  Ketter  el  al, 
Appl  Phys.y 
86,  5850(1999} 


At  high  feruperatam, 
desorption  of  gallium  becomes 
significant  and  results  in 
depletion  of  the  growing  layer 
with  galisiu? 
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X' 

Double  PL  peaks  ar< 
observed 

in  the  temperature  range 
700-350cC  where 
die  strain  effect  on  InGaN 
composition  is  mostly 


500  600  700  800  900 

Temperature,  °C 


The  appearance  of  In-rich  arid  In -depleted  zones  is  related  to 
co-existence  of  strained  and  unstrained  fe.Ga.N  Islands 

rather  than  to  phase  separation 


Effects  of  surface  chemistry’  in  AIGaN  MO  WE; 

effect  of  the  growth  rate 
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80 

Extrcrimenf: 

S.  Keller  el  al, 
J.  Appl.  Pliys., 
86,5850(1999 } 


Intensive  gallium  desorption  may  be  suppressed  by  increasing  total 
group-ill  flow  (grow  th  rate). 


Particle  formation  during  chemical  vapor  deposition 
ofSiC 

(A.N.  Vorob’ev  et  al..  Materials  Science  and  Engineering  B61-62, 
172(1999)) 

Irradiant  layer  observed  over  substrate  in  experiment  at 
' Siemens  (.Erlangen,  Germany) 


Particle  formation  during  chemical  vapor  deposition 
of  Si  (cooperation  with  NIST  - R.  Davis,  J.  Maslar) 


Layer  of  silicon  particles 
(modeling), 

A.N.  Vorob’ev  et  al., 
Electrochem,  Soc.  Proc.. 
2001-13.462  (2001) 


Layer  of  silicon  particles 
(experiment), 

R.W.  Davis  et  al„ 

AJ.P  Conference 
Proceedings, 

550.  292  GOOD 
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Summary 

♦ The  main  difficulties  in  MOVPE  modeling  are  related  to 
development  of  predictive  surface  chemistry  models. 

♦ Quasi-thermodynamic  approach  allowing  description  of 
the  surface  chemical  processes  and  growth  on  basis  of 
thermochemical  data  of  gaseous  and  solid  species, 
represents  the  good  alternative  to  detailed  multistep 
kinetic  models. 

♦ Quasi-thermodynamic  models  of  nitride  cry  stal  growth  and 
epitaxy  require  reliable  thermodynamic  properties  of  solid 
phases  and  vapor  species 


Wliat  we  need? 

♦ Reliable  thermodynamic  properties  of  GaN,  AIN,  lisN, 
AJGaN,  InGaN,  AlGahiN,  AUnN 

♦ Thermodynamic  properties  of  vapor  species  and  products 
of  decomposition 

♦ Thermodynamic  properties  of  different  adducts  between 
metalorganic  species  and  ammonia 

♦ Kinetic  parameters  of  rate-limiting  processes  (diffusion 
coefficients,  rates  of  gas  phase  chemical  reactions,  sticking 
probabilities,  etc.) 

♦ Experiments  and  models  of  adduct  and  particles  formation 
in  the  gas  phase  during  CVD 


Use  of  Valence  Force  Field  (VFF)  modeling  for 
estimating  thermodynamic  properties  of  multi” 
component  compounds 


Use  of  Valence  Force  Field  (VFF)  modeling  for 
estimating  thermodynamic  properties  of  multi” 
component  compounds 


lattice  constant 


d - lattice  constant 

A 


Elastic  energy  of 
lattice  distortion  is 
described  in  terms  of 
stretching  aid  bending 
of  tetrahedral  covalent 
bonds 


Elastic  energy  of 
lattice  distortion  is 
described  in  terms  of 
stretching  and  bending 
of  tetrahedral  covalent 
bonds 


-Cation 
0 - Anion 


Enthalpy  of  mixing  for  ternary 
compounds 


Comparison  with  experimental  data 

Interaction  parameter  W is  defined  by  the 
expression 


II 


: W-x(l-x) 


0.0  0.2  0.4  0.6  0.8  1.0 

InAs  mole  fraction  in  InGaAs 


t 


Agrees  well  with 
available  data 


Compound 

V EE  -modeling 
(kj  Mole) 

Experiment 
(kj  Mole) 

InGaAs 

12.985 

11.76-12.56 

InGaP 

16.112 

13.36—1 8.45 

InGaN 

30.464 

□ 

GaAsP 

4.283 

4.19 
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Contributors: 

R.  Talalaev,  E.  Yakovlev,  I.N.  Przhevalskii,  S.Yu. 
Karpov.  A.N.  Vorob'ev,  A.V.  Lobanova,  N. 
Podol'skaya  and  I.Zmakin 


Additional  points  made  by  speaker: 

Group  HI  nitrides  offer  the  opportunity  to  develop  super-bright  LED  lights  with  a potential  market  in  the 
$ 20  B range. 

Metal  organic  vapor  phase  epitaxy  (MOVPE)  is  the  preferred  method  of  fabrication  for  these 
semiconductors.  Current  drawbacks  to  development  are  associated  with  a lack  of  good  process  models 
and  empirical  data.  The  quasi-thermodynamic  model  presented  here  offers  a way  to  overcome  this  lack 
of  data. 

The  semiconductor  nitrides  of  greatest  interest  are  GaN,  AlGaN,  and  InGaN.  Progress  has  been  made 
using  a quasi-thermodynamic  approach  to  model  the  epitaxial  growth  process  at  the  growing  surface. 

The  ability  to  accurately  model  the  growth  process  is  complicated  by  the  process  variables  of  strain, 
temperature,  pressure,  surface  chemistry,  and  kinetics.  The  quasi-thermodynamic  model  takes  the  host 
of  process  variables  into  account  and  has  been  applied  to  group  in  nitrides.  Results  are  encouraging  to 
date.  Opportunities  exist  in  the  need  for  accurate  surface  chemistry  models,  data  on  thermodynamic 
properties  of  EH-nitrides,  and  kinetic  parameters  for  various  reactions. 
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NSF  Supported  Computational  Education  Program  at  Penn  State 
Zi-Kui  Liu,  Penn  State  University,  University  Park,  PA 


An  Intearated  Education  Program 
on  Thermodynamics,  Kinetics, 
and  Materials  Design 

Supported  by  NSF  DMR-0073836 

Zi-Kui  Liu*  (liu@matse.psu.edu) 

Long-Qing  Chen,  Karl  Spear  and  Carlee  Allison 


Objective 

Improve  the  student  learning  experience  and 
educator  teaching  experience  on  two  of  the 
core  contponetiis  in  the  curriculum  of 
materials  science  and  engineering,  i.e. 
ThermcxJynamtcs  and  Kinetics, 
emphasizing  on  graduate  study. 


Department  of  Materials  Science  and  Engineering 
The  Pennsylvania  State  University 


Materials  Science  and  Engineering 


PERFORMANCE 


f rom  hup:  twvv.erc4mse.org  i rosn  hup;-  vvww.mv  n.edu 


C onven  t i ona  1 S itualion 

* Common  stereotype:  Therm odynatmes.  and 
kinetics  are  problematical  to  learn  and  difficult  if 
not  impossible  to  apply  in  the  real  world. 

» Traditionally,  the  fields  of  materials  science  and 
engineering  have  been  predominantly  on 
processing  of  materials,  establishing  structure- 
property  relations,  anti  measuring  properties. 


The  core  is  not  reflected. 


Current  Graduate  Curriculum 


Materials  Science  and  Engineering 


Research 
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Development  in  Recent 

* Computational  Thermodynamics 

- Coupling  of  phase  equilibrium  and  thermochemistry 
(CALPHAD) 

~~  Firsl-principles  calculations 
» Phase  transfonnation  simulation 
— Sharp  interface 

- Diffuse  interface  (phase-field) 

* System  materials  design  based  on  the  above 

- Northwestern  University  and  Questek  Innovation 


Faculty  Expertise  at  Penn  State 

The  Project 

* MSE:  34  faculty  members.  J 50  graduate  students. 

♦ Time  Period:  09/0 1/00-08/3 ! /03. 

* Karl  Spear 

- Ttemuthnands'  modeling  of  ceramic* 

* Financial  support:  NSF,  University 

- lead;  graduate  thc*mi<>sly«»mk-s 

matching  and  in-kind. 

* Lona-Otna  Chen 

- Phase  Ifcki  stouMen 

♦ Build  a new  PC  lab  in  the  department 

~ Teach  undergraduate  therniadynamies.  graduate  felaetics 

* Zi-Kut  Liu 

• Web  site:  http://nsfedu.metsce.psu.edu 

- rhetmudynanuc  modeling  *4  alloys 
feadt  undergraduate  phase  SfanstonnatioB 

- Knowledge  of  system  material*  design 

3<> 

Modified  Undergraduate  Courses 

* MatSE40  f : 1 hemodynamics  of  Materials 

— introduction  of eoiapuiaiksnal  thermodynamics:  3P-T 
diagrams  of  pure  element  . 

* Metal  405:  Phase  Transformations  in  Alloys 

— Calculation  of  phase  diagrams  and  equilibrium  phase 
fractions 

* Metal  435:  Metallurgy  Lab  !l 

l.xperijneiiuif  and  calculated  phase  diagrams 

* Metal  436:  Metallurgy  Lab  111 

A computational  tab  on  avoiding  dogging  in  a 
continuous  casting  process 


The  Computer  Lab 
with  20  PCs  for  students 


mmm 


An  Integrated  Education  Program  on 
Thertnodvnaraics,  Kinetics,  am!  Materials  Design 
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Modi  lied  G rad uate  C ou rses 


New  Courses 


* MatSCSOI:  Thermodynamics  of  Materials 

— Addition  of  one  session  evesy  week  on  computational 
thenitodyoasnics 

- Homework  on  computational  activities 

- Exams  based  on  the  output  of  calculations 

• MatSC503:  Kinetics  of  Materials 

introduction  of  mobility  database 

- Homework  oo  diffusion  simtslMfem  ffikmt) 
introduction  of  Montc-Carto.  microscopic  diffusion  and 
phase-field  simulation  tools. 


* MatSC59?C:  Intmditction  trf Computational 
Thermodynamics 

i itermoilj  uatnic  model  ing,  Database  development 

» VlatSC 597k:  Computational  Materials  Science  11: 
Continuum,  Mcsocalc  Simulations 
Dscfns,  Monte-Carlo.  Phase-field 

* MatSC  59SA:  System  Materials  Design 

- Engineering  design  principles.  Materials  design 
procedure.  Design  projects 


U 


Example  of  MatSE401 


Metal  436:  Computational  Lab 


Example  of  Metal405 


Homework  of  MatSCSOI 


* A manufacturer  wanted  to  increase  the  Cr 
content  of  a material  from  18  to  25  weight 
per  cent.  However,  clogging  starts  to  occur 
during  the  continuous  casting  of  this 
material  because  solid  Cr,0.  is  formed.  By 
calculating  the  equilibria  in  the  sted/slag 
system,  a simple  correction  could  be  found 
through  the  adjustment  of  the  steel  melt 
composition  or  temperature. 


♦ Answer  the  following  questions  using  Fact  Sage. 

- Can  a Tsii(s)  crucible  he  used  to  contain  te305  with  no 
reaction  at  1000K?  As  IbtiOK?  If  a reaction  occurs, 
what  ate  the  solid  products? 

- Can  a Ni{s>  crucible  be  used  to  contain  I'c.O,  with  no 
reaction  a!  I300K?  If  a reaction  occurs,  what  arc  the 
solid  products? 

Can  a Nits)  crucible  be  used  to  contain  1 cO  with  no 
reaction  as  1 BOOK?  li  st  reaction  occurs,  what:  arc  the 
solid  products? 


IS 
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Homework  of  MatSCSO  1 


Homework  of  MatSC503 


A tknvusg.  1UM  ! • gas  mixture  Is  «|uiiihntSed  with  n krrgs> 
q®tttiiy  of  8 Tit  xlTiOts)  mixture  at  1900K  before  ii 
passes  over  small  separated,  samples  of  V,  Mu  Or,  Zn. 
Cm,  t o,  and  -Ni  heated  at  I HOOK  Use  SM  Is  passed  over 
these  samples  until  they  finally  equilibrate  with  the 
i 1,0  I K ratio  produced  by  equilibnuirat  with  tire  HJ0OK 
1 ilMl  iOts)  mixture 

- Wist  tv  tit  the  dscmicaS  xtskc  of  the  islsvvc  meet!*  he  when  the 


1-abd  the  phase  radars?,  its  the  attached  tmt  S e-St-C  rsrsptdh  diagrams 
wish  0.45%<  ■ ami  tt.4«C.  wsp«#vef> . ! he  ssssmte  irdi cate  ste 
phase  Ixiustdaries  wish  site  mown  phases;  sis  equstsbrhiai  with  each 
taker. 

Sunnlaie  the  an  hilf  tltfltssiim rtf  the  K«*Si-C  diffusion  couple, 
t. xpkitii.  t.lsei«»dytiamh-,iSly  .bsP  tasselkaisy.  why  cat  Mi  disTyses  tram 
the  tow  easlrson  camsaitisttan  reg sea  (A)  so  the  high  carters 
dmwratfon  region  sits 

Simulate  the  time  needed  to  h(«yiswem/c  the  ea *hm  eoneerKfassost  m 

nuk. 


utece  ticscrihetS  equilibrium  is  reached  its  the  sy  stem? 

What  « i ;?  the  lirsai  oxy  gen  potential  he  at  equilibrium  in  the 
KWOK  furnace 


Make,  asasthar  simutetton  at  a higher  tempemu#*  where  the  two  alloys 
are  still  In  the  slag le  tec  phase  region  and  find  out  the  time  needed  to 
homegetstee  tl»  embem  ®accmtsahm. 


Software  in  MatSC503 


Therm odynam ic  M odeling 


MatSC59?C:  Computational  I hennodyrtamlcs 

» iaimtteton 

Review  th*MV5>#tmais  pinslgfe 

Rets  lie.  ehemieoi  iesicihes,  sstni  pistsve  ybyts.ssv  is  v t besmear  seism  s- 
phiicipk's 

* Hienrasd}  ouniie  Modeling  of  Pose  Sssbskuicvv  tJuiary  ant! 
XteticesnpoeetH  Systems 

AtetfeiaistOsiefeg 

ivtetetsssg:  -sf  im&m  satailw*  fSawtes  m$  sslttettee  soiBSissrs  phases 

* Camjnsteiiomil  1 echmque*  grid  Applications 

J.M«ta>e  of  gate  selrshsswes 
Dtshsisase  of  feitiar-y  systems 
i Cf.sksse  ( <s'e-m:;f>  smsUkssisssvisscnf  systesWS 
Case  studies  oftcai  prsfctesns 

* Mmiprajeet  of  the  tmm 


MatSC597K:  Continuum.  Mesocale  Simulations 


| Thefftiochetttkai  tltda: 

Pltase  equllibrktti!  data: 

1 enthalpy,  entrapy.  hear 

liquidity,  srtltdus.  phase 

1 capacity,  activity. 

boundary. 

SfH'fW  lsv«e  ftf  ( M HI  \i>  im  € tuttppWwwi  1 state 
Es.TH.EEM*  FACT*  MTfMTA,  PAMIMT, 
THiRilCiCAlTAfHFHMOSCiflT:  » 


* Or,  Uu 

fKv-  k<v  \hi;i\n<x}yvw<k  imj  kir^ik 

MtOx'tno.-  o3  timOwy 

* Dr.  D*h Roy 
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IWSd  wjvcs^r^i raa&:jVsfft>a3$i3S! 

* Da  . Oim 

Kincliar  ;trl<>  xjmu  kmof r <sc<fer<K£  ml  hvr«srali«>n 

Mkmeopfc  <8$ i*m  pix^m 

3 i »ipNaf'>J.  piBaim  md  S m4m 

i af»-»4l»iUi8-4^84tKwv,  f>1w.v  s <$&&&&.  prvviptinlc  mxKttiiK 
C.Vted'Rvs  Mid  iippiivi&  fliChiC:  mlxvivsH 
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MatS€$9A:  System  Materials  Design 


MatSC598A:  Project  One: 


* Principles  of  Design 

- Materials  selection 

- Materials  as  a system 

Sy  stem  approach  to  design 
— Cost  and  time  management 

* Design  Tools 

— ThmmKiynansic  design  of  unemstruentre 
— Computational  thermodynamics  (Thenno-Cale) 
— Computational  Kinetics  (Dictat) 

* in-Cfass  Example 

* Team:  Project 


MatSC598A:  Project  Two: 

i oday  s commercial  high  strength  loo  alloys  (HSLA)  have 
typical  yield  strength  of  80ksi  (about  56kgi?mm-  or  550MPa). 
The  figure  hekm  shows  the  development  of  f IS!  A in  the  past 
half  century  With  the  ocean  exploration  of  gas  and  petroleum 
going  to  deeper  and  deeper  water  and  the  increase  of  pipeline 
pressure,  a.  yield  strefigt.lt  as  high  as  120ksi  is  desirable,  Design 


Research  data  demonstrate  that  HO  to  of  air  pollutant 
due  to  exhaust  gas  of  a vehicle  is  generated  w (thin  200 
seconds  between  a vehicle  start  and  the  time  the  catalyst 
reaches  its  optimal  operating  temperature  (around  540  C). 
This  is  particularly  a problem  for  vehicles  used  for 
delivery  as  they  are  sw Itched  on  and  op  often  everyday. 
You  are  asked  to  make  a materissl-adMcd  design  to 
develop  an  optimal  solution  for  this  ease  using  the  system 
materials  design  approach. 


Evaluation  Procedures 

* Focus  groups  to  gather  student  feedback 

* Surveys  of  attitudes  and  knowledge 

* Classroom  observations 


such  an  alloy  using  the  System  Materials  Design  approach  with 
other  properties  balanced  based  on  your  application. 


Analysis  of  student  presentations  and 
reports  through  use  of  rubrics 


Survey  Results 


2002  Summer  Workshop 


I he  majority  of  students  indicated  that  course 
packets  containing  notes,  problems,  and  guides  for 
using  the  technology  would  aid  their  learning. 


May  20-21  Computational  Thermodynamics 

May  22-23  Processing  Simulation  and  Dietra,  and 
Introduction  of  Phase-Field  Simulation 

.May  24-25  Thermodynamic  and  Kinetic  database 
development 

http ://  n sfedu . m etsce.psu.edu 
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Vision 


Additional  points  made  by  speaker: 

Examples  of  course  problems  and  materials: 

• Clogging  in  continuous  casting 

• Phase  diagram  and  free  energy  curves 

• Thermodynamic  modeling 
Materials  design  projects: 

• Catalyst  design 

• 830  MPa  (120  ksi)  HSLA  (high-strength  low-alloy)  steel 
Questions/comments  from  workshop  attendees: 

Are  course  materials  available? 

• Workshops  are  free  to  participants  from  education 

• Materials  can  be  distributes  by  CD  or  e-mail 

• Web  site 

How  many  students  are  involved? 

• 30,  mostly  PhD  students 
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Thermodynamic  Measurements 

Philip  Nash,  Illinois  Institute  of  Technology,  Chicago,  IL 


ILLINOIS  INSTITUTE  ▼ 

OF  TECHNOLOGY 

Thermodynamic  Measurements 

Philip  Nash,  Hsm-Ning  Suh 

Thermal  Processing  Technology  Center 

Illinois  Institute  of  Technology 

10  W 32nd  St,  Chicago.  IL  60616 

Ackrtott&cigmente 

Oie  Kleppa.  K.  Ftaner.  ft  heddy,  H Jacohsan 

i -HsKf  mmMwt&s 

Outline 

• introduction 

• Expenmental  Methods 

* Experimental  Thermodynamic  Activities 

* Summary 

introduction 

• Aiioy  and  process  development  requires  a thorough 
knowledge  of  phase  equilibria  and  thermodynamics 
micFOStruciurai  development,  nan-isothermal 
processing,  distortion,  quenching 

• Available  experimental  data  is  limited  particularly  for 
ternary  and  higher  order  systems 

• Calphad  technique  otters  a method  of  determining 
phase  equilibria  and  thermodynamics  in  an  alloy 
system  in  much  shorter  time  than  by  experimental 
methods 

Introduction 

• Thermodynamic  modelling  is  dependent  on  accuracy 
of  input  data 

• First  principles  calculations  offer  another  method  of 
determining  phase  equiiibna 

* Some  thermodynamic  and  phase  equilibria  data  are 
essential  for  accurate  modeling 

♦ Accuracy  of  experimental  data  is  critical  to  reliable 
modeling 

9M&S&  i S * S & 
or- 

Data  Needed  for  Modeling 

Introduction 

* Gibbs  energies  of  individual  phases  as  a function  of 

composition  and  temperature 

- The  results  of  different  types  of  thermodynamic 

* Heat  capacities 

measurements  are  interrelated  through  the  equations 

* Activities 

of  thermodynamics 

♦ Enthalpies  of  formation,  mixing,  transformation 

• Entropies  of  formation,  mixing,  transformation 

* Different  types  of  measurements  complement  one 

* Phase  equilibria  (temperatures,  compositions) 

another 

* Physical  properties  (density,  thermal  expansion,  etc.) 
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Experimental  Methods 


Calorimetry 

Direct  reaction 
Soiaion 

Solaie-sotveni  Sr op 
Differentia!  scanning 
Vapor  pressure 

Gas  phase  equiUfcria 
Rate  of  evaporafioo 
EMF 

Other  methods 


Experimental 


The  heats  of  formation  are  determined  using  a high  temperature 
reaction  calorimeter  with  a typical  accuracy  of±ikJ/mole 

The  measurements  are  generally  made  with  the  calorimeter  set  at  1473&2K, 
using  an  argon  atmosphere 

The  calorimeter  is  calibrated  using  pure  copper 

Samples  are  produced  by  mixing  elemental  powders  in  a mortar  in  the 
required  molar  ratio  and  pressing  them  into  a small  pellet 

Typical  sample  weight  is  about  Uhl  mg 
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Kieppa  Calorimeter 
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AI*Ni-Y  Isothermal  section  [3] 


Heat  of  Formation  of  Ni  Al, 
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Results 

ttv  using  direct  synthesis,  the  standard  enthalpy  of  formation,  AtttJ®1K , is 
eah'id;tfvd  from: 

»AHs,  298KI  * hNifs,  2<>ftK)  ♦ cYfs,  2>MSK)=  A»aN\Vr{  147$  K)  Afik„tte„  f ! ) 
AI.NVUS.  2W  KJ « AiaNibV;.fl473  k)  AHfka, (2) 

from  reaction  (I)  ami  (2)  we  get 

aAi{*,  24>KK)  4 b\i|s*  29»K)  4 eVfv  29.Sk)  - 298  k) 

t he  vkutihn  tf  enthalpy  of  fmm.ilmn  is  thit*  uhmiiud. 

All  ^ sb  AH  * All 

_il  I f 1 ttfcjN-ri:**  ^ * i iras  < :*«*«« 
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, are  molar  enthalpy  changes  for  reaction  ( l ) and  (2|. 
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Enthalpies  of  Formation  of  Al-Ni-Y  intermetallics 

Ni 
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Advantages  and  Limitations  of 
Direct  Synthesis  Calorimetry 
Advantages 

Drily  two  heat  effect?,  need  he  measured 
Absolute  temperature  measurement  not  needed 
High  accuracy  (typically  * IfUfmote} 

Experiment  directly  yields  enthalpy  o?  formation  as  29SK 

Limitations 

Some  systems  give  incomplete  reactions,  particularly  very 

fts§b  melting  paint  compounds 

Temperature  drift 

fiedttefors  in  sensitivity  over  time 


Heat  of  Mixing 


Heal  of  Mixing 


<«■ 


Schematic  of  calorimeter  set-up 
for  mixing  measurements.  Solute 
reaches  thermal  equilibrium  with 
solvent  before  being  dropped 
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tfNY-S  = NC<‘  ::v;  i 


Differential  Scanning  Calorimetry 

the  temperature  difference  controls  Hie  electrical  power  to  the  sample  anil 
reference  iu  order  to  maintain  the  same  temperature.  The  peak  area  directly 
rorres|Hinil«  to  the  heat  consumed  or  produced  In  the  sample. 

Temperatures  are  measured  in  thin  plates  in  contact  with  crucibles,  therein 
mt-WH  to;  the  difference  in  heat  How  from  them.  St"nal  proportional  to  (he 
difference  in  heat  capacities  between  the  sample  and  reference  and  thus  the 
instrument  works  as  DSC  , 


r^ww. 


^^ismmsu  **«*«*?** 


• ptets 


1600C  DSC-TGA 


Simultaneous  DSC-TGA  to 
1GQQC.  DTA-TGA  to  1750C  in 
Oxidizing,  reducing  or  inert 
atmosphere 


EMF  Measurements 


Applications  of  Calorimetry  and  Sources  of  Error 

* Integral  and  partial  enthalpies  of  formation,  mixing  and 
transformation, 

♦ Heat  capacity-error  for  scanning  DSC  -3%,  for  step 
DSC  ~2%,  for  drop  calorimetry  -5% 

• Errors  can  arise  from 

mass  and  temperature  measurements 

incomplete  reactions 

reactions  with  crucible  or  atmosphere 

impurities 

eaiisraitai  errors 

beat  effect  errors 


AG  = — nr  t 

i M I measured  as  a ftmetion  at'  temperature  tor  fixed  composition 

E - a + bT 

AG.=RT\na.=AH--TASl 
= -nF{a  + bT ) 


Application  of  Gibbs-Duhem  permits  calculation  of  the  activity 
of  the  second  component  and  hence  integral  Gibbs  energy 


Experimental  Setup 


: saasitMio* 


l ...i 


R.  G F&cioy.  A.  M.  Yabya  avxi  L.  &r&#i>r,  'TiiWin&iynQmrc  Pfop&tfc&GfTi' M 
int&r'()Zl&:.:cs\  Juvrml  of  Ai^ys  and  Compounds.  Y <il  321,  up  223-22?.  2031 


TsmperaiUire  (K) 

Measured  EMF  of  T^Al  intermetalHc  vs  T 

8.(7.  Xftfity  et  af.  Journal  AXms  aat  Ompataufs  W (Ml} 


Workshop  Presentations 


69 


Activ ily  coefficient  vs.  Xv  at  S5(t  K 

R G ft&Zdf  st  ol.  Jauntat  cf  Jiitrcs  and  CfmptnaxZs  Ml  (2*HH t 223 -22? 


Facilities  for  Experimental  Studies  of 
Alloy  Thermodynamics— 
Inputs  to  Calphad  Phase  Models 


Applications  and  Sources  of  Error 

* Activities,  partial  and  integral  Gibbs  energies, 
enthalpies  and  entropies 

• Major  limitation  is  in  maximum  temperature 

♦ Errors  can  arise  from 

mass  snd  temper^ere  measurements 
reactions  with  atmosphere 
impurities 

multiple  valance  states 
thermal  gradients  in  cell 


Knudsen  Cell  Technique 


* Evan  Copland.  Nathan  Jacobson — Materials  Kvfeon , RASA  Glenn 
Research  Center,  Cleveland,  OH 

• High  Temperature  Knudsen  Cell  Mass  Spectrometer  tor  thermodynarmc 
act-Vity  measurements 


Btefibrate  aStoy/va cor  in  smaS  cel! 

Sample  vape>  wife  mass  spectrometer. 
Pressure  ft 

Best  accuracy — repwvs  Enter^-a;  sta nuard 

DovSe  or  Triple  cob  system 


... — Gas 
Solid 


• High  Temperature  Vacrngn  Furnace  for  tiguidus  lines. 
» DTA 


• Modified  cnm-Tterea:  magnetic  serdar  instrument 
Windows  based  Data  asguisrSas  system 


* Standard  metaSographic  preparation  + opbca:  electron  microscopes 


Measurement  of  Thermodynamic  Activities 
Measure  vapor  pressure  of  element  in  alioy- 
compare  to  pure  materiai-a(Ai)  = P^fAlloyJ/P^fPure  Al) 


Fe-AI  System  at  1573  K 

(Ion  Current  Ratio  Technique) 
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Al  Activities  in  oc2  Ti3Ai  + y TiAl 

(Double  Cell  Technique) 


<4>3£  Sf XXJK, 


Ti  Activities  in  a2  Ti3AI  + y TiAl 

(Double  Cell  Technique) 


1700  1600  1500  K 


TI45IC.8FW 


High  Temperature  Vacuum  Furnace — 
Liquidus  Data  for  Refractory  Containing  Systems 


Other  Methods 


OPTICAL 

PYROMETER 


E ( 

I 1 


HEAT 


W-SRa/W-SSH* 

THERMOCOUPLE 


TUNGSTEN  MESH 
HEATING  ELEMENT 
2*2”  HOT  ZONE  — " 


- TEMPERATUHE  TO  3080  “C  OIFFUSON 

PUMP 


Z.lYJSVl 1. 

Experimental  Thermodynamics  Research  Activities 
in  USA 


» Equilibrium  phase  composition  by  WDS  or  EDS 
analysis 

* Metallography-phase  amounts  by  lever  rule 

* XRD-phase  boundaries,  site  occupancy 

♦ DTA-phase  boundaries 

♦ Diiatometry-phase  boundaries,  CTE’s 

• Diffusion  couples-phase  boundaries 


if. : SS TOW W 

n iwwxxzxi-? 

Experimental  Thermodynamics  Research  Activities 

■ France 


* Ofe  Kleppa-Direci  Synthesis  Calorimetry,  XRD 

binary  slims 


Gachon  & Hertz-U.  Nancy,  Ccfne?-U.  J.  Fourier, 
Castanet  & Bros-U  Marseille.  Legersdre-U.  Paris 


Philip  Nash-Direct  Synthesis  Calorimetry.  DSC.  XRD 

Jsnary.  ternary  alloys,  higher-order  systems 

Ramana  Reddy-EMF 

binary  alloys 

A,  Navrotsky-Direct  Synthesis  Calonmetry 

oxides,  nifties 

N.  Jacobson.  E.  Copland-Knudsen  Effusion.  High  T 
liquidus 


Germany 

Sonwer-Max-Planck-Stuttgart . Schalier-U  Kiel,  Hitpert-fMES 

Poland 

FferterS.  Moser-MMS-Krakow 

Japan 

YamagucftHwate  U,,  Jacob-Tohoku  U,.  Shoji -Nagoya  U,. 
Yamsna-Kyote  U„  Sheng-ONRI 

Austria 

A Mikula  & H.  ipser-U,  Vienna 


binary  ternary  alloys 


Italy 

Borzsne-U  Rome,  Ferrs-U.  Genoa 


Workshop  Presentations 


71 


Summary 


« Thermodynamic  and  phase  equilibria  can  be  measured  by  a 
variety  of  techniques  and  with  varying  accuracies 

« AH  techniques  have  inherent  limitations 

* CaJonmeiry  is  the  most  versatile  technique  but  requires  different 
types  of  calorimeter  to  obtain  ail  of  the  data 

* Most  data  results  from  studies  of  individual  systems  as  part  of 
some  broadef  project  objectives 

• Systematic  investigations  are  needed  to  yield  consistent  data  that 
can  be  used  for  understanding  alloy  behavior  and  as  benchmarks 
for  first  principles  calculations 

♦ Very  few  research  groups  in  US  involved  in  experimental 
thermodynamic  or  phase  equilibria  studies  of  alloys 


Additional  points  made  by  speaker: 


Concern  was  expressed  that  the  number  of  centers  not  only  in  North  America  but  also  world  wide  with 
expertise  in  experimental  determination  of  thermodynamic  quantities  is  very  small  and  that  abolishment 
of  one  of  these  centers  jeopardizes  the  availability  of  quality  experimental  data  for  thermodynamic 
quantities. 

Improved  accuracy  for  measured  data  is  an  important  goal. 
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Configurational  Entropies  in  Real  Alloy  Phase  Diagram  Calculations 
Fan  Zhang,  CompuTherm  LLC,  Madison,  WI 


* is  Kcal  Ukn  fbaw  Wagreip  < akwhufryB* 


Configurational  Entropies  iti  Heal  Alloy 
Phase  Diagram  Calculations 

Fan Zhang 

ux:.  4S~a.  mi sm««  n-mm.  vs.a 


Ci 

Coll  ahoraiors: 

W,  A.  Oates 

Stimx  Hexeanh  batiste,  ! nmmfy  -lil't  « A 

V A.  Chang 

Ikgpanmtn  Satnetumd h>zimsmi;z  fmimstivof 

ifmofmaXkisSmi.  Afo<,%<>«,  jr/.yrnj.  i sj 
S.-S ..  Cb$t 

CstmOhemiM:  O's  rt&ntmm  \Um;m  m$3Tf».  f.vs 


gjf  < omptiTbi'rni  LL< 


< <}8fi2srmi<«M38  i ut  ^ 


2 tt&gnMtt  < staalatton 


Why  Phase  Diagram  Calculations  ? (1) 

• 8«w  alleys  with  desired  properties  are  constantly  needed 
in  today’s  world. 

* Phase  diagrams  are  road  maps  For  alloy  design  — they 
are  the  starting  point  to  understanding  any  phenomenon 
or  process. 

* Multicomponent  Phase  diagrams  are  needed  since 
commercial  alloys  usually  cuntain  more  than  lour 
com  ponents. 

• It  is  impossible  to  determine  such  information  purely  by 
experiments. 

^ CampuThemi  U.C 


Why  Phase  Diagram  Calculations  ? (2) 

L Large  smother  of  sot>-syste» 


IPS 


i otai:  31752  * 


15 


I ■ *’  Its:  ■«».' 


i Z'.  ymlawiL- 


2.  Multicomponent  phase  diagrams  can  fee  extremely  complex 

I hcrefore,  computationai  approaches  are  needed  to  reduce 
the  experimental  effort. 

< omjniTt«*!  ni  U.( 


Computational  Approaches 

• ’ First  Principles’ 

* Pltetiomettological  fCalphatli 

Successful  aspects  of  Caiphari  approach: 

* Many  self  consistent  thermody  naroic  descriptions  have 
been  published  for  binaries,  ternaries,  even 
quaternaries. 

• I hermody  tiumic  databases  for  moliicom ptmeot  systems 
base  been  developed. 


CempoThenn  U.C 


< I <rt«vpfer«  la  Real  Wins  Phase  < a&eufcrtfeww 


Objective  of  Calphad  Approach 

I dea:  • Describe  the  Gibbs  energy  of  each  phase  in  a 
sy  stem  using  a simple  empirical  model  equation. 

* Optimize  the  model  parameters  using  She  available 
experimental  inihnnation. 

* I xtrapoiaie  the  use  of  the  model  parameters  to 
temperatures  and  compositions  where  the 
experimental  data  are  not  available. 

Objective:  to  “predict"  unknown  properties 
(multicomponent  alloys)  by  known  information 
(binaries  & ternaries). 


< onfuorta i ntroftfc*  an  Real  \fo>  Ffeaw  Ducom  < akntrifem- 

Criteria  for  the  Successful  Models  and 
M odel  Fa ra  meters 

♦ How  good  they  are  in  describing  the  known 
information  (phase  diagrams,  thermodynamic 
properties,  and  so  on)? 

♦ How  good  they  are  in  predicting  the 
unknown  properties  (extrapolation)? 


$ CompuTherm  LLC 


Qjf  e«inj!iiTb(’n»  \\X 
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Ctoifisttralmnat  turnon*'*  m Real  AS&y  }*fea*e  Dussraa*  < ahs&Suto* 


i ootisorattonat  l m Hc*l  \Stov  !*#«***  UtagnMo  < 


The  Calculated  ISi-Al 
Phase  Diagrams 


{.  aiphad  Assessment 


f;mt!|HlThtTM)  UX  » > 


< totmwtio;  in  ilmlMUsy  I'havc  l.*ea£T»ttt  i.  ::*toufc»tto8<e 


A Realistic  Approach  for  Phase  Diagram 
Calculation  at  the  Present  Time 

i 1 1.  Modeling  equation  should  include  a dash  more  «f  physics 
so  that  model  parameters  arc  more  reliable  in  the  'prediction' 

{2).  Modeling  equation  should  he  simple  so  that  it  is  readily 
applied  to  multicomponent  systems. 

Improve  ihe  Description  of  Configurational  Entropy 
Currently  l sed  hy  ( alphad  Approach  (BfV) 
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Cuat&aimtottol  £mref>k>  to  Real  \a»<  ?•&«**  Diagram  < akalattou* 
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i*  CompuTherm  IJLC 


* ftwiroiHs.*  m ttrrf  \Jtoy  b*tia*r  ttosgram  < atouSatesa* 


i 8»t«s*erattt>aa8  Lnt  rv^R-s.  ia  Real  AS«a  Pfcs'*-  Diagram  < a too  teutons 

The  Modified  CSA  |99f)al] 

Original  CSA: 

Modified  CSA: 

My  using  t]  as  an  adjustable  parameter,  more  fteeibility  is 
obtained  in  describing  the  Gibbs  energy -composition  curves. 

By  introducing  Ff , the  cluster  interference  is  approximately 
taken  into  account  by  a two-term  repression  in  the  modified 
CSA.  fesshi-wrm  evpre sskm  is  used  is  the  € VM ) 


^ Compufhomi  U.C  s’1 
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i fa  Kc-d  Vffos  Pbaw  jfrttUffrsfl* 

Fee  Coherent  Diagram  in  MC,  O M.  CX  SA  amt  MC'S  \ 
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< owpitThi-rnt  U.,( 


< «at«»ntt*0a&t  m Real  Phaxe  l>u$r»m<  atfo4»?K-^« 


Application  of  the  CSA  to  Binary  Systems 

* Ati-Cu  system 

• Ati-INi  system 

* Cd-Mg  system 

• Al-Ni  system 


< XstnptiTlirmi  IJX 
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Cd-Mg  Phase  Diagram  by  FP  and  CSA 
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Apply  the  CSA  Model  to  the  M-Ai  System 
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CompuTherm  LIT 


& CompuTherm  LLC 
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ff>  R«at  Mlaj  Pits'*  Diagram  < a’orttialion* 


< Xafeuraliwtwl  Fairvpics  in  Real  Aftay  Phase  Dost  am  i ak-«Ut»»»» 


Summaries 

• We  art-  NOT  looking  for  iht  most  accurate  ratettianon 
ofeonfigstmtiioital  entropies  (MC  or  large  cluster  CVM) 
as  in  PROTOTYPE  phase  diagram  calculations. 

• We  arc  looking  for  a calculation  method  suitable  for 
RF.AL  alios  phase  diagram  calculations  which  giscs 
correct  topologies  (Bragg- W illiams  doesn’t  for  no  SRO). 


The  CSA  mode!  is  very  promising  for  this 
purpose. 


t ompitTh '(  to  U.C 


C-omptiThmn  S.  W.: 
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Additional  points  made  by  speaker: 

Improved  models  are  needed  to  avoid  physically  unreasonable  Gibbs  energy  functions  in  regimes  where 

the  phase  is  not  stable. 

Example  metastable  y/y  in  Ni-Al: 

• Description  of  both  phases  as  two  independent  phases  with  no  order/disorder  transition  resulted  in  an 
unreasonable  diagram  in  regimes  of  the  phase  diagram  where  these  phases  are  metastable. 

• Treatment  of  y and  y as  one  phases  and  the  use  of  Bragg-Williams  approximations  for  the 
calculation  of  the  entropy  of  mixing  do  not  generate  the  correct  fee  metastable  phase  diagram  due  to 
the  neglect  of  short  range  order. 

♦ Diagrams  obtained  from  first  principles  and  CVM  (cluster  variation  method)  or  MC  (Monte 
Carlo)  give  the  physically  most  probable  metastable  diagram  (work  from  Pasturel  and  Colinet). 

Example  metastable  fee  diagram  in  Ti-Al: 

• Results  from  different  model  treatments  show  diagrams  with  similar  features  to  Ni-Al. 

Questions/comments  on  various  models: 

• Bragg-Williams  as  a point  approximation  method  cannot  describe  all  order/disorder  transformations 
correctly. 

• The  expression  for  the  entropy  in  CVM  is  a too  complex  for  practical  purposes.  Although  these 
phase  diagrams  are  topologically  correct,  the  accuracy  of  results  is  insufficient  for  practical 
applications. 

• Cluster/site  approximation  (CSA)  in  its  pure  form  has  the  problem  that  the  disordered  phase  is  stable 
at  0 K.  This  problem  is  overcome  by  the  modified  CSA  with  an  adjustable  parameter.  This  model 
gives  results  that  are  in  good  agreement  with  results  obtained  from  MC. 

Questions  & Answers: 

What  about  the  vibrational  entropy?  - The  model  includes  additional  parameters  for  the  entropy  to 

account  for  this. 
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Objectives 


• Understand  Trends  in  Ordering  Phenomena. 

- Ordering  trends  in  bcc-based  alloys  (with  V.  Drchal  and  J.  Kudrnovsky) 

• Predict  Phenomena: 

- Transient  Ordering  (with  L.  Reinhard) 

e.g„  Ti-V,  Fe-Cr, ... 

- Ordering  in  Complex  Phases  (with  A.  Finel) 

e.g.,  A1S, ... 

- Pressure-induced  Ordering  (with  P.  Singh  and  G.  M.  Stocks) 

e.g.,  Al-Ge,  Al-Si,  ... 

~ Provide  Thermodynamic  Data  for  Ab  Initio  Phase  Diagram  Construction 

~ Provide  Estimates  of  Heats  of  Formation  (from  total  energy  vs.  volume  results) 

~ Provide  Input  for  Phenomenological  Approaches  such  as  CALPHAD 

~ Challenge  the  validity  of  some  assessed  phase  diagrams,  and  the  validity  of 
ab  initio  results 


Background  and  Motivation 


Phase  Diagrams  are: 

• The  basic  road  maps  for  alloy  designers 

• The  graphical  results  of  battles  for  survival  between  competing  phases. 


Orders  of  Magnitude: 

• Total  Energy  1<P  Ry/atom 

• Formation  Energy  10-3-10-2  Ry/atom 

• Ordering  Energy  104  - 10-2  Ry/atom 

Units: 

1mRy~  13.6  meV 
1mRy~  157  K 

ImRy-sy  1312.76  J/mol  (or  J/g.at) 


Ab  initio  Electronic  Structure  Methods 


Constraints  within  Dersssty  Functional  Theory  ©FT) 

• Local  Density  Approximation  (LDA) 

Exchange  Potential,  e.g.: 

- von  Barth  & Hedin 

- Vosko 

- Ceperley  & Alder  (Perdew  & Zunger) 

' Beyond  LDA 
LDA+U 
GGA 
SIC 
GW 


Other  Constraints 

• Shape  of  the  Potential  Function: 

- Muffin-Tin 

- ASA  (equal  or  non  equal  sphere  radii) 

- Full  Potential 

• Core  versus  Valence  Electrons 

• Scalar  Relativistic  versus  Fully  Relativistic 

• Magnetism:  Co-linear  (Ising-like),  Non  Co-linear  (Heisenberg-like) 

• Atomic  Positions  and  Site  Occupancy 

• Bom-Oppenheimer  approximation  (electron-phonon  coupling) 


Methods  Limited  to  Ordered  Structures: 

(Wave  Methods) 

- Pseudo-Potential 

- LAPW  or  FPLAPW 

- LMTO  or  FP-LMTO 

- ASW 
-FP-LASTO 

Methods  that  cart  also  handle  Disordered  Structures  (via,  e.g.,  the  CPA): 

(Green's  function  Methods) 

- KKR 

- TB-LMTO 

Properties  : 

• Equilibrium  Properties 

- lattice  parameter(s) 

- Structural  energy  differences 

• Elastic  Properties 

- bulk  modulus 

- elastic  constants 

- phonon  spectrum 

• Electronic  Properties 

- densities  of  states  ( y of  heat  capacity  data) 

- spectroscopic  data  (ARPES,  UPS,  XPS,  PAS) 

• Transport 

- conductivity  (DC  and  AC) 

- reflectivity 

• Magnetic  Properties 

- magnetic  moments  and  magnetic  ordering 

and  Alloying  Effects  on  Properties 


Software  Packages 


is 


Pseudopotential-based  Methods: 

VASP 

Siesta 

ABINIT 

Full-potential  Methods: 

Wien97/Wien2000  (FPLAPW) 
FPLMTO 


78 


Workshop  Presentations 


Order  versus  Disorder  for  a Binary  Alloy  Ai-c  Bc 


From  tbs  Disordersd  State 


eb 

£»b/ 

eI 

£o« 

J 

From  the  Ordered  State 

3 


(1-=C)eI  - cCEb  = *u«(ec)»  £.i,s 
Mixing  En«rgy 

^uix(c)=  ETcp*  (c)-(i-c)El  -cEb 
Ordering  Energy  of  Configuration  C 

-&:0:d  - ^2* ipk-s  $pU-c j 

/Eord  Xq,  V,  where  qS  = cc(n“- ccn,)/2 
V.=  +V“  -2Vj* 

V,  >0  — w (A-B),  pairs  favored  — ► ORDERING 
V,  <0  — w (A-A)«  or(B-B).  pairs  favored  CLUSTERING 
‘Ordering*  Energy  associated  with  Segregation  if  n»  = n, 

• Ef 


Formation  Energy  of  Configuration  C 

j£fo:w-  Vv-  <£...»«*  V«  <fc  <.2pr.-1.-t 

[V]  Volume  =[  .[4EForm]Volume 


^ tnGnmtmGmn  c 


CALPHAD  Thermodynamics 


Mutti-coroponent  Phases  . 

G* 1 2 * 4({6  },T)  = f*fG  *({c  i },T)  +ld  G *({c , },T)  +“G  4({c  i},T) 

™'G4({Ci},T)=  Zc,  Q+({Ci},T) 

(1)  ld  G 4({ci  },T)  = RT  Xci  In  c, 

“ G*«Ci  },T)  = See,  X’  L*(T)(ei  - c if  p Liter)  = u + vT  + ... 

|,j>J  p=0 

Two-subtattice  Model  of  a Binary  AHov  . 


(2)  ld  G *({yf  },T)  = RT  ZNs  Z yf  In  yf 

s=J,ll  i 

" G*({yis  ).T)  = yi  yk  [ Ey!  S’ La.BiI (T)(yi  - yB)P  ] 

i p=0 


Ordering 


+ yi  y*B  [ Eyi  (T)(yS  -y6f  ] PLi,j:k  =PLitJ:k* 

i p=0 

G c ({c,  },T)  = G 0l5«c , },T)  + AG°,d  «yf  ),T) 


Eqs.(1)‘ 


ACt3'**  (T)  = xsG°rd((yf  ),T)-**G°rd((C|}.T)  Eqs.  (2) 


The  Role  of  Ab  Initio  in  CALPHAD 


Heats  of  Formation 
Heats  of  Transformation 
Thermodynamic  Functions 
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Assessment  of  Ab  Initio 
Phase  Diagrams 


Thermodynamic] 

Kinetic  ] 

Database 

Database] 

4 N 4 

=t>teDl«>i 

CALPHAD] 

Kinetics  j 

From  the  Disordered  State 

<£c>,i  -1.X*  { PivC*'5(pEh.«J1'>Vw(e''; 

Requirement 

Green’s  function-based  method 


From  the  Ordered  State 

Requirement 

Wave  (or  Green’s  function)  method 

Advantages 


Easy  to  implement 
nited  tc 


Advantages 

* Decoupling  between  the  energy  of  the 

disordered  state  and  the  ordering  I * Not  limited  to  simple  structures 
energy  j • In  some  instances,  no  constraint  on 

* Fast  convergence  of  the  expansion  in  j the  shape  of  the  potential  (FP-LAPW, 

terms  of  effective  interactions  FP-LMTO) 

* Concentration-dependent  interactions 

« Actual  prediction  of  the  ground-state  | Limitations 

properties  « Selection  of  the  number  and  type  of 

| clusters 

Limitations  j * Selection  of  the  ordered  configu- 

* Implementation  so  far  within  ASA  and  rations 

Muffin-tin  j * The  interactions  are  concentration- 

* Questions  related  to  charge  transfer  j independent  except  via  volume  effect 

effects  i ♦ Slow  convergence  of  the  expansion  i 

In  an  ab  initio  context,  only  simple  terms  of  effective  interactions 
structures  (fee,  bcc,  hep)  are  dealt  with  9 Very  questionable  predictability 

Unified  Approach  to  Alloy  Stability? 

Can  thephysicsof  the  chemically  disordered  state  and  of  any  ordered 
configuration  of  an  alloy  be  described  within  a unique  framework? 


Ab  initio  to  CALPHAD 


Ab  Initio(M) 

AGAUc.T)  = AEaMc)  + AEord(c,T)  - Tdsfflc.T)  where  AEMitc)  = c(1-c)[°L  + (2c-1>  L] 

CALPHAD  (CP) 

AGPfl[c,T)  = G(c,T)  - ZcjG? 

i 

= +RT  Scjln  cj  +.  Zcjcj  XPLjj  (cj-cj  )**  where  PLjj  = u + vT  + ... 
i ijx 

Determination  of  the  {u,v}  parameters  by  minimizing  the  difference 

£GW  - ASP^c/T) 

The  assessment  can  be  performed  with  the  following  input  fresh  initia 
• Gibbs  energies  and  enthalpies  of  formation  as  functions  of  temperature  and  alloy 
composition. 

• Phase  diagram  information 


Ab  Initio  Methodology 

ir 

• Tight-Binding  Linear  Muffin-Tin  Orbital  (TB-LMTO)  method  for  the  description 
of  the  electronic  structure 

• Atomic  Sphere  Approximation  (ASA) 

• Unequal  sphere  radii,  zero-charge  transfer  (i.e.,  local  charge  neutrality) 

• von  Barth  Hedin  exchange  potential 

• fully  relativistic  treatment  of  the  (core  and  valence)  electrons 

• Coherent  Potential  Approximation  (CPA)  for  the  treatment  of  the  chemically 
random  state  of  the  alloy 

• Generalized  Perturbation  Method  (GPM)  to  determine  the  Effective  Pair 
Interactions  (EPIs) 

• Cluster  Variation  Method  and  the  Newton-Raphson  technique  for  the 
determination  of  the  equilibrium  free  energy 
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Input  from  Ab  Initio:  Heats  of  Formation 


[=■ 


Ab  Initio  Results:  Formation  Energy  of  the  Chemically 
Random  Solid  Solution 

• *“>■'-  »+ 

46tS*(c')=£  (1-6  It  L*(4  -t)  L] 


Interface  between  Ab  Initio  and  CALPHAD:  The  Ni-Cr  case 

Ni  2Cr  Heat  of  Formation  obtained  from  ab  initio  was  used  to  include  the 
ordered  phase  Cl  1b  (or  OP6)  in  the  phase  diagram  of  Ni-Cr. 

OP6  (Cl  1 b)  Ni  & (LMTO-ASAr«ullt) 

’ I 1 ’ 1 1 I 1 1 »'■  1 '■’T~^rrT 


Ab  Initio  Results:  Ground-state  Analysis 

iSoed-  V* 


\ Vs-W  j. 

. 

fa-W 

■ \ / 

• 

\ / 

i 

i 

Comparison  between  Ab  lnitio+  CVM  and  Ab  Initio 
•f  CALPHAD  Results  for  Ta-W  and  Ta-Mo  Alloys 


m 


Comparison  between  Ab  Initio  + CVM  and  Ab 
Initio  + CALPHAD  Results  for  Ta-W  and  Ta-Mo 

Molar  Gibbs  Energy  versus 
alloy  composition  at  various 
temperatures 


02  0.4  0.6  0.8 

Mole  Fraction  Ta 


0.2  0.4  0.6  0.8 

M ole  Fraction  Ta 
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Ab  initio  + CALPHAD  Predictions 
of  Phase  Diagrams  for  Mo-Ta-W  Alloys 


Future  Directions  for  Ab  Initio 


• Electronic  Structure  Methods 
Electron-electron  correlations 
“Charge  transfer"  effects 
Relativistic  effects 


Ab  initipv+  CVM 


Conclusions 


Ab  initio  approaches  provide  information  on: 

• ground-state  properties  and  ordering  trends 

• the  origin  of  stability  and  local  order  in  alloys 

• thermodynamic  functions  and  phase  diagrams 

• short-range  order,  validated  by  experiments 

The  ab  m/f/oCALPHAD  approach  can  be  used  to: 

• supplement  thermodynamic  databases 

• reproduce  the  main  features  of  predicted  ab  initio  phase  diagrams 

• predict  the  impact  of  ab  initio  predictions  for  binary  alloys  on  multi- 
component  phase  diagrams 

• validate  theab  initiodata  by  using  the  CALPHAD  predictions  on 
multi-component  phase  diagrams 


Examples  of  ‘'still”  challenging  systems: 


> Liquid  state 


• Hydrides  (e.g.,  Pd-H) 


• Transition  metal  carbides,  nitrides,  and  oxides 


• Unified  Approach  to  Alloy  Stability:  SRO  vs.  LRO 

Can  theqhvsicsof  the  chemically  disordered  state  and  of  any  ordered 
configuration  of  an  alloy  be  described  within  a unique  framework? 

• More  accurate  prediction  of  the  free  energy 

Relaxations 

Vibrations 

Magnetic  contribution 


• Magnetic  alloys  (e.g.,  Ni-Cr) 

• Shape  memory  alloys  (e.g.,  Ni-Ti) 

• Correlated  systems  (e.g.,  rare  earth  and  actinide-based  alloys) 


Issue  pertaining  to  lattice  stability 


•Treatment  of  Multi-component  Alloys:  Connection  to  CALPHAD  approach,  and 
user-friendly  interface  to  software  applications  such  asThermo-Calc 

• Description  of  the  Liquid  and  Amorphous  states 

• Kinetics  of  Transformation 

Under  isothermal  conditions  (TTT  diagrams) 

Under  non-isothermal  conditions 

• Dynamics  of  Transformation 

Fundamental  understanding  of  complex  phase  formation  and  transformation 


Even  in  the  simplest  description  of  the  vibrational  properties  within  the 
Debye  Model,  the  Free  energy  is  given  as: 


F(r,T)  = E*i(r)  - kBT.[D(eo/T)  + 3ln(1-e  -0d/T)] 


• zero-point  energy 


Instable  Phase: 

0d  is  not  defined,  hence 
the  total  energy  is  not  defined 


Additional  points  made  by  speaker: 

Background:  Phase  diagrams  are  graphical  results  of  the  survival  battles  of  competing  phases. 

Ab  initio , CALPHAD  and  experimental  phase  diagrams: 

• Comparison  of  ab  initio  and  experimental:  Ni-V,  Pd-V,  Cu-Zn 

• Coupling  of  ab  initio  and  CALPHAD:  Ni-Cr 

• Comparison  of  ab  initio  and  CALPHAD:  Mo-Ta,  Ta-W  (very  good),  Mo-Ta-W 

Ab  initio  can  provide  good  information  for 

• Ground  state  properties 

• Stability  and  ordering  in  alloys 

• Short  range  order 

• Supplement  thermodynamic  databases 
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Limitations  and  challenges  for  ab  initio : 

• Liquid  phase 

• X - H systems 

• Transitionelementcarbides,  -nitrides,  etc. 

• Magnetic  properties 


See  also  “Some  Notes  on  First-Principles  Electronic  Structure  Methods  and  Calculations,”  page  17. 
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Database  Development  and  Industrial  Applications 

Arthur  D.  Pelton,  Centre  de  Recherche  en  Calcul  Thermochimique  (CRTC),  Montreal,  Quebec 


Thermodynamic  Database 
Development  and  industrial 
Applications 


Arthur  Psitoti 

CBCt 

(Cesstre  d«  RecbM.Vsa  »rs  Catoyl  T her rmsc h i m > y y e ; 
Ecote  PotyseobBiqu®  de  Mantes!, 
tterstriat,  Caw* * 


FactSage 


riMoeparar 
CM  S3}<?<  AD,! 
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Copyrght  Thermfael  1S7S  2002.  AS  rights  r esen/wt 

»;>**<  -'m  i i wp 


FactSage  5. 1 


Tn;>;  & MmWV 
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Steps  In  Developing  Multicomponent 
Solution  Databases 


* Search  Ite&tvre  for  all  available  thermodynamic  and  phase 
fiquiShnnm  date.. 

« Select  appropriate  model  for  each  phase  (Gibbs  energy  as  a 
function  of  T and  composition). 

8 Critically  evaluate  and  optimize  each  hutary  subsystem  to  obtain 
model  parameters. 

• Use  model  to  estimate  properties  of  ternary  phases.  Add  ternary 
model  parameters  where  necessary. 

« Use  model  with  binary  and  ternary  parameters  to  predict  Gibbs 
energy  of  myfUsempeRehf  solutions, 

8 Calculate  multicomponent  equilibria  try  Gibbs  energy  minimization 
using  databases. 


fcvafuafed/Optimized  F*A*C*T  Solution 
Databases 


Choice  of  models 


Single  subtattice.  random  mixing: 

- Simple  alleys 

~ Simple  comroom-ion  ionic  systems 
Sublattice  models,  random  mixing 
{Compound  Energy  formalism}: 

- interstitial  solutions 

- Complex  solid  alloys 
~ Long  range  -ordering 

- Ceramic  solutions 

~ Simple  molten  ionic  solutions 
Sublattice  models,  short-range-ordering: 

Molten  slags  and  glasses 
■ MoRen  sulfides 

- Reciprocal  salt  solutions 


The  predictions  ere  only  as  good  as  the  models  The  proper  choice  of  a model 
which  correctly  reflects  the  structure,  ordering  and  speciation  of 
a given  solution  Is  essential. 


Evaluated/Optimized  PA*C*T  Solution 
Databases  f con  lift  tied  f | 


® Alte/s-Carbonitrides: 

SGTE  Solution  Databases  (ferrous  and  non-famsus} 

* Molten  Oxide  Slag/Glass: 

Oxides  at  Al.  As, 

8,  Ca.  Cd*.  Cte.  Cu,  Fe**.  Fe»*,  K.  Mg.  Mrs,  Na.  Ni.  Pb, 

Si.  W- 

Tts>.  It,  Zr>,  Zr. 

~ Dilute:  S.  SO^ 

PG^CO^.F  CUHjOiOH 

8 Scad  Oxides  (ceramic}; 

- ilsnenite 

(Fe.  Mg,  Mrt.  Ti>  TiOs 

- Psandebrooteie 

(Fe  Mg  Mrs  Ti)  Ti?Oj, 

- Monoxide 

(Ca,  Cu.  Fe2*,  Fe*”.  Mg,  Mn.  Ni.  Ti.  Zn.  Zr.  Al)0 

~ Garnet 

(C-3  F^Vj/Cf  A!. 

- Spinets 

(Fe  Zn.  As,  Mg  Co  Nr.  Cr  Mn)j04 

- Forsteote 

(Mg.  2n)2St04 

- QtivfCte 

(Ca  Mg.  Fe)zSsOa 

- Corundum 

<Ai  Cr,  Fe),Os 

8 Solid  Oxides  (ceramic,  continued): 

- Melite  (Ca,  P&yZn,  fe*».  Fe3*,  AStfFi^,  At.  Si)jO? 

- Pyroxenes  (Ca.  Mg.  Fe)(Mg.  Fe)SijO«  (ortho,  cirno.  proto) 

- Spinels  (Fe,  Mg,  MnsTIjO* 

- Perovakite* 

- ZfOj-based  solutions 

- Feldspars 

- WoiiasicRite 

- Pseudowoilastsnite 

- Etc. 

* Sate  (liquid  and  solid) 

- U.  Na,  K Rb,  Cs,  Mg. Ca.  Si.  Ba.  At  Fe  [2*  3f.  la.  Ce  fid/?.  Cl.  Sr  I.  OH 
MQ.j.  CO;i.  SOs 

- Cryolite-CaF.rAl20,i-AM'ta 

- FeC^-FcCi/CoCS.,-NlCl<-MnCI,'-h1gCb-AiCir.. . 
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Evaluated/Optimized  PA*C*T  Solution 
Databases  (continued! 


« Sulfide  Mattes: 

C«-f  e-Ni-Co-Pb-2fi-S-A8-. 

» Moten  Aioys: 

- Fe  + 32  salutes 

- Pb  > 12  solutes 

- Sn  * 18  solutes 
~ Zn  * 4 solutes 

- AS,  Mg  * J solutes 

* Solutions: 

Concentrated  (nors-idsab 

« Semiconductors  {solid  and  liquid): 

Hg~0£i~Te"Sb~Se~Zn-')N<3''Pb"T:“G3~Q  i 


F*A *C*T  Consortium 


B 


N 

tmmis 

A 


Mm mm 


SCHOTT 


teck 


IIS  Materials 

Umvfersiiy  Tokyo 


mm 

1 &££> 


ft 


'tsSBltl  iwro 


CORNING^ 


fftfWn 

mm-'mmm 


Mg0"Fe0-Fe203"AI20s  system 


■ Accuracy  : good 

« Optimized  region  : Safe-Solidus  region 


Haiti  Solutions 


* Slag: 

C«gO<APOj#eO-  Fe,©.) 

i * Monoxide: 

PgO,  A§*0*  HO,  Fe,Oj| 

( Magi  msiowti stile  (MWii 

! * Corundum: 

{A ?A>  Fe.O-i 

a Spinel; 

im,  ?«**,  «*♦,  AFP  pg,  F<P\  fe3*,  Al3\  ¥a|,  04 

(Magneiite(Fe  }OJ.  Magnesiofemte(MF.  MgFe:OJ. 

hietvym>e(Ff)AF,Oj  Spm&i(MyAi,Oj) 

Oxygen  Partial  Pressure  over  the  MgAij 
FeAl:€L  Spinel  Saturated  with  A!,03  and 


Fe 


m 


e.a 

m&e  fw« im 


IfcV^O* 


Enthalpy  of  Mixing  of  the  Metasfahie 
Wi$&ltCA~¥eM20A  Spinel 


{• 


X 


69  o i a?  a 3 


r$ 

MsilA 


The  Miscibility  Gap  in  the  Fe304-SSgFe204 
m§ALQ4*F£M2Q4  Spinel' at  BOCFG  * 


m 


MgFe?G4 


Fe304 


FeAUO, 
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j Ca0^g0“Fe0“Fe203"Si02  System 

m 

* Accuracy,  good 

* Optimized  region:  Whole  system  exsept  Ca,St04  region  j 

Main  Solutions 


- Siag : 

« SfefsoxM© ; 
« Spinel : 


« OSieine : 


|C3O-Mg0-F«0-F«.0rSIO;4 
iCaO.MgO.fto.r^jO,!  (Om.  wusm) 
|f«*%  f®>\  Mg)  W&\  F#,  Mg.  Vsj?  0« 
(Msgrwtitei^OJ.  Msgimiotetm(MF. 
(Ca,  Mg,  Fe>}  [Ca,  Mg,  Fsm  Si.04 
{ForstHftlH(hkj2SiOJ.  FayDlitefFe^iOj, 


MoniiceiSitetCaMgSiOJ,  yCajSrCy 

- Pj»at«  CC*  Mg,  Fr*f  * ! f»  if  Q6 

• proto  pyroxene,  orftsts  pyrayafsa.  cSsiSi  pyroxene.  io»  diao  pyroxene 
» WoHastonite : fCa.MglSiol 

- cs'-CajSsO* : {CfcMgijSiG* 

| » SrCa^SiO,  ; fC&,Mg,Fe*')gSsO<) 

| » MeUNte : Ca?  [Mg,  Fe^jSijO ?(Akerm3fiite,  Fs-Aksrmmiie) 


f » Fe-Ugutd  : PA'CTT 


The 

in  the  Ca?Si04-Fe,S: 


Gap 

gj5i0„  Olivine 


Ca.SiO. 


MgJSiOt 


mote  finctio ft 


The  Ca$i0,-Fe5f0.r&lgSiQ 
Phase  Diagram  at  1 2SG°C 
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System  CaMgSi,Ofi-CaFeSi20§  Section 


Saturated  with  metric  Fc 


l aivSUh 


System  Caf^gSijOg-Mg-SiO^-FeO  in  Air 


CaMgSips 


/ 

| gpjftgj 

\ \/_^\ 

< ■ ■ v ....•>  . . . ^ 

$$ 0 ^5  $ 04  vtisioffit 

^ **  FeO 

' *Jif»  yycicxerc 


0 


r ^ 


S«SF»fe>«->Spir»i  ~’-Xs 


</$■*%»  CfeJSS*  «&$  ${<xa&t?  f«si  tfcS&CSssS  «&t$  Ssr  r8<y  &»*&  *&»&£  ${  $&&$&£&*» 


System  CaMnSl,0,,M§.MQ4 
at  PfO2)=10'6  atm 


CaMgSiPs 


Mgpio,  - - ~ •* 


FeO 


Ltgwidus  of  CafVtg* 


*t206*F®0  Section 

IS4  ate 


* 
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The  Application  of  FactSage* 
Software  to  SteeSmaking  Processes 


MfhctSage" 


Feb..  2002 

In-Ho  Jung,  Youn-Bae  Kang  and  Arthur  D.  Pelton 


| • inclusion  modification: 
j Al203  modification  by  Ca  treatment 
\ « minimizing  the  losses  of  Ti  In  low  carbon  steel 
i * inclusions  in  stainless  steel: 
j p|,  fit]  and  f€j  


CRCT,  Ecole  Polytechnique  de  Montreal,  Canada 


1,  Inclusion  modification: 


AI203  modification  by  Ca  treatment 

\ * AijOj  inclusion  is  formed  after  Al  deoxidation  process  f 
i and  very  harmful  to  steel  quality  and  nozzle  clogging 
| * So  Ca  treatment  is  performed  in  the  secondary 
| sfeelmaxing  process  to  modified  Ai,Q,  to  low  melting 
j temperature  composition 

I » To  understand  Ca  treatment,  the  inclusion  diagram  of  Ca  j 
j and  Al  is  necessary 

! « The  effect  of  S to  the  Ca  treatment  is  calculated 
| ® Real  Ca  treatment  process  is  simulated 


Simulation  of  Ca  treatment 
fCaSi  wire  infect  ionf 


j 2,  Minimizing  the  tosses  of  Ti 

MM\ 

| In  low  carbon  steal 

JESLi 

AA&taxt  sst  £*&  *>  (9*  i*a*  of  itiei  j 


\*  Ti  is  aided  In  the  low  carbon  steel  os  alloying  element 
I in  the  RH  process  after  Al  deoxidation 
| » Bui  because  Ti  is  strong  oxidizing  element,  Ti  can  he 
j consumed  to  formation  of  TiO*  inclusion 
j « To  minimize  the  losses  of  expensive  TI,  the  inclusion 
j diagram  of  Ti  and  A I is  necessary  
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i 3.  Inclusions  in  stainless  steel;  [M|,  pi]  and  [C] 


f « Stainless  steel  contains  a lot  of  nitrogen 

I*  Wdh  addition  of  Ti  alloying  element,  Ti  east  fee  consumed 
as  TIM  inclusion 

* To  minimize  the  losses  of  expensive  II  and  reduce  the 
inclusion,  the  inclusion  diagram  of  Ti  and  N is  necessary 
» The  effect  of  carbon  to  the  TtN{TsN-Ti€}  formation  is  also 
I calculated 


Inclusion  Diagram  in 

S3W 

Fe-18  wf%Cr-8wf%  Ni-TI-N-O-C  at  188D°C 

***** 

| Flowsheet  for  copper  extraction  from 

{HIM 

sulfide  concentrates 

Databases  a nd  Software  for 
Thermodynamic  Simulation  of 
Copper  Smelting  and  Converting 


Sergei  A.  Degterov  and  Arthur  D.  Pelton  I 

Ecole  Polytechnique  de  Montreal 
Centre  tor  Research 
in  Computational  Thermochemistry 

Manuel  ZaraaSloa 

Nor8u<ta  Technology  Centre 


~1250VC 
P(0.)  = tOMO”  atm 
P(Sy)  ~ 10s  atm 


"N 

J 


r 


Copper 

Sulfide 

Concentrate 

h 

t\ 

I \ 

r \ 

25-30%  Cm  ; 
28-34%  Fe  j 
28-34%  S | 
2-6%  SiOj  } 
Smell 

amount  o’:  I 
Pb.Zn,  Nl.  j 

Ca.  etc  . j 


SMELTING 

Blast  Furnace 
Reverberatory  Flash 


-1250  C 
P(Oj}  = WMO’-1  atm 
P(S);~  10-*-10-?atm 


Mitsubishi 

1 K 

Noranda 

j:  \ \ 

(40-73%  Cu  1 

30-6%  Fe  | 

| 

I 30-20%  S | 

t 

Slag 

mai.ii  ■*  CONVERTING  r* 


Blister 

Capper 


IV 


0.2-2% 

S.2-1.3% 


Cu 


sv - 

\ 2-8%  Cu 
: 25%  SsG, 


Slaa  Discard 
or  Treatment 


Slag 

Cleaning 
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Sulfur  content  of  silica^aturatod  slag  In 
equilibrium  with  matte  and  liquid  Cu  or  v-Fe 


1:200%'  ( 22) 
iiWC(22i 
l2$<f:C<22) 

1300*0  OW 
120</C<31 » 

tC(2<ri 
ttWCCl) 
i 2 $i¥Cy  calcma^i 


0 ••>• 

0 


40  50  50 

Wojrh?  % C»  m Matte 


Copper  content  of  slag  in  equilibrium  wilt 
matte  at  125CFC  and  P{G,)-10S  atm 


90 


Workshop  Presentations 


Sulfur  content  of  liquid  Co  and  partial 
pressure  of  S->  for  the  Cu-matte  equilibrium 


1 Tgrstpersnsrsl  Kola  fraction  si  S 
| *G  1 in  fette 

I K»1a  fraction  of  S ! 
in  Csj 

, ates 

| &pjKS 

I Calculated  ] 

I Exp.***  | Calculated  | 

[ £xp.0*i  j 

| Calculated  § 

| 1150  | $.3280 

0.3280 

0.S2S1  { 0.0235 

1 4.0;  104  j 

j s.m  • kh  j 

1 1200  | 0.3240 

| 0 3245  | 

1 0-0288  1 00222  j 

3.75  ■ 10*  l 

0 87  < 104  | 

j 1250  1 0.323S  | 

| 0.3230  j 

| 0.0340  | 0.0337  j 

1 2 16  < 104  5 

1 2.C5  > 105  \ 

Additions  of  CaO  on  solubility  of  Cu  in  slag  in 
equilibrium  with  Cu-Au  alloy  and  St02 


tette/Slag/SiOj/Cu  equilibrium:  distribution  of 
slag,  matte  and  liquid  copper  at  1300°C. 


Additional  points  made  by  speaker: 

Choice  of  models  - although  model  choices  may  be  identical  in  subsystems,  the  use  of  different 
extrapolation  methods  is  likely  to  result  in  different  extrapolation  behavior. 

Cluster  pair  approximation  and  quasichemical  models  are  used  in  the  slag/glass  and  sulfide  mattes 
databases. 

F*A*C*T  consortium  includes  15  sponsors  for  database  development. 

FACT-app  (Chem-app)  software  interface  for  inclusion  into  other  software  programs. 
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Development  of  Thermodynamic  Databases 
Philip  J.  Spencer,  The  Spencer  Group,  Ithaca,  NY 


Development  of 
I Thermodynamic  Databases 


Background 


Thermodynamic  calculations  and 
simulations  based  on  critically  evaluated 
data  are  widely  used  as  a basic  tool  in  the 
development  and  optimization  of  materials 
and  processes  of  many  different  types. 

Such  calculations  can  result  in  considerable 
savings  in  costs  and  time  associated  with 
experimental  development  work. 

For  these  reasons  thermodynamic  database 
development  is  an  important  and  increasing 
scientific  activity. 

; Tic  < -rtwjjt 


Outline  of  presentation 


Available  databases  and  their  content 
Importance  of  awareness  of  experimental 
problems  in  carrying  out  data  assessment 
Importance  of  model  choice/compatibilty  in 
progressing  from  binary  to  higher-order  systems 
Funding  problems  associated  with  database 
development 

Practical  applications  of  thermodynamic 
calculations 

uV^l  <&  Ti* Sptnecf  i 1 


Available  thermodynamic  databases 


itc  substances 

- JANAF  (U.S.A) 

- IVTAM  (Russia) 

- SGTE  (Europe) 

- F*A*C*T  (Canada) 

- MALT  (Japan) 

Th««?  database*  all  contain  cnthalpv.  cnfrop> . Gibfe*  vnen*)  anti 
heat  capacity  value*  im  condensed  and  gaseous,  stoiehrometrie 
compound* 

Ik 


Available  thermodynamic  databases 


Useful  References 

to  Computational  Thermodynamics 


- SGTE  - range  of  systems,  compr  ehensive  for  steels 

- COST  507  - light  metal  (AI-,  Mg-,  Ti-5  alloys 

- ThermoTech  - Ni-base,  Ai-base,  Ti-base  alloys 

- NIST,  Metallurgy  Div.  - Ni-base  superalioys,  solders 

- F*A‘C‘T  - Fe-,  A}-,  Cu~.  Pb-,  Zrt-rich  alloys 

- Tohoku  University  - solders,  Cu-base  alloys, 

sulphides  in  steels 


| Th&£&?£  Cs&&b&i>ki  Kansas s a?  Wads  «l.  K.Haek,  Institute  af 

5 SgatefiaUt.  London,  1998 

| kimiMtem  TH& Hn&dyn ami t Dais: 

$ &>x>f h-c'-rr  arid  Davabaamaa  e<-l  P J Sptntvtx.  MRS  8u8«ta«  24  (1$£9) 

j Hq.4 

$ Wark«h&i>  or.  Wiv*r*s:&<*tr  for  ZlvmorH  or:!  Oforr 

| ofSolkiom  CALPHAD  19  (1995}  433 

| Wwksho-p  ©in  Th «if modyrzmic  Modelling  Solution's  &nd 

Allays,  CALPHAD  21  (1997)  155 

Rmggs#r$  VVarkaaop  c>n  A o&iz.zxtons  ri 

of  Thermodymms:  in  Fresco*  fdsdonr><?  md  hto 

AprAsoMtorr  fd  Credit CALPHAD  24  (2000)  55 

Workshop  or  Yhzvfrodyr.&mo  M-rcottirg  <?nd 

Zertscbrrft  fur  Metallkunde  92  (2001)  513 

hx '?  a?  Ttv  < irp/tp  « 
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Need  for  awareness  of  experimental 
problems  in  data  assessment 


The  critical  assessment  of  thermodynamic  data  is  a 
skilled  task  which,  preferably,  requires  experimental 
s experience  as  well  as  experience  in  the  use  of  the 
sophisticated  software  available. 


Because  of  the  very  regrettable  decrease  in  the 
number  of  laboratories  carrying  out  thermodynamic 
measurements,  there  is  a real  danger  of  assessment 
work  being  carried  out  by  talented  scientists  with 
sound  theoretical  skills,  but  ar*  inadequate  firsi- 
band  knowledge  of  the  experimental  'difficulties 
associated  with  the  numbers  they  are  assessing. 


(Hi  . Knliiitlptc*  <sf  xolulyirt.  &01  iin.xi 
dr  Cwmjxwww  activiues,  aiflU) 

HI  Hcct  content,  Cj, 

IP  Enthalpies  of  A,  1 

9 Enthalpies  «i  imusfojoiatioa,  A H 
Tf-tjpcnrturc®  of  ttansfonnation 
' 7 Paaxc  «jus!ihtt»,  J.  >:!  x.  pi  }>” 

'Hix-  Omup 
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Important  considerations  in  thermodynamic 
database  development  for  alloys 


( '<  >S 1 Action  507:  IVIa&sat'CsView  arid  evaluation  onSiermbehemkai 
and  ihennophysk'al  properties  to  provide  a database  for  the 
development  of  new  Itgin  alloys 


| Many  commercial  alloys  contain  several  constituent 
I elements,  e.g.  steels,  superalloys,  Al-base  alloys. 

! The  provision  of  a reliable,  self-consislent  database 
| for  calculations  over  a wide  range  of  compositions 
I requires  careful  preliminary  preparations  with  regard 

I to  data  and  model  compatibility. 

I An  example  of  such  work  is  the  COST  507  database 
I for  light  metal  alloys. 


COST  Action  507:  Light  aiioy  systems 
of  potential  interest 

Kev  svslw; 

Additional  elements 

AopSication 

Af-Mg-Mn-Fe-Si 

CU,  Cr,  Tl,  C 

Cans 

Al-Mg-SFCu 

Fe 

Extrusions 

Al-Zn-Cu-Mg 

Zr,  Cr,  Si 

Aerospace 

Al-Li-Cu-Mg-Zr 

H 

Ail  flames 

Ti-A!-Mo-Zr-Sr*-Nb 

ON  C 

Gas  turbines 

Ti-Ai-V-O-N-C 

Gas  turbines 

. V 

tiW  <!»>*» 

’ '' 

Signatory  Countries: 

Austria.  Belgium,  France.  Finland,  Germany,  Greece,  Italy. 
Netherlands,  Norway,  Portugal,  Spain,  Sweden.  Switzerland, 
United  Kingdom 

Participant:  Baikov  Institute  ot  Metallurgy,  Russia 


Co-ordination  Groups: 

A.  Experimental  thermodynamic  investigations 

B.  Phase  diagram  evaiualions  and  compilations 

C.  Thermodynamic  assessments  and  calculations 

D.  Measurement  and  compilation  of  thermophysicai  data 


COST  507 : y afrix  of  binary  systems  representing 
building  blocks  for  alloys  of  potential  commercial 
interest 
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COST  507:  Application  of  evaluated 
thermodynamic  data  for  light  metal  alloy  systems; 
to  casting  and  heat  treatment  processes 


Key  system  for  investigations 


Partners  in  measurement  and  evaluation  work: 

. InM.  f.  Anorganteche  Chnmte,  Untv  Vienna,  Austria  (acf&lty  & ph. 

Dept  of  Met  and  Materials,  Untv.  Leuven.  Belgium  I phase  ti&grem 

MPI  f.  Watatlforschimg.  Imt.f.  Wertatoff..  5tuttgart.  Germany  tr^lorlmefryl 
* U?hfs.tuhl  f Theoret  Htfettenkunde.  RWTH  Aachen.  Oerm.my 

{xsm*.$(rrtsnt) 

MPI  t Metatlfon*chimg  PM L.  Stuttgart.  Germany 

Centre  de  Thermochlmle  e\  Mfcrocatortm..  Marseille.  France  fcaS«r*m#t*p 
^ Uni  vefs  Made  do  Mlnbo.  Gulmaraea.  Portugal  {activity 


«v;»  a' 


1%:  8p**is<*  <irrnp 


14 


SYSXfcM  « Ti 


- ,,r  “^vinrsts : is~ 


fo  ■»  «jf»  ‘ 


TU-  xpr^er  wi.p 
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<?% 

''fils? 


COST  507 


Bcfcrjfrfjn  l*?rmwht?ro«5ft8  fsnrt 
ttverm-wshyaicai  pr^pe°tte5  to  provide  a 
estates**  w th*  esv**©****** 
tsl  atm  fcyfcl  istktft. 


Therutoehemical  database  tor 
tight  metsl  aiioys 

Volume  2 


Tfe<  < Uv-isp 


Development  of  a database  for  Al-Cu-Mg-Si-Zn  allays 

Data  compatibility 


; Consider  the  example  of  a hypothetical  database  for 

Al-Cu-Mg-Si-Zn  alloys 

I A fuli  thermodynamic  description  for  this  system  requires 
saif-consistont  assessed  data  for 


ta&m&rtrs  (5): 
Binary  aiioys: 

fit!} 

Corn  ary  riir-y--.’ 

net 


At.  Cu,  Mg,  Si.  In 

Al-Cu.  At -Mg  Al-Si.  At-Zn,  Cu-Mg.  Cu-Si,  Cu-Zn, 
Mg-St.  Mg-Zn,  Si-Zn 

Al-Cu-Mg.  Al-Cu-Si,  Al-Cu-Zn,  AI-Mg-5i.  Af-Mg-Zn, 
At-Si-Zn.  Cu-Mg-Sr.  Cu-Mg-Zn.  Cu-Si-Zn.  Mg-Si-Zn 
.Qsiaterey  allays:  Al-C«-Mg-Si,  AhCtr-Mg-Zn,  Al-Cu-Sr-Zn. 

|5}  Ai-Mg-Si-Zn,  Ctt-Mg-Si-Zn 


Tic  >/*msu‘! sr  i ttxeip 


Development  of  a database  for  Ai-Cu- Mg-Si-Zn 
alloys  Model  compatibility 


The  work  of  SGTE  in  database  development 


SGTt  is  4 iwvKcR4owcrt  •***««*  «r «a 
et»  Sf» 

3jan*:>  r v<-S>V!«.  svvi  >>Cc.< 


z; 

SGTE 

Scientific 
Group 
Thermodata 
Europe 

Our  aims  arc: 

♦ TN?  i,*WOt41  WO  WttS&IftKtK  UJ  tftj*  W Vf&tf  CSfUiJ.WK'*  fhW  m#*** 

<***  •<-  ■*°**?<*>*  C-CT'STXM  <-5>  l»5  JywMfrtvZT  «.*t  Z-*V*y  <ZXUcC*«iSy 

* USOSTKK*  <d8firtoa»C*-<Js  tftat  <i>  s*>#y  tsemssavArm;,  «Ur»  ur*S 


• TXg  5CE  6iJ8t-5S*8  w\»  iJsimwt 

&S2JCC? 


SGTE  Member  Organizations 


SGTE  application-oriented  databases 
currently  in  preparation 


Trainee 

INPG  - LTPCM  fOnwrobto) 
Thftrm&ddta  {Grawiobte} 

IRSID  (R8at2iftf6£*lfts.WettE} 

Untv  da  Paris  Sud  (Chatenay) 

Sweden 

Royat  Inst  Tecfmol.  $Sto<ikh<rfm} 
TherrooCalc  AB  ^Stockholm) 


ill!  HUllMAack  ca> 
Max-Ptank-lnsL-PML  {Stuttgart} 
GIT  {Herzoganrath} 

United  Kingdom 

AEA  Technology  {Harwell) 

Natl.  Physical  Lab.  {Teddington) 


»t  process 

C.maxL*  ThermFact  {Montreal)  ?;*$>>&  Tohoku  tfotv.  - Mat.  Sci.  & Eng.  {Sendai) 
I >,;  NIST.  Metall.  D*v.  {Gaithersburg} 

- 'TvS,  C..'  \ \w  i j»>r:gr  > 


Bmms&  s«  any  u sers 
sarpiSrsmsms  far  fhgrmadynas 
sufcsat*  of  tire  SGT?  SoMtien 
feefPB  gP5?ared  fsr  ysfesssie 

- Aluminium  aiioys 

- Nickel-base  aiioys 

- Solders 

- Titanium  alloys 

- Copper  aiioys 


rsic  ;iv;c,  the  feUoWlftcj 
Oatahasa  are 


Noble  metal  alloys 
Steels 
Ceramics 
Hard  metals 


Us*  ivjvc<.vy 
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LANDOLT- 
BORN  STEIN 


i *****  ' 'v  : 

! Thenr*ody»amlc 
Properties  of 
laurgaaic  Maleriafe 

$Crrt' 

mm to&M. 


’ 

// 

THEStSTE  CASEBOOK 

I 

Thermodynamics 
at  Work 

Edtced  by  Dr  K>  Hack 

1 

.;•••  ■ 

«)A  ii  ii 

i >««««*»<  t*#m 

- 

&-v\ 

Difficulties  encountered  in  obtaining 
funding  for  database  development 


Applications  in  materials  design  and  development 


Assessment  work  requires  costly  «*pefbse  and  can  be  vasty 
time-consuming 

Database  development  requires  assessment  of  systems  other 
than  those  of  obvious  direct  relevance  to  a practical  problem 
With  increased  Industrial  funding  tor  a project  being  a 
requirement  of  government  agencies,  the  resulting 

corifitfemtaiity  restriction^  are  mete-ssinrsiy  limiting  genera! 
availability  of  assessed  data 

Database  development  is  greatly  assisted  by  use  of  already 
available  assessed  data  for  relevant  sub-systems.  However, 

the  confidentiality  reguire merits  resulting  from  industrial 
funding  are  leading  to  costly  eSispNcadon  or  assessment  work. 

with  resulting  danger  of  incompatibility. 


The  t >ftM fp 


design  and  development  of  new  stainless  steel  grades 
including  prediction  of  amounts  and  formation  temperatures 
foi  carbide  and  nitride  phases 

Light  metal  allots 

prediction  of  solidification  ranges  and  phase  formation  in 
casting  processes 

Coatings 

prediction  of  comjxtsitwm  end  temperature  ranges  for 
formation  of  metastabte  costing  phases  m PVD  processes 


ign  and  development  of  new  stainless  steel  grades 


Design  and  development  of  new  stainless  steel  grades 


- Stainless  steel  production  continues  to  grow  by  about  3 to 
5%  per  annum. 

. As  a result  of  customer  demand  for  specialized  properties, 
and  the  need  to  reduce  costs  associated  wit!)  expensive 
alloying  elements,  new  stainless  steel  grades  are  being 
developed 

- One  development  dire chon  involves  combining  the 
properties  of  ferritic  and  austenitic  steels  in  an  optimum 

way  in  the  well-know  duplex  stain-rose  sleols.  Another 
direction  aims  to  replace  .nickel  by  manganese  and/or 
nitrogen. 


.<•»<'  tv^rtcitr  ‘>-T-vp 


Before  such  new  steel  grades  can  be 
produced  and  marketed,  their  mierostructura 
properties  must  be  optimized.  This  allows 
definition  of  suitable  conditions  for  casting, 
hot-rolling,  best  treatment,  cold  rolling,  etc. 
The  desired  properties  of  the  steels  can 
thereby  be  produced. 


jCjviiOf  txRJhjS 
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Design  and  development  of  new  stainless  steel  grades 

Existing  grades; 

Ferritic:  X2  Cr  11 

Austenitic:  X5  CrNi  18  10 

Under  development: 

Duplex:  X2  CrNiMoN  22  5 3 

Austenitic:  X5  CrNiMo  1712  2 

Proposed: 

Austenitic:  X5  Cri^nN  17  12 

,'}'#■  5¥.  ' ' ?"  ■ ' '">  ..  %/ 

(SS  2 i tC  Tk< 


&at«^s?e<t  pwawt  f&rmxteo  m wa  $&$*«*«&  seaialros 


Design  and  development  of  new  stainless  steel  grades 


Calculations  of  phase  constitution  using  assessed  thermo- 
dynamic data  provide  reliable  guidelines  tor  selecting 
suitable  compositions  for  new  stainless  steel  grades. 

Calculated  liquidus  and  solidus  temperatures,  formation 
temperatures  for  precipitated  phases,  and  composition  and 
temperature  ranges  of  stability  of  deleterious  phases,  enable 
time  and  expense  to  be  saved  in  establishing  suitable 
conditions  for  casting,  hot-rolling,  heat  treatment,  and  other 
operations  for  the  alloys. 


Casting  of  light  metal  alloys 


JS#»«8fSS8$$8M**S8i8S8$8$SS«*«M*8!5S888«4M«iS8S8S^^ 

I Thermodynamic  calculations  have  been  carried 
| out  to 

| provide  the  following  information  as  basts  for 

1 controlled  casting  operations 

f x temperature  range  of  solidification 

| - crystallization  paths  for  equilibrium  and  non-equilibrium 

1 conditions 

1 ■ nature  and  amounts  of  phases  formed  during  solidification 

1 

% «*  heat  capacities  and  enthalpy  effects  associated  with  phase 
| transitions 

I 

| o.vlt  ft."  •«} 

1 

Casting  of  light  metal  alloys 

Calculations  have  been  carried  out  using  the 
following  models:  i 

Equilibrium  cooling 

r Schell  Model  (no  diffusion  in  solid  phases 
formed) 

* Diffusion  in  the  solid  phases 

■ 

:C.?  1V<  S5V»».:.f  <>»vip  4* 
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Casting  of  light  metal  alloys 


Schell  model  applied  to  alloy  AA2024 

The  Scheif  model  has  been  applied  to  simulation  of 
the  solidification  of  the  commercial  AI-Cu-Mg-base 
alloy  AA2024.  This  alloy  contains  the  additional 
alloying  elements  Fe,  Mn  and  Si. 

Important  missing  data,  e.g.  for  the  system  AI*Cu*Mn, 
were  evaluated  as  part  of  the  work, 


Schei!  simulation  for  Alloy  AA2024 


• tvo  diffusion  sis  the  sols!  phases 
f Homc>*j»ne<S!JS  iitjtfki 

• Negiigiiite  unosreootmg 

« No  kinetic  3nd  curvature  effects 

• PiaSDlixc  dendrite  rncrohoioqv 


omparison  of  Schsil  simulations  with  experimental  results 
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Prediction  of  composition  ranges  for  metastahie 
coating  structures  produced  in  PVD  processes 


Certain  coatings  produced  by  PVD  techniques  display  superior 
properties  fot  wear  and  corrosion  resistance  applications.  This  is 
leading  to  the  systematic  development  ot  new  high-performance 
coating  materials. 

PVD  techniques  allow  required  coating  compositions  to  be 
achieved  relatively  easily.  Variation  of  the  substrate  temperature 
may  afso  allow  the  coating  structure,  and  hence  its  properties,  to 
be  changed. 

Thermodynamic  calculations  allow  prediction  cl  conditions  for 
the  appearance  of  passible  arMastalrte  phases  btsfing  WB 
coating  processes  and  thereby  assist  in  the  selection  df  optimum 
coating  pawsssfers  to  pnsdooe  cast tags  with  desired  |tw#orties. 


Prediction  of  composition  ranges  for  metastable 
coating  structures  produced  in  PVD  processes 


I (N.Saunders  and  A.P.  Hftlodownik.  J.Matar  Res.  (13SS>  1 38J 

I "If  an  alloy,  that  in  equilibrium  contains  a multiphase  structure,  is 
f Co-deposited  at  low  enough  temperatures,  the  surface  mobility  is 
l insufficient  fot  the  breakdown  of  the  initially  fully  intermixed 
i depositing  atoms,  and  the  film  is  therefore  constrained  to  be  a 
!pi«*phss«  structure.  This  structure  should  reflect  the  most 
| energetically  stable  single-phase  form  available  ?o  it  at  the 


s This  phase  can  be  directly  found  from  the  Gibbs  energy  versus 
i composition  {G lx)  diagrams  of  the  ailoy  system  of  interest...” 


life  <>tvU2 


Prediction  of  composition  ranges  for  metastable 
coating  structures  in  the  (Ti.AI)N  system 


Prediction  of  composition  ranges  for  metastable 
coating  structures  in  the  (Ti,AI)N  system 


Metastable  (Ti,AE)N  coatings  with  the  cubic  NaCI 
Structure  are  superior  to  TiN  for  cutting  tool 
applications  due  to  their  superior  oxidation 
resistance  and  hardness  at  elevated  temperatures. 

However,  the  deposition  of  |Ti,Ai)N  films  with  Ti/AI 
ratios  below  about  30/70  leads  to  a hexagonal 
coating  structure  which  is  not  suitable  for 
tribological  coatings. 


Thermodynamic  calculations  allow  the  ranges  of 
composition  and  temperature  in  which  the  metastable 
cubic  and  hexagonal  coatings  are  formed  to  be 
predicted. 

Of  prime  importance  for  the  calculations  are  the  energy 
differences  for  the  transitions  from  the  stable  to  the 
metastable  structures,  i.e.  TiN,.— >TiNh„  and 
AIN^^AIN,^. 

Since  no  experimental  data  are  available  for  these 
transformation  energies,  a combination  of  calculation 
and  estimation  techniques  must  be  used  to  obtain  them. 


Prediction  of  energies  of  transformation  between  the  cubic 
NaCI  and  hexagonal  structures  for  transition  metal  nitrides 


Thermodynamic  calculations  for  (Ti,AI)N  coatings 


<8  Hi  i&  *:  W 53  S4  % 8 tt  8?  £5  85  8*  «•  6J  SS  88  «8  8J  f 


Yfc  Ntr-inver  > '.r<  up 


A thermodynamic  analysis  of  the  (Ti,AI)N  system  was  made, 
which  combined  calorimetric  masses  arrsaMs  of  the 
metastable  to  stable  phase  transition  tor  3 compositions  in 
the  section  AIN-TiN  with  estimated  vaiuasfot  Gibbs  energy 
differences  between  the  NaCi  and  Zr  S wurtsste  structural 
forms  of  pure  TiN  and  AIN  This  allowed  Gibbs  energy 
curves  to  be  calculated  for  the  hexagonal  and  cubic  phases 
for  various  temperature  in  the  AIN-TiN  system 


The  aernt  of  intersection  of  these  curves  at  each 
temperature  defines  the  composition  at  which  there  is  a 
transition  from  one  structure  to  the  other.  This  composition 

was  round  to  be  nearly  temperature  independent  with  s 
calculated  value  around  0.7  mole  fraction  AIN 
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Gibbs  energy  of  formation  of  the  cubic  and  hexagonal 
phases  in  the  AlN-TiN  system  at  different  temperatures 


Experimental  measurements  for  (Ts,AI)N  coatings 


I Experimental  studies  of  the  extent  of  the  metaslable  cubic 
I range  in  the  section  AIN-TiN  have  been  carried  out  by  Cromer  et 
| sir  using  electron  probe  micro-analysis,  X-ray  photoelectron 
I spectroscopy  and  thin  him  X-ray  diffraction. 

| ■ Their  study  showed  that  tor  co-deposition  temperatures 
N Sivcee  300  and  T88  ‘C  a two-phase  region  exists  in  the  roeta- 
| *tahl«  A IN -TOO  phase  theorem.  This  has  a composition  around 
I x{AIN;  = 0,86  at  lower  temperatures,  but  its  boundaries  expand  to 
Sower  and  higher  AIN  concns,  with  increasing  temperature.  As 
ip  VC  substrate  wi/u-Mf  no  two-phase  rt<v>rs  was  observed. 

'14 r rerun’  MPrthsit  DNwxhiitl  m ta  A r-div.  Sfe  i'O'Jj  bits  w& bboard 
KiVScfys  The M/ieiA-.  Metafe;A Sodsry  {te&hroAIA 


Experimental  measurements  for  (TlAIJN  coatings 


Prediction  of  composition  ranges  for  metastable 
coating  structures  in  the  Al203-AiN  system 


TIpx  f>} Ksicsy  wujj 


There  is  no  reciprocal  solubility  below  -1627K  in  the  solid  state 
between  AI.Q,  and  AIN.  However,  a metastable,  homogeneous 
oxynitride  coating  piiase  coaid  have  potential  applications  in 
many  areas  because,  e g . 

* the  refractive  index  of  a pains  ASH  film  Is  high  whereas  die 
sefr active  index  of  At  ,0.  is  tew 


the  thermal  r.osteuafivity  of  AiN  is 
<ii  AijO., 

is  a semte ondaetdc . Af.,0 , 


comparer!  t»  that 

3 '".'A1 

insulator 


If  Al-O-N  phases  with  variable  composition  can  be  stabilized  by 
vapour  deposition,  it  may  be  possible  to  combine  the  different 
properties  of  the  binary  compounds. 

riV2J  ft?  TV  Xf&xxae  iiv'/ttp  ;X 


Experimental  measurements  for  A12Q,-A!N  coatings 


Thermodynamic  calculations  for  Ai^O^-AiN  coatings 


filetastabie  At-O-N  layers  have  been  produced  by  reactive  magnetron 
sputtering  ion  plating  using  an  At  target  and  an  Af-CytV.  gas  mixture/. 

For  a substrate  temperature  of  ISO  C and  variable  partial  pressures 
of  Cb  and  N3.  a variety  of  coatings  with  different  0 and  N contents 
worn  prepared  and  analysed  by  means  of  X-ray  photoeleetrori 
spectroscopy,  Auger  electron  spectroscopy,  high-resolution 
transmission  electron  microscopy  and  high-rasolution  scanning 
electron  microscopy. 

C;.yosr.-/ch  At-Q-N  phsso.  were  found  to  be  nanooryBlatiine  with  a 

cubic  gamma-A|./3,  {cpiael)  structure.  U-rasg-cn-rch  p'n&tiesa 
displayed  a vsonenai  Aifr  u.vurttitei  structure.  For  O ff  ratios,  around 

t 3,  the  deposited  phase  was  TEfS  . 


■a  vcr.  /sreoler.  / Pcow.!.  /.tn  Sst-t  FiOc  /r'A;  tttt  At 


. Its  >pes;r-  i'mtp 


An  available  thermodynamic  assessment  of  the  AI2Oj-AfN 
system'  was  amended  to  incorporate  descriptions  for  tee 
wurtzite  and  amorphous  phases  for  the  entire  composition 
range.  Using  these  data,  GiW»s  energy  curves  were 
calculated  tot  the  different  matastafele  phases  for 
tathporatees.  aetwoon  room  tempo: aiute  site  AOO-  K.  The 
experimentally  observed  composition  ranges  tor  the 
phases  at  463 K were  used  to  help  define  the  location  of  the 
Gibbs  energy  curve  for  the  amorphous  phase 

> The  metastable  phase  diagram  was  calculated  and 
comparison  made  with  experimental  results. 

-and  5S  Jorwsatt  Z Mytfflk&t.  $%  ??4 
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Gibbs  energy  of  formation  of  the  spinel,  amorphous 

and  wiirtefte phases  in  die  Al,OrAiN  system  at  463K 


Calculated  metastebfc  phase  ranges  in  the  A1..0  .-AIN  system 


fcrf  If&ctfOn  AIM 


Calculated  metastable  phase  ranges  in  the  Ti.B  :-TiC  system 


Conclusions 


- Thermodynamic  calculations  ate  able  to  provide 
information  of  significant  practical  importance  on  tte 
composition  and  temperature  ranges  of  metastable 
coating  phases. 

Estimation  or  calculation  techniques  are  required  to 
provide  transformation  energies  for  metastable  structures 
of  potential  coating  compounds  (carbides,  nitrides, 
borides,  oxides,  etc.},  but  the  present  work  shows  that 

with  limited  experimental  support,  satisfactory 
psedlcrtaes  t>?  metasfafete  phase  ranges  can  te  made. 


Summary 


There  have  been  rapid  advances  in  software 
packages  for  simulation  of  the  thermodynamics  and 
kinetics  of  materials  development  and  processing. 

These  advances  must  now  he  matched  through 

increased  support  for  the  -measurement  and  critical 
■mm&mrmtft  of  "missing  data”  required  for  the 
simulations. 

Database  development  is  essential  for  full 
implementaton  of  thermodynamic,  eafetiteisoji 
-potential  in  many  areas' of  materials  technology. 
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Additional  points  made  by  speaker: 

Needs  for  database  construction: 

Experimental  experience:  independent,  first  hand  knowledge  of  experimental  difficulties: 

• Mn-Ti  phase  diagram  used  to  illustrate  some  of  the  difficulties  that  may  be  encountered  in 
experimental  work  and  their  effects  on  the  thermodynamic  assessment. 

• Zr-B  shows  temperature  differences  up  to  500  K between  data  sources. 

Self-consistent  databases  require  careful  preparation:  Example  COST  507 

• Number  of  components  very  large 

♦ Need  to  break  down  into  smaller  tasks  focussing  on  important  subsystems:  Example  Al-Cu-Mg- 
Si-Zn 

• Modeling  of  the  intermediate  phases  and  the  subsystems  must  be  compatible 
SGTE  and  its  aims: 

• Identifies  focus  areas  for  database  development 

• Landolt-Bornstein  publications:  unary  and  binary  systems 

♦ Funding  problem:  experts  are  required,  time  consuming. 

Predictions  for  PVD  (physical  vapor  deposition)  processes: 

Metastable  phases  are  important,  since  phase  which  has  lowest  Gibbs  energy  at  any  composition  will 
form  first. 
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Application  of  Thermodynamics  to  Alloy  Design  and  Development 
Charles  J.  Kuehmann,  QuesTek  Innovations  LLC,  Evanston,  IL 


Application  of  Thermodynamics  to  Alloy 
Design  and  Development 

Charles  J.  Kuehmann 
QuesTek  Innovations  LLC 
Evanston,  Illinois 
March  21  2002 


pUEsTe  Ht 

INNOVATIONS  t-LC 


by  Dtsiffii 


PROCESSING 


UHS  STEELS 
STRUCTURE 


EBflEEBIffiS 


Hierarchy  of 
Design  Models 


Process  - Structure  Modelma 


TC<CohJlQ!CTRA  - (C-  j ppT.H 
ABAQUS'EFG 


FlAPWl 

loVM 


CpUiKsTHs  W 

iSSQVATiOHS  U.S 


***$»$*  B&ti&t* 


INNOVATIONS  LUG 


How  CMD  Programs  are  Used? 


Available  Multicomponent  CMD  Models  for 
Ferrous  Alloys 


I1  mcess-St ructit re  Models 

Structure-Pronertc  Models 

❖ Displume  I run  sformn  lion 

❖ Basic  Physical  Properties 

* Martensite  Ms.  tM(  1) 

• Liquid  Density  and  Solid  Lattice 

• Bainile  B,.  T I1VCC I' 

Parameter  and  Thermal  Expansion 

❖ Tempering  Microsimciurc  — 

■ BGC.  R C Shear  Moduli 

■ B<  (.  l'e3t  Paracouilihritsm 

■ Electrical  Resistivity 

• (Carbides.  Inter  metal  lies  and  T'CC 

• Thermal  Conduct!  vh\ 

Dispersion  Formation  Models 

❖ Strength  Model  including  Solid 

❖ Precipitation  Kinetics 

Solution,  (drain  Size,  Precipitation 

* PrecipiCak1^  3D  General 

(Shearing.  Bypassing)  Strengthening, 

Precipitation  Program 

and  G TS  Y8  Prediction 

* Coarsening  Rate  Model 

❖ Toughness  Model 

v Microsegregation  (Schell)  and 

❖ (train  Boundary  Cohesion 

Macrosearegafion  (buoyancy) 

❖ Multicomponent  C orrosion  Resistance 

❖ Grain  Growth  Mixiel.  Grain  Refining 

Model 

Dis|x.rsion.  Homogenization  and 
Austenitization  Designs 
❖ Grain  Boundary  Segregation 

Note:  not  all  of  these  models  arc 
implemented 

INNOVATIONS  LLC 
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CpUEsTEH 
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M* uerixte  by  Ili’fii'fsy* 
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More  About  TCIPC 


1 scr  Operation  Environment 
{e.q.,CUI,  iSlCiHT) 


1_L 


\ ser  Application  Program  {('  ) 

Including  Other  Numerical  C<xles 

— t — t — i— 

H IK 

(Implcmciiit'd  in  ANSI  ( ( ■ i is  jth 

I M\ Si. stem  Calls 


Themio-Culc  Thema'A  ak 


As  < Ml)  programs  currently  use 
TCIPC  ti>  incorporate'  Thermo- 
Cnlc/DK  IRA,  this  implies: 

* Program  can  not  run  without 
Thermo-Calc  and  DICTRA 

• User  need  to  have  proper 
privilege  and  setup  for  running 
Thermo-Calc  and  OJCTRA  in 
order  to  run  ( Ml)  programs 

u It)  tc  r-  P roc  css  ( cm  mutt  icalion 
T Intra-Process  Communication 


R :i;iviin>2  Ditk'ie.'i!  tVocCsSCS 
Ictruinnied  when  User  P-oamm 
I i j Visaed 


CpUEsTeK 
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Materials  by  Des>ig 


CMD  User  Interface  — Pre-Processor 


CpuesTen 
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Matt  mi*  by  Dcnigri 


CMD  User  Interface  — Post-Processor 


PrecipiCalc™  Input/Output 


t*»  t*a 

N-fttae  Pfcii  ' --lev** 

k 5**1  •••»«•>{  {ftoavs# 


w,i».  * • tv  S*v»: 


P)»f  Aiuw 
Sc  aUn  •'.Shiftittg 
Itrsnliv. 

Tuuisf<>rma$i»H 


<»>  UxU  Noany  rtaH  » 

fat  O 

(itV|jt»tt)*w'H> 
|t|r(Mti»«|) 

fntt.r<nt«>oj  o 

n»>«fsooa.nno 
p O «nv»«s  o 


Input  — 

• Alloy  (ompo.sitkm 

• Prior  Heat  Treatment 

Condition 

"'" t 

• Thermal  Cycle 

wvj 

• Fundamental  <>ttamiiies 

f l.iuik  and  SuiTiicc 

Thermodynamics,  Lattice 

Parameters.  Mobiliiy.etc.) 

ttiiPsii 

• NimicHcal  Paramek'rs 

Output  — 

Particle  Size 
Distribution  versus 
Time  ( Mean  Size, 
Volume  Fraction,  and 
Number  Density) 

Compositions  of  Matrix 
and  Precipitate 
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Material*  by  Design 


IpLJErsTet-! 
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Alum  till*  by  Destig. tf’ 


Modeling  Structure  of  PrecipiCalc™  Software 

Dynamics 


Basic  Multicomponent  PrecipiCalc™  Equations 

J JR) 


ih  transient  nuclcanot;  Uncub.ihOJ-  lime  r ) 


Cofumuiiv  — — 

a SH 

JK  2*n\ 

th  kl  I ""  if  | 


r ~ P4L rJ-iPuT - j«rf  ~ -•> <?. • c:> j 

I,  '-  f j 

[v  “f  «|<r?  ••<?> ~c?)ja»d|*{J  A f'~! (it'  - < “'}] 

with '<>  i,  .n and  >.  ;.k  -•  2 u e>,  ptxK&i. v omtvxient  1 in balsoced 

Mass  fiaiancc  if  - — - 1 JTl'l.'f  (R)tl ft  t I - '/■)  where  tj>  - — i ftPtlR 


]i  pnrticlc  rudum  / •— iiftic 

JiKO  particle  frequency  iUstritottion  (l-rcr1) 

:V,  u>taJ  number  density  of  particle  f ! th1) 
fi  — - chemical  potential  C Gibbs  fret-  energy 
i \ solute  coiicemrsiiojj  f\ . — • d-ITusivity  main* 
j?-— surface  energy  • „ Jic.*ku  volume 
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Future  CMD  Extension  — Integration  with 
iSIGHT/DEFORM  to  Achieve  Component-Level 
Materials  and  Processes  Design 


aimSIGI  IT 


Process  Integration  an 

d Design  Optimization 

1 1 

DEFORM™-HT 

CMI)  Software  System 

TlKMinal- Mechanical  Finite 

• Process-Structure  Models 

Element  Component  Level 

• Structure-Property  Models 

Simulation 

PrecipiCale1'9*  software,  a component  of  CMD  system,  has  been 
successfully  integrated  with  this  manner  in  DARPA  AIM  (Accelerated 
Insertion  of  Materials)  initiatives  for  Ni- based  sttperailov. 

Cp  LJ  E:  sTfe  M Materials  by  Design* 

INNOVATIONS  LLC 


The  Materials  Genome  Initiative 

• Priorities 

- Thermodynamic  Databases 
— Mobility  extensions 

- Mot ar  Volume 

• Treatments  for  accuracy  and  variation  and  their 
propagation 

• Efficient  database  development  methods 

- New  experimental  techniques 

- Optimal  integration  ot  quantum  mechanics 


Cp  LJ  EE  sTe  W Mateads  tn  Design- 

INNOVATIONS  L.LC 


Additional  points  made  by  speaker: 

Performance  properties  <h>  structure  processing 
Modeling  employs  thermodynamic  quantities. 

Computational  materials  design: 

Entire  suite  of  software  is  used  depending  on  dimension  level  (atomic,  nano,  micro,  ...)  to  model 
properties  =>  "faster  science"  & "better  engineering" 

Software  includes  a series  of  process- structure  and  structure-properties  models. 

Software  has  no  direct  interface  to  Thermo-Calc  - runs  as  independent  process. 

PrecipiCalc:  uses  nucleation  conditions,  growth  law,  full  diffusivity  matrix  and  fitting  parameters  of  the 
precipitate  surface  energy  to  predict  particle  size  distribution  and  composition  of  precipitates. 


106 


Workshop  Presentations 


Construction  and  Application  of  a Diffusion  Mobility  Database  for  Ni-Base  Superalloys 
Carelyn  E.  Campbell,  NIST,  Gaithersburg,  MD 


Hie  Caristntc&m  amiA^aMai  of  a 
DiffimonMdM^D€Aix^ 
jbrNi-Base  Supetx Map 


t\JtST 


C.  E.  Campbell 

National  Institute  of  Standards  and  Technology 
Metallurgy  Division 
Gaithersburg,  Maryland  20899-8555 
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y precipitation  during 
non-isothermal  cooling 


Heat  treatment  Protective  coating 
optimization  applications 


c 

.§ 

V J 

r 

! 

\ f 

V 

Cooling 


v+Y 


Distance  across  denckrte 

'J|p  Heating  risers  difusicn  during  extended 

service  times  at  high  temperatures. 


Distance  across  denckrte 


MIST 


Diffusion.  in  aBmanjABoy 


»«  »«  © 

•••+••«««  »©© 
a® 


| Heat 


Position,  * 


Pick's  first  law  for  Flux,  J 

J = -D± 
dx 

Fick's  second  law  (conservation  of  mass) 

ac=_a_f'Z)ac 

dt  dx  1 dx 

S?  D is  independent  of  composition 

dt  dx- 

ff  Semi-infinite  sample.  Error  function  solution, 


cix.t)  = A+Beri\  , ' 
l V4Dr 


rstlST 


Temperature  dependence 


r Inputs: 

- Thermodynamics 

- Diffusion  experiments  (unary,  binary,  ternary  systems) 

* Tracer  diffusivity, 

* Intrinsic  diffusivity, 

» Interdiffusion  coefficients 

>'  Output:  Diffusion  mobilities 
~ Composition-dependent 
— Temperature-dependent 

'S  Validation 

- Diffusion  correlations  (Activation  energies) 

- Extrapolation  to  higher  order  system 

• Diffusion  coefficients 

• Multicomponent  diffusion  profiles 


ch 


Experimental  phase  diagram.  Determine  Gibbs 
Thermochemical  data  G = f(x,T,P) 


Calculated  phase 
diagram 


Liquid  4 

V l i t<Te  ; 

Liquid 

e 

J. __ & P 

1 

3 

! 1 

cj  ' e 

1 

/ \ P 

A 

CD 

i i 

B 

Assume  all  partial  molar  volumes  are  equal 
Assume  a vacancy  diffusion  mechanism 


• Fluxinl-D  _„-gr -t 

J1  —■  --  I x 
£1  3: 

* Interdiffusion  coefficient  (chemical  diffusivity) 

Dj  = ZA  — 


® Intrinsic  diffusivity 


D>)  = ^ jSp'  —x  jpMp 


djuf 

dxJ 


y — i — 0 —ideal  —excess 

U = U + (j  + Lr 


Binaries  lip  Ternaries*p>  Quaternaries  I 


Mobility 


nth  order 
systems 


M'  - ^ exp|  ^ | whereAf?  = , 7’j  and  AT'  = f (a, T j 


NET 


NIST 


A Q = ^ xpQ  + , 


0 pq 
- XPX1  Q 


Workshop  Presentations 


107 


A^essnwntqfDiffijsonM(M!M^ 


[Estimate  Mobility^ 


Compare  experimental 
and  calculated  D 


Simulate  diffusion  process 


4 [Adjust  Mobility] — *1 


Experimental 
diffusion  data 


Mobility  M=f  (c,T) 


Calculate  diffusion 
Coefficients  D = f(c,T) 


1 

'• T a 1050  "C 

T-850-C  ' 

1 OS 

& 

S " 

atukon  pc  at' to 

v „ -*Bflu*lonCo«flttant 

X 

\ 

\ 

w 

WIST 


M,  = cxpl — 1 v.  here  AG,  = f(c,,T)  and  M’  =1 

RT  ^ RT  I . . 

For  a binary:  G*  =c,Q'.+ciQl  1 He, -c,)#-’  + (c, -c,)2C,'J  +...) 


Ezmmpies  of  Fit; ferr  Binary  Interactions 

Ni-M-  Cr-C&Hf-Nb-Mo-Re-  Tcs-TS-  W 


Ni-Co  Interdiffusion 


Interdiffusion  with  Ni 


mist 


Atomic  Percent  Ni  Weight  Fraction 

Previous  assessments  Ni-AI-Cr  Engstrom  and  Agren.  Z Metallkd.  87  (1996)  92. 

Ni-AI-Ti  Matan  et  aJ.,  Acta  mater.,  46  (1998)  4587. 

Current  assessments:  NhCo,  Ni-Hf,Ni-Mot  Ni-Nb,Ni-Re,  Ni-Ta,  Ni-Ti,  Ni-W 
Co-Cr,  Co-Mo 

C.  £ Campbell,  W.  J Boettmger.  U.  R.  KaOner,  Acta  Mat,  SO  (2002)  775 


Assessment  of  M~  W 


Interdiffusion  data 


■JT* s » j s 

W%-rrrf-t^rrrrrrrTT?T 


sTiTiTsTt  t J f 7TTTTTT 


fTT^rrrrrr 

800  “C 


t+1 


)UU 


V -Q 

r For  a pure  metal 


= 17 


Activation  energies  in  the  fee  phase 

ag;  = *mg«  +x».qZ  +*kxw°aZ‘w 

a Ql  = xmQ£‘  +xwQ*  +xmxvl°A^w 

Self  activation  energies 

, Optimized  parameters  . 


Element 

Crystal 

Structure 

Ty,  K 

TM,  K 
(fee) 

Ty,  K (fee) 
(Kaufman) 

Activation 

Energy 

(Jtinole) 

-OmM  (tec) 
(SGTE) 

-QfRTy  (fee) 
(Kaufman) 

Ni 

fee 

1728 

1728 

1725 

-287000 

20.0 

20.0 

Al 

fee 

933.5 

933.5 

931 

-142000 

18.3 

18.4 

Cr 

bcc 

2133 

1475 

860 

-235000 

19.2 

32.9 

Co 

hep 

1770 

1768 

1768 

-286175 

19.5 

19.5 

Hf 

hep 

2504 

1952 

2076 

-235350 

14.5 

13.6 

Mo 

bcc 

2895 

1740 

1530 

-254975 

17.6 

20.0 

Nb 

bcc 

2468 

1300 

1170 

-274328 

25.4 

28.2 

Re 

hep 

3459 

3084 

2830 

-382950 

14.9 

18.3 

Ta 

bcc 

3296 

1416 

1540 

-268253 

22.8 

20.9 

Ti 

hep 

1946 

900 

1421 

-256900 

34.3 

21.7 

W 

bcc 

3695 

2229 

2230 

-311420 

16.8 

16.8 

Mtsr 


Conpcraan  uMiNi-CxxCr-Molkms.al.  130D°C' 


Diffusion 
coefficients 
calculated  using  tt 

thermodynamic 


Composition  (Atomic  Percent) 

Ni  = balance 

Djfxl(t4m2/s 

Cr 

Co 

Mo 

Measured 

Calculated  j 

NIST  Thermotech 

24.2 

24.1 

7.4 

75±15 

10.2  L 

10.7  [ 

22.7 

24.5 

7.4 

9.7  ± 1.9 

10.’. 

io.e 

20.8 

25.0 

7.4 

95  ± 20 

9.35 

10.3 

18  4 

25.6 

72 

10.1  ±Z0 

9J56 

10.0 

15.2 

25.8 

7.4 

82  ± 1.6 

9.35 

9.74 

io.e 

26.2 

7.4 

65  ± 1.4 

855 

925 

6.4 

27.1 

7.7 

6.4  ± 1.3 

8.4 

859 

3.2 

47.9 

7.7 

65±  14 

454 

5.03 

Ni 

26.8 

1.7 

6.6 

8.9  ± 1.8 

10.3  1 

10.3 

26.5 

4A 

6.4 

6.0  ± 1.2 

9.61 

9.78 

26.3 

7 A 

6.6 

4.8  ±1 .0 

856 

922 

25.8 

19.8 

7.1 

3.7  ±0.7 

7.01 

753 

25.8 

21.4 

7.1 

42  ±0.8 

653 

758 

25.9 

13.2 

7.1 

3.3  ± 0.7 

7.47 

753 

-Co 

22.2 

3.7 

6 2 

-20  ± 0.4 

IS  1 

-4.69 

65 

23.9 

7.6 

-1.7  ±0.3 

-2.37 

-2.27 

R^M3s(rtI)D^imaiM€inx^I^B0C 


/Ri 


ene-N4  (x1&,4m2/s) 


+ 1195 
-1137 


+ 13.93  + 34.83  + 3434  + 42.43  + 5150  + 4951  +53.22 
+17.00  -8.25  -5.67  - 555  -133  -7.10  -9.69 


VL 


-426  -537  +13.67  -331  +8.93  +9.91 

-8.33  —0.280  -0.426  +7.57  -055  -036 

+031  +0.25  + 0.66  + 0.27  + 24.05  + 0.74 

-0.68  + 0.33  + 053  + 0.24  + 0.26  + 0.76 

+ 1.63  +1.35  +4.94  +4.94  +6.25  +657 

■aLS 


+ 2.49 
-0.45 


+8.25 
-0.17 
+ 0.85  +031 

+ 0.50  +0.23 

1-23.62  + 5.41  / 

i&S 


Ni  = dependent  variable 


Rene-N5  (x10-14m2/s) 


A 


Al 


Co 


Cr 


Hf 


Mo 


Re 


To 


W 


Mtsr 


Al  +93.16  +13.92  + 33.46  - 651  +33.42  25.44  + 48.63  + 50.87 
Co  -651  +27.22  -856  - 27.64  -4.95  -5.11  +3.87  -9.21 

Cr  +4.15  -4.23  +21.02  -635  -032  -0.78  +13.81  +6.89 

Hf  0.S6  +0.07  +1.70  +262.1  +152  0.87  +237  +1.84 

Mo  —0.35  - 0.30  -030  -1.905  +7.71  -035  -0.13  -0.19 

Re  -0.75  -032  -036  -259  -035  +0.08  -051  -032 

0.03  +033  + 0.98  -4.17  + 0.64  + 0.86  +7.75  + 0.87 
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Rene-N4/Rene-N5  d I2SB  °CJbr  lOOh 




Databases  used 

- Thermodynamics: 
Thermotech 

- Diffusion  mobilities: 
NIST  Ni-mob 


Rene-N4 

Rene-N5 

Y 

Y 

6.35  mm  6.35  mm 
Double  geometric  grid:  200  points 


f\asr 


Expen  mental  work  performed  by  T.  Hansen,  P Me  re  wether,  B.  Mueller, 
Howm^  Corporation,  Whitehall,  Ml. 


Hovwnst  £3 


Rene-N4/Reri&N5  Diffusion  Couple  ct  15SB  °C 


Cr  profile 


Al  profile 


Rene-N4 


Rene-N5 


Howroe!  2 
Castings  “»> 


Ren&N4fflen&N5DiffiGkxi  Cotffe  ct  22S8  °C 


Ren&N4men&N5D$uaon  Goupie  ct 1293 °C 

■ ■ • • '•'•'"•'■•'"r— — r — 


Ti  profile 


Ta  profile 


Re  profile 


W profile 


0065 


7 o 


-Tw*  lOh 
— Tine  = 10C  h 


Rene-N4 


Rene-N5 


INJ1ST 


Hoamnet  0 
Castings 

x-:XyX'XX-»X>:-x 


Rene-N4 


Rene-N5 


Howmet  2 


NIST 


Porosity  Prediction 


Rene-mUN-im  1000hctll90°C 


Back  scatter  image;  100  h 


Distance  (m) 

1 Predicted  position  for  f 

f Maximum  gives  location  of  pore  formation.:  I maxlrnum  P°re  formation  j 

L j 


At  1 1 50  °C  equilibrium 
phase  tractions 

• Rene-88:  ^=1 

• IN-100:  fv=  0.638  fy  = 0.362 


Rene-88 


IN-100 


Howmet  Q 
Castings 


-400  -300  -200  -100  0 100  200  300  400 

Rene-88  Distance  (pm)  IN-100 


NIST 


Experimental  data  from  J.  C.  Zhao,  GE-CRD,  Schenectady.  NY 
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Ren&SmN-im  lQ00h€tll50oC 

"":•  1 


H aatTxeafei  &a.tQ pdxt  i^.tx>n : 
SoWffiisatb&i  ofN  &-11A  1-4  ST  a (afc%  } 


Ntsr 


Rene-88  Distance  (pm)  IN-100 


Solidification  assumptions 


TN 

itions 


* Cooling  rate  - 8K?s; 

* a/2  -95  pm 


* Microstructure  evolution 
j t-.  1 t - Os 

j '-f  ^ I 0*<l«** 


J 


Solidification  Path 


Phase  FmAm,EMxMm 


Distance  (pm) 


• Used  composition  and  phase 
fractions  from  solidification 
calculation. 


* Assume  linear  heating  rates 
beginning  at  800  K 

» Assume  y fraction  is  in 
equilibrium  with  y matrix  at 
each  grid  point. 


• Assume  incipient  melting 
occurs  at  the  center  between 
dendrites. 


KMAJ 


hi 


Ni-11AI-4.5Ta  (at.%) 


: So|i<)us 

0.01  (K/s) 

10  (K/s) 

; 1000  (K/s)  100  (K/s) 

: 

; "^soKrus(Vy) 

te 

0.0001  0.001  0.01  0.1  1 10  100  1000 
Yr  (K/s)-i 


Mtsr 


For  Ni-11AI-4.5Ta  alloy 
solidified  at  8 K/s 

Heat  from  800  K to 
1655 K all  K/s 

> Hold  at  1655  K until  t = 
1 325  s (At  = 470  s) 

grnA  Homogenous 
microstaicture 


MIST 


r Database  Development 

- Ni-base  superalloys: 

• Ni-AJ-Co-Cr-Hf-Mo-Nb-Re-Ta-Ti-W  (Fe) 

• FCC  (y)  phase,  y\  B2 

- Comparisons  with 

• Diffusion  correlations  at  melting  temperature 

• Experimental  quaternary  diffusion  coefficients 

• Multicomponent  diffusion  couples:  (Ren6-N4/Ren£-N5;  Rene-88/1N-100) 

- Experimental  evaluation  of  selected  ternary  and  quaternary  systems 


r Applications 

- Heat  treatment  optimization 

- y precipitation  during  non-isothermal  cooling 

- Coating  deterioration 


isasr 


Ntsr 
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Additional  points  made  by  speaker: 

Simple  criterion  for  validation  of  mobility  database  is  that  the  ratio  of  the  diffusion  activation  energy  of 
a pure  component  to  the  melting  point  of  the  pure  component  is  roughly  constant:  -Q/RTM  =17  . 

Porosity  prediction:  Position  of  maximum  in  vacancy  flux  gives  position  for  maximum  pore  formation. 
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Incorporation  of  CALPHAD  Calculations  into  Phase  Field  Modeling 
William  J.  Boettinger,  NIST,  Gaithersburg,  MD 


Phase  Field  Modeling  using  Thermodynamic 
and  Mobility  Databases 

• Sharp  vs.  Diffuse  Interface  Modeling 

• A few  examples  of  Phase-Field  Modeling 

• Phase  field  variable:  artificial  or  natural  order  parameter? 

• Free  energy  functions 

• Gradient  Thermodynamics 

- Interface  Energy 

• Mobility  Parameters 

- Phase  field 

- Diffusion 

• Addition  of  other  fields 

- Elasticity 

- Electrostatic 


Traditional  Materials  Science  Approach: 
Example:  Solidification 


Tip  Kinetics 
L • Tip  Temp.  < TL 

• Core  composition  > kC0 


Microsegregation 

• Core  composition  > kC0 

• Solid  Diffusion 


Dendrite  Arm  Coarsening 

• Side  arm  shrinkage/growth 

• Side  arm  bridging 


Sharp  vs.  Diffuse  Interface 


mist 


Nisr 


Microsegregation  fWarren  and  Boettinger,  1996) 


Solid 

Liquid 

Sharp  Interface  Model 

C 


Solve  diffusion  equations  in  solid  and  liquid 
Match  solutions  at  interface;  interface  conditions 
Difficult  ‘moving  boundary’  problem 


Diffuse  Interface  Model 

C 


Add  new  variable  and  equation;  (J)  <x,t) 
Solve  new  eqn.  and  modified  diffusion  eqn. 

over  entire  domain 
No  boundary  conditions  at  interface 
Easy  to  program  (new  large  computer) 


Distance 


MIST 


MIST 

Simulation  of  f Precipitation  w/Coherency  Strains 


Binary  Eutectic  with  Dilute  Ternary  Impurity 
M.  Plapp  and  A.  Karma  1999 


Y.  Wang,  D.  Banetjee,  C.C.  Su,  (OSU)  & A.  Khachaturyan,  (Rutgers) 
Acta  Mat.  46  (1998)  2983 
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Active  researchers:  Phase  Field  Modeling 
attempting  “Calphad’Thermodynamic  Databases 


• Y.  Wang,  D.  Banjeree,  C.C.  Su  and  A.G.  Khachaturyan,  Acta 
Mat.  46(  1998)  2983.  Ni-Ni.Al polynomial  4 sublattice  stress 

• U.  Grafe,  B.  Bottger,  J.  Tiaden,  S.G.  Fries,  Mod.  & Sim. 

Mat.  Sci.  & Eng.  8(2000)  871.  TQ  interface Ni-AI-Cr-Co  Ti 
solidification 

• I.  Loginova,  G,  Amberg  , J.  Agren,  Seattle  2001  TMS 

• Miyazaki,  Calphad  25  (2001 ) 231.  Discrete  lattice  Fe-Co-Al 
and  InGaAsP,  w/stress 

• L.Q.  Chen,  C.  Wolverton  V.  Vaithyanathan,  Z.K.  Liu,  MRS 
Bulletin  26(2001)  197 . Al-Cu  0'&  0 ‘ ppt.  First  principles 
free  and  interface  energy  used  to  get  phase  field  paratneters. 


Phase  Field  vs.  Order  Parameter(s) 

•If  thermodynamic  database  treats  the  phases  together 

• e.g.  chemical  miscibility  gap,  order-disorder 

• Natural  order  parameter!  s) 

• Thermo,  benx'een  the  phases  is  not  artificial 

• Can  get  gradient  energy  coefficients  from  bulk  thermo! 

• Need  interatomic  dimension  / interface  coordination  #'s 

• If  thermodynamic  database  treats  the  phases  independently; 

• e.g..  liquid  - solid, 

• Need  to  introduce  artificial  ‘phase  field  parameter,  0 

• Introduce  energy  hump  (W)  between  the  phases  and 
gradient  energy  coefficients  to  get  interface  energy/ 
adsorption  consistent  with  experiments. 


WIST 


Ntsr 


What  is  the  Phase  Field 
variable  for  solidification? 


Rough  Interface 


Gradient  Thermodynamics 


F=J 

J v 


Equilibrium 


f(c,0,T)+^-\V<p  f+^|Vcf 

{ 


dV 


^£=2X-e;V:0=  0 


80  00 


^ ^ ^ - £ c V c = constant 

8 c 0c 


M e = Interface  Kinetics 
M c = Diffusion  Mobility 


Dynamics 


Coupled  CdhfrFfnfiard  & Calm-Alien  Equations 


Mtsr 


NIST 


Combined  Alloy  Free  Energy  Function  for  independent  a & (3  phases 


Wheeler,  Boettinger.  McFadden: 


f(<t>,C,T)  = (1- p((f)  ))fa(c,T)+  p(<t))fp(c,T) 


Alloy  Free  Energy  Function 

For  equilibrium  far  from  the  interface 

^ = 0 and  ^ = constant 

00 


Alternate  method: 

Kim.  Kim  & Suzuki,  Phys.  Rev  E 60(1999)  7186. 


mmp  common’ 

j "tgf  TAKCrST  PLANS 


NJissr 


LIQUID 


NIST 
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Interface  Energies  via  Thermodynamic  Functions;  natural  order  parameters 


/ ' ' r ' 


Dregia  & Wynblatt,  Acta  Met.  39(1991)774 


vsssss.wyAWMQMMmamtM» 


Diffusion  Mobility  Parameters  for 

Phase-Field  Calculations  for  Multicomponent  Alloys 


Ternary  Regular  Solution  Ag-Au-Cu  with  heats  of  mixing 
F = J {/  (c, , c2 , T ) + (Vc,  r + kI2  (Vc,  ■ Vc2 )+  *a  (Vc2  f )iV 


kn  =2  Zvrf2Q13 
k22=2Zd2£l2i 
kl2=2Zd2{Cln+Q.23-Q.u) 

d = interatomic  distance 
Zv=  interface  coordination  number 
= 3 for  { 1 1 1 } 

= 4 for(  100) 


Siinilarly,  gradient  energy  coefficients  Braun,  et  at.  Phil.  Trans,  r.  soc  Lon.  a U997i 
cart  be  obtained  for  order/disorder  cases!  wangetai.,ActaMat.<i998> 


• Phase  field  models  are  usually  derived  in  a volume  fixed  frame. 

• A particular  component,  say  ‘n\  is  picked  to  be  the  solvent. 

• One  needs  mobilities,  Ltj  "(<])),  for  the  dynamic  constitutive  law. 


xl’5F  Yi4a/(*n--*-> 


8X, 


ax, 


-1 


• Using  the  mobility  matrix,  M.  obtained  for  each  phase, 
the  L"  matrix  for  that  phase  is  given  by 


L’  = PMP 1 
where 
P =5  -X 


(matrix  multiplication) 


• Then  £."(<[>  )=p(<|))Z./'+(  l-p(<|>))L:" 


where  L,  and  L,  are  for  each  phase. 


Elastic  effects:  Cubic  Tetragonal  Transformation 


Electrochemistry 


3 field  variables, 

■■•'A ti 

♦ 


1 

1 

ti1  = 1 

tii  = 0 

T),  = 0 

TJ,  - 0 

t)2=0 

ri2=  1 

T}2=0 

r|2=0 

t|3=0 

r|3=0 

ti3=  1 

1)3=0 

t 

t 

t 

* 

Domain  1 Domain  2 Domain  3 


F(tp,Cj,yf)=  J 


••  "ss  w vsss.  ysssss* 


fv  (<P’C,)+  “ PW  + -f  IV<p[ 


dV 


F Helmholtz  free  energy 

fv  Helmholtz  free  energy  per  unit  volume 

<p  phase  field 

Cj  concentration  of  species  j 

\\>  electrostatic  potential 

P~  lLzjCF  charge  density 

N 

Ky  phase  field  gradient  energy  coefficient 


Stress  free  strain,  function  of  1);  add  elastic  energy  to  free  energy  functional 


Application  to  Electrochemistry,  Jon  Guyer  at  NIST 


Conclusions 

• Phase  field  calculations  are  based  on  thermodynamic  functions  and 
a few  simple  dynamic  postulates. 

• Calphad  evaluations  are  ideal  input  data  for  phase  field 
calculations. 

• For  thermodynamics  involving  natural  order  parameters,  gradient 
energy  coefficient  and  hence  interface  energies  can  be  calculated. 

• Otherwise,  gradient  energy  coefficients  and  interpolation  between 
individual  free  energies  must  be  obtained  by  fit  to  measured  surface 
energy 

• Diffusion  mobility  assessments  can  be  implemented  in  phase  field 
calculations. 

• Numerical  difficulties  still  limit  application. 
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Additional  points  made  by  speaker: 

Diffuse  interface  model  helps  to  overcome  the  moving  boundary  problem  of  the  sharp  interface  model. 
Phase  field  method  is  ideal  for  modeling  ordering  processes  since  it  naturally  implements  an  order 
parameter. 
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Current  and  Future  Applications  of  CALPHAD  Technology 
Larry  Kaufman,  Brookline,  MA 

Abstract 

During  the  past  year  the  author  has  had  the  opportunity  of  participating  in  two  international  symposia 
devoted  to  exploring  recent  examples  of  the  development  of  methods  for  predicting  multiphase 
equilibria  in  diverse  materials.  The  first  of  these,  held  at  the  130th  TMS  meeting  in  New  Orleans  was 
organized  by  Zi-Kui  Liu  and  entitled  “Computational  Thermodynamics  and  Materials  Design”  (1).  The 
second,  organized  by  Patrice  Turchi  was  entitled  “CALPHAD  and  Alloy  Thermodynamics”  was  held  at 
the  131st  TMS  meeting  held  in  Seattle  (2).  Most  of  the  papers  presented  at  these  symposia  are  published 
and  present  a graphic  record  of  the  significant  accomplishments  made  in  this  field  in  recent  years.  This 
workshop  provides  a fitting  epilogue  to  these  meetings  in  identifying  those  areas  where  development  of 
data  bases  which  are  needed  to  foster  future  progress  and  reporting  on  techniques  and  advances  in  such 
efforts.  The  author  has  chosen  several  examples  of  recent  work  in  the  analyses  of  multicomponent 
systems  used  in  lithium  battery  anodes  (3,4),  giant  magneto  resistance  (5),  zirconia  ceramics  (1,6), 
niobium  alloys  (7),  metal-carbon  systems  (1,8),  metallic  glasses  based  on  Al-Fe-Ni-Gd-Y  alloys  (9,10), 
transformation  kinetics  (1,11)  and  corrosion  (1 1)  to  illustrate  the  variety  of  applications  to  which 
databases  for  multicomponent  systems  can  be  applied  to  describe  equilibrium  and  kinetic  behavior. 
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Recent  Developments  in  Anode 
Materials  for  Lithium  Batteries 

M.M.  Thackeray,  J.T.  Vaughey,  and  L.M.L.  Fran&son 


Liihium-ion  batteries,  preferred  for 
their  high  energy  and  power,  also  present 
several  challenges.  Of  particular  concern 
are  unsafe  conditions  that  can  arise  in 
lithium-ion  cells  that  operate  with  a 
fully  lithiated  graphite  electrode.  If  the 
cells  in  those  batteries  are  overcharged, 
especially  in  large-scale  applications, 
thermal  runaway,  venting,  fire,  and 
explosion  could  result.  This  paper 
examines  research  into  alternative , 
interne talhc  electrode  materials. 

INTRODUCTION 

State-of-the-art  lithium  ion  batteries 
can  pose  a danger,  particularly  for 
large-scale  applications,  such  as  elecinc 
vehicles.  Safety  concerns  arise  from  the 
following  cell  chemistry  features: 

• At  the  top  of  charge,  the  lithiated 
graphite  electrode  operates  very 
close  to  the  potential  of  metallic 
lithium,  which  is  highly  reactive 

• ThedelhhiaiedLi^CoO.  electrode 
can  release  oxygen 

• The  electrolyte  usually  contains  a 
flammable  organic  solvent. 

Cells  that  are  overcharged  can  lead 
to  thermal  runaway,  venting,  fire,  and 
explosion.  For  these  reasons,  the  voltage 
of  every  lithium-ion  cell  in  a battery 
pack  has  to  be  controlled  by  electronic 
circuitry  to  prevent  the  cells  from 
becoming  overcharged.  Extensive 
research  efforts  are.  therefore,  under 
way  worldwide  to  find  alternative 
electrode  materials  and  non-flammable 
electrolytes  for  these  batteries 

Lithium-ion  bartencs  use  a carbon 
negative  electrode  and  a lithium  transi- 
tion metal-oxide  positive  electrode,  the 
prefened  materials  being  graphite  and 
LiCoOj.  respectively.  A fully  lithiated 
graphite  electrode.  LiC4.  provides  a 
specific  capacity  of  372  mAh/g  that 
translates  to  a volumetric  capacity  of 
818  mAh/ml  based  on  a density  of 
12  g/mL  Li/Tj/Li^CoO,  cells  operate 


typically  between  4.2  V and  3,5  V.  thus 
providing  high  energy  and  power. 

METAL  ALLOYS  AND 
INTERMETALUC  COMPOUNDS 

With  respect  to  negative  electrode 
(anode)  materials,  metal  alloys  or 
intermetallic  compounds  are  attractive 
alternatives  to  graphite  because  they 
can  be  selected  to  operate  between  0 and 
1 volt  above  the  potential  of  metallic 
lithium.  Many  intermetallic  compounds, 
particularly  the  binary  systems.  Li.M 
(e.g . M = A).  Si.  Sn.  Sb.  Pb)  have 
already  been  thoroughly  investigated  .' 
Unfortunately,  metals  have  very  dense 
structures,  so  reactions  with  lithium  tend 
to  be  accompanied  by  large  changes 
in  volume  and  by  major  structural 
rearrangements.  The  severe  crystal- 
lographic changes  that  occur  during 
charge  and  discharge  of  the  cells  lead 
to  the  mechanical  disintegration  of  the 
intermetallic  structures,  insulating  layers 
on  the  surface  of  tbe  electrochenucally 
pul  verized  electrode  particles,  and  a loss 
of  electronic  particle-to-partkle  contact 
These  effects  seriously  compromise 
the  cycling  efficiency  and  cycle  life 
of  lithium  cells.  A good  example  is 
aluminum,  which  reacts  with  lithium 
at  approximately  300  mV  vs.  lithium 
according  to  the  reaction: 

Li  + Al  -*  LiAl  (I) 

Both  aluminum  and  LiAl  have  face- 


centered- cubic  (fee  I structures  ( Figu 
However,  in  Li  Al.  the  aluminum  a 
no  longer  assume  the  same  posi 
as  they  do  in  the  parent  alumi 
structure,  the  phase  transition  req 
the  diffusion  of  both  aluminum 
lithium  atoms,  which  leads  to  a stru 
that  may  be  defined  in  terms  of 
independent  fee  aluminum  lattice- 
intersect  with  two  independen' 
lithium  lattices  The  volume  expa 
that  accompanies  the  transitu 
nearly  200%  per  aluminum  3 
It  is  not  surprising,  therefore, 
aluminum  electrodes  do  not  ma: 
their  structural  integrity  onelectrochei 
cal  cycling 

Considerable  effort  has  been  madt 
overcome  the  limitations  of  imermeu 
electrodes  by  using  composite  man 
als,  in  which  an  electrochemica 
active  phase  is  embedded  within 
eleciroc heroically  inert  phase  "*Awi 
known  example  is  SnO,1  which  rea 
with  lithium  m two  distinct  stages: 

4 Li  + SnO  -»  Sn  + 2 LuO 
x Li  + Sn  *-*  Li.Sn  (x  < 4.4) 
The  first  reaction,  in  which  SnO 
reduced  to  metallic  tin,  is  irreversi 
because  it  is  difficult  to  remove 
lithium  clecirochemically  from 
insulating  Li.O  matrix.'  The  sect 
reaction,  which  can  involve  the  ton 
tion  of  several  phases  in  the  Li 


stem,  is  reversible.  Although  the 
pacity  retention  of  SnO  electrodes  is 
;mficantly  superior  to  a pure  Li^Sn 
item,  the  composite  electrodes  display 
arge,  irreversible  capacity  loss  on  the 
st  cycle,  which  has  precluded  their 
plication  in  lithium-ion  cells.  Olher 
mposite  electrodes,  such  as  FeSn2,  in 
iich  Li.Sn  operates  within  an  inactive 
« (or  FeSnC)  matrix,  have  been  found 
show  similar  limitations  * 

Graphite  and  LiCoO;  electrodes  show 
od  cycling  behavior  because  they 
erate  as  insertion  electrodes  (i.e., 
hium  can  be  inserted  in  and  out  of  their 
uctures  without  significant  structural 
mage  or  decay).  With  this  in  mind, 
gonne  National  Laboratory  recently 
gan  to  seek  intermetallic  structures 
it  could  accommodate  lithium  with 
nimum  structural  change,  that  is.  the 
rent  and  lithiated  structures  should 
ow  a strong  structural  relationship 
ward  one  another  Several  com- 
-unds  with  interesting  structure- 
xtrocbemical  property  relationships, 
ve  been  identified.  Those  compounds 
elude  Cu^Sn,.’  InSb.*"1*  Cu,Sb,” 
njSb.  and  MnSb  15 


INSERTION  ELECTRODE 
OPTIONS 

Cu,Sns 

Cu4Sn,  (Cu^Sn)’14  has  a nickel - 
arsenide-type  structure  in  which  20% 
of  the  interstitial  sites  arc  occupied  by 
copper;  it  reacts  with  lithium  to  form  a 
stable  lithiated  zinc-blende  type  structure 
LijCuSn  according  to  the  reaction 
10  Li  + Cu0Sn,  — * 5 Li,CuSn  + Cu(4) 
The  ideal  nickel-arsenide  to  zinc -blende 
phase  transition,  shown  in  Figure  2. 
may  be  regarded  as  topo tactic  because 
the  CuSn  zinc -blende  framework  can  be 


obtained  if  50%  of  the  un  atoms  in  the 
nickel  arsenide  (CuSn)  component  of 
the  parent  Cu4Sn,  structure  are  displaced 
into  neighboring  interstitial  sites.  The 
displacement  of  the  relatively  large  tin 
atoms  during  the  transition  limits  the 
kinetics  of  the  reaction  and  leads  to 
hysteresis,  particularly  during  charge 
The  copper  atoms  that  reside  in  the 
interstitial  sites  of  the  nickel-arsenide 
framework  are  extruded  from  the 
structure  dunng  Reaction  4 The  oickcl- 
arsemde  to  zinc-blende  phase  transition 
is  accompanied  by  a 59%  expansion 
in  crystallographic  volume.  Further 
lithiation  of  Li2CuSn  displaces  the 
remaining  copper  from  the  structure 
to  yield  Li44Sn 

2.4  Li  ■*-  LijCuSn  — » Li„So  -*■  Cu  (5) 
Cu4Snt  electrodes  have  a rechargeable 
capacity  of  approximately  200  mAh/g 
when  cycled  between  1.2  V and  200  V 
vs.  lithium  (i.e.,  if  the  reaction  is 
restricted  essentially  to  Reaction  4) 
Although  this  does  not  compete  with  the 
specific  capacity  of  graphite  electrodes 
(theoretical  value  = 372  mAh/g).  the 
volumetric  capacity  that  is  obtained 
experimentally  from  Cu6Sn,  dunng 
Reaction  4.  that  is.  1.360  mAh/ml 
(based  on  an  average  electrode  density 
of  6.8  g/ml)  is  significantly  superior  to 
the  theoretical  volumetric  capacity  of 
graphite  (818  mAh/ml). 

InSb 

The  recognition  that  a zinc-blende 
(diamond-like)  framework  such  as 
CuSn  in  Li2CuSn  provides  a cubic, 
three  dimensional  interstitial  space 
for  lithium  immediately  prompted  the 
investigation  of  zinc -blende  structures 
such  as  AlSb.  GaSb.  and  InSb  as 
possible  electrode  hosts  for  lithium.*-11 
The  performance  of  InSb  electrodes  is 
particularly  noteworthy  The  overall 
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Stable  and  metastable  evaluation  of  the  Cu-Co-Fe  Phase  diagrams 
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M Bamberger.  A Munitz*.  L Kaufman**  . and  R Abbaschian"** . 

Materials  Saence  and  Engineering.  Technion.  Techmon  City.  Haifa.  ISRAEL 

•Nuclear  Research  Center-Negev  P.O.  Box  9001.  Beer-Sheva,  ISRAEL  5848 
^140  Clark  Road.  Brookline.  MA  02445-5848.  USA 

""Department  of  Materials  Science  and  Engineering,  University  of  Florida.  Gainesville, 
Flonda,  32611  USA 

The  Thermo-Calc  SSOL  database  for  the  Cu-Co-Fe  ternary  system  was  modified 
according  to  expenmental  data  obtained  by  electromagnetic  levitation.  Using  the 
modified  database,  the  solidus  and  the  liquidus  temperatures  were  calculated  Very 
good  agreement  was  obtained  between  the  calculated  and  the  expenmental  data.  The 
metastable  two-melt  separation  and  the  compositions  of  the  two  melts  were  calculated 
as  well,  and  compared  with  the  expenmental  results. 

1)  INTRODUCTION 

Recently,  interest  has  been  growing  in  extension  of  solid  solubility  of  alloys 
exhibiting  metastable  miscibility  gaps,  such  as  Cu-Co  alloys.  Such  materials  are  known 
to  exhibit  giant  magneto  resistance  (GMR)'1,Z).  namely  a large  drop  in  the  electrical 
resistance  under  the  influence  of  a magnetic  field.  For  example,  a supersaturated 
CusoCoio  solid  solution  annealed  at  440°C  for  different  durations  exhibits  an  increased 
GMR  of  up  to  1 1 % at  room  temperature.  It  is  believed  that  the  heat  treatment  caused 
solid-state  spinodal  decomposion.  which  is  responsible  for  the  enhanced  GMR 
Previous  studies  showed  that  supercooling  of  Cu-Fe.  Cu-Co  or  Cu-Fe-Co  alloys  beyond 
a certain  limit  results  in  metastable  separation  of  the  melt  into  two  liquids,  one  Cu  rich 
(called  L2),  the  other  Co  or  Fe  rich  (called  LI)*3-6*.  Liquid  phase  separation  and 
formation  of  metastable  phases  were  observed  during  rapid  cooling  of  these  alloys 
However,  there  is  only  meager  information  on  the  stable  ternary  Cu-Fe-Co  diagram 

Recently,  some  experimental  data  for  the  melting  temperatures  were  published 
concerning  the  Cu-Co-Fe  system'11’.  Based  on  the  Fe-Cu(7,B>  and  Co-Cu(88)  binary 
phase  diagrams,  approximate  isothermal  sections  have  been  constructed'10’. 

The  present  work  is  aimed  at  obtaining  a reliable  database  for  the  Cu-Co-Fe 


system,  and  to  perform  consistent  calculations  of  the  stable  phase  diagram. 
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Thermodynamic  Modeling  of  the  Zr-0  System 


Raymundo  Arroyave,  Larry  Kaufman  , and  Thomas  W Eagar 

Department  of  Materials  Science  and  Engineenng.Room  4-134 
Massachusetts  Institute  of  Technology,  Cambridge,  MA  02139  USA 
rayroundo@mit  edu, Larry kaufman@rcn.com, tweagarOmit.edu 


Abstract.  In  this  study,  the  complete  zirconium-oxygen  system  has  been  critically  assessed  at  1 at  from 
300°C  to  liquidus  temperatures.  Thermocbemical  measurements  and  phase  diagram  data  were  used  to 
model  the  Gibbs  free  energies  of  seven  phases  Additionally,  the  ordered  interstitial  HCP-based  solutions 
were  included  and  considered  as  simple  line  compounds  By  using  the  PARROT  module  in  Thermo-Calc, 
it  was  possible  to  optimize  the  parameters  of  the  models  used  to  describe  the  Gibbs  free  energies  of  the 
HCP.  BCC.  liquid,  t ZrOi-z.0  ZrOj-z  and  a ZrOj_£  phases.  The  Gas  phase  was  considered  to  behave 
ideally  Although  phase  diagrams  including  the  stoichiometric  zircorua  phases  have  beer,  assessed , this  is 
the  first  time,  to  the  best  of  our  Itnowledge  that  a complete  assessment  of  this  system  is  published 


SCC  BCC/HCpTm*«Cp'  ,|'  11 
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Figure  2.  Calculated  Isothermal  Sections  in  Nb-Ti-Crat  1800C, 
1600C.1500C  and  1300C.Two-Phase  Regions  are 
identified  by  Tie-Lines,Three-Phase  Regions  by 
Triangles  and  Single  Phase  Regions  are  labled. 
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Figure  3. Calculated  Binary  Phase  Diagram  Components  of 
the  Ti-Al-Cr  System.The  Liquid  Phase  is 
designated  by  L 


Figure  4.  Calculated  Isothermal  Sections  in  Ti-Al-Cr  at  1 500C, 
1300C,1200C  and  1 OOOC.Two-Phase  Regions  are 
identified  by  Tie-Lines,Three-Phase  Regions  by 
Triangles  and  Single  Phase  Regions  are  labled. 


Figure  13. Calculated  Isopleth  constrained  by  the 
Conditions  that  weight%Cr=4wt%Al 
and  wt%Ti=3wt%Cr. 


Figure  17.Calculated  Isothermal  Section  in  the  Nb-Cr-Mo 
System  at  1200C. 
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Validation  of  Thermo-Calc  Database  with  Experiment 

Experimental  versus  Calculated  T eman/  Phase  Diagrams  (Ni-Cr-Mo) 


Equilibrium  and  Scheil  Results  tor  Solidification 
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Kinetics  ot  Alloy  Phase  Transformation:  TTT  Diagrams 


Pretsminary  results  Irom  OICTRA:  application  to  temporal  evolution  of  ordering  rn  Ni-Cr 
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99  SETUP  FILE  FOB  CALCULATING  A Ni2Cr /AUSTENITE (fee) 

60  TRANSFORMATION  IN  A BINARY  Ni-Cr  ALLOY. 
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90  THE  INITIAL  STATE  IS  AN  AUSTENITE  OF  lOOmicron  THICHNESS 
03  THE  COMPOSITION  OF  THE  AUSTENITE  IS  Ni-33at%Cr. 

99  APTER  AUSTENITISATION  THE  SPECIMEN  HAS  BEEN  QUENCHED  DOWN 
09  TO  273K.  THE  SYSTEM  IS  ASSUMED  CLOSED,  SO  WE  NEED  NOT  SET  ANY 
00  BOUNDARY  CONDITIONS  (CLOSED  SYSTEM  IS  DEFAULT) . 

09 

00  Ni2Cr  IS  EXPECTED  TO  GROW  INTO  THE  AUSTENITE.  FOR  THIS  REASON 
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Fig. 14  Experimental  measurment  of  TCP  phase  formation 
in  C-22  by  Tammy  Summers  and  Sharon  G.  Torres 
compared  with  calculation  of  P_Phase  formation  in 
a Ni-21.1-13.5Cr  Alloy. 


Figure.  16  Examination  of  the  stability  and  phase  diagram  of  the  fee  and  bee  phases  in  the  Cr-Ni 
System.Top  row(95Cha),Middle  row  SSOL,Botiom  row(OlKau). 


DISCUSSION  OF  THE  STABILITY  OF  THE  FCC  AND  BCC  PHASES  IN  THE  Cr-Ni  SYSTEM 

As  indicated  in  the  INTRODUCTION  the  LS  values  used  in  current  CT  work  have  been  criticized 
Because  ab  initio  values  for  some  elements  are  substantially  larger  in  magnitude  that  the  LS  values 
that  are  employed.Imtially^Kos^SCha  and  98Gri)  it  was  suggested  that  these  differences 


Isoatructural  Enthalpies  of  Formation  for 
tha  Ni-Cr  Syatam 


Cr  Ni 

XNi 

Calculated  using  tha  phenomenological  (solid  lines) 
and  the  first-principles  TB-CPA  (dashed  lines)  models. 
Curves  (A),  (B),  and  (C)  correspond,  respectively,  to 
Skriver's  (36.66  kJ/mol),  an  intermediate  (23.01 
kJ/mol),  and  Kaufman's  (10.46  kJ/mol)  values  of 
„ fcc/bcc 
H Cr 
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Additional  points  made  by  speaker: 

For  the  assessment  of  the  thermodynamic  quantities  of  the  TCP  phases  in  Ni-X  systems  results  from  ab 
initio  calculations  were  taken  into  account. 

The  thermodynamic  description  of  the  Ni-Cr-Mo  system  was  used  to  generate  Pourbaix  diagrams  for  the 
prediction  of  the  corrosion  behavior  under  various  conditions. 


. 
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Thoughts  on  the  Acceleration  of  Phase  Equilibria  Research 
J.-C.  Zhao,  General  Electric  Company,  Schenectady,  NY 


Acceleration  - Thoughts  for  Discussion 
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Acceleration  - Thoughts  for  Discussion 


Diffusion  Multiple  Approach 
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Acceleration  - Thoughts  for  Discussion 

Diffusion  Multiple  Approach 


Acceleration  - Thoughts  for  Discussion 
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Acceleration  - Thoughts  for  Discussion 

Example:  Nimonic  263  Testing  of  ihermnsiyoormc  <jamoa» 


___  Acceleration  - Thoughts  for  Discussion 
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Acceleration  - Thoughts  few  Discussion 

<if  iiwjxmipmivio 

Improvement  Needs  for  Ni-data 
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Acceler  ation  - Thoughts  for  Discussion 
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Additional  points  made  by  speaker: 

Predictions  of  TCP  phases  for  Ni-Base  systems  with  currently  available  thermodynamic  databases  are 
not  very  reliable  =>  lack  of  experimental  data  needs  to  be  overcome  to  facilitate  more  accurate 
modeling  of  these  phases.  Accurate  ternary  phase  diagram  data  are  needed  for  the  improvement  of  the 
model  descriptions  of  these  phases  for  multicomponent  databases. 


